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The important thing is not to 
stop questioning. Curiosity has 
its own reason for existing. 

Albert Einstein 



Et voila ... est fait ! 



Vous avez ete nombreux, ces demieres annees, a vous poser la meme question. Certains d'entre vous 
allerent meme jusqu 'd la prononcer d haute et intelligible voix. Parfois avec malice, mais le plus 
souvent avec compassion. U inquietude etait meme perceptible dans la voix de certains. Pour etre 
franc, je me la suis souvent pose moi-meme, cette question ... 

Et alors cette these, ga avance ? 

Je sais, je n'ai que trap rarement repondu a vos interrogations. La raison en etait tres simple: la 
seule chose queje savais, c'est qu'ilfaudrait de toute fagon qu'elle soit terminee a temps. Mais voila, 
le temps passe et avec lui, les annees defilent. Le travail accompli prend forme. La demiere annee 
de these pointe le bout de son nez avec, dans ses bagages, la redaction tant redoutee. L'hiver est 
rude, le printemps accueilli d bras ouvert. Mais celui-ci cede dejd sa place a I'ete qui lui-meme se 
termine. Cette these qui a suscite tant d' interrogations prend forme. Petit a petit, les corrections et 
modifications defilent. Jusqu' a ce que ... Aujourd'hui, je suis enfin pret a repondre a cette question. 
Alors, n 'hesitez pas ! 

Vous avez ete nombreux a contribuer d'une maniere ou d'une autre a V elaboration de cette these. 
Certains directement, d'autres par leur seul presence. II ne me reste plus qu'd trouver les mots pour 
vous remercier tous. 

La premiere fois que j'ai franchi la porte du bureau de Peter, je cherchais un sujet de memoire. 
J'ignorais a I'epoque qu'entrer dans ce bureau me conduirait a terminer une these cinq ans plus 
tard. Mais a posteriori, si tout etait a refaire, je franchirais a nouveau cette porte! Durant toutes 
ces annees, Peter a toujours ete disponible, que ce soit pour repondre a mes questions toujours plus 
nombreuses (et pas forcement bienformulees, voire meme parfois denuees de sens!) ou pour discuter 
de nouvelles idees d explorer Meme s'ilm'est souvent arrive d'etre encore plus perplexe apres avoir 
discute avec lui, je pense avoir beaucoup plus appris aupres de lui que durant n 'importe quelle autre 
formation que j'aipu suivre. Je lui exprime des lors toute ma gratitude. 

Comme toute these, celle-ci a ete realisee au sein d'un groupe de recherches, le Service de 
Physique Theorique de I'Universite Libre de Bruxelles. Que de bons moments passes durant ces 
quatre annees a discuter avec ses membres (durant les breaks bien sur) ou a f aire les cent pas dans 
cefameux couloir du N07. Merci a vous tous, Chiara, Frangois-Xavier, Gilles, Hylke, Isabelle, Jean- 
Marie, Jonathan, Manu, Paola et Thomas. Je souhaiterais tout particulierement remercier Michel 
pour avoir ete mon promoteur « offtciel » pendant trois ans, avec tout le travail administratif que cette 
charge implique, et Laura pour ses conseils qui tomberent toujours a point ainsi que pour sa patience 
sur les pistes de ski ( et oui, le travail c 'est dur!). Merci egalement a Josep dont j 'ai partage le bureau 
pendant pres d'un an et qui m'a sorti a plusieurs reprises de mes problemes sous Linux. Merci a 
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vous, Quentin et Sabrina, pour tous cesfous rires qu 'on a portages dans ce meme bureau et qui m 'ont 
permis de m' evader un instant. Merci Fu-Sin pour ces quelques breaks au Soleil quifurent trop peu 
nombreux mais qui, j'en suis sur, seront plus frequents a I'avenir. 

Ces huit annees passees a I'Universite Libre de Bruxelles furent une experience unique. Je 

souhaiterais des lors remercier V ensemble du corps professorcd du departement de physique pour 
I'enseignement regu ainsi que les etudiants qui m'ont accompagne dans ce qui s'apparentait parfois 
a un veritable parcours du combattant. Merci tout particulierement a Nassiba et Vincent. 

Cette these, je la dois en grande partie aux professeurs du College Don Bosco de Woluwe-Saint- 
Lambert qui m'ont enseigne des disciplines aussi varices que I'histoire et les mathematiques. Je 
voudrais remercier en particulier pour leur passion communicative Albert Demelenne, Ingrid t'Kint, 
Yves Monin, Philippe Moreau, Veronique Bouquelle ainsi que Frangoise Damien a qui je dois mes 
premieres observations de M 31, Jupiter et Satume. 

J'adresse egalement mes remerciements a tous les membres de ma famille pour m' avoir tou- 
jours soutenu a chaque etape de mon existence. Je remercie en particulier ma « petite » soeur pour 
I' illustration placee sur la couverture de cette these. 

Je ne sais pas de quoi demain sera fait. Mais je sais d'ou je viens et a quije dois d'etre ce que 
je suis. Pour vos valeurs que vous m 'avez enseignees, pour votre patience, pour m 'avoir toujours 
soutenu, pour tout cet amour que vous m 'avez prodigue et pour tous ces moments de joie que vous 
m 'avez offerts. Papa, Maman, merci. 

Ces remerciements ne seraient pas complets si j'oubliais celle qui supporte mes humeurs et mes 
interrogations depuis si longtemps. Merci pour ton amour et ton soutien sans failles. 



Michael, 
Septembre 2009. 
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Abbrtviaiions & Conventions 



Throughout this thesis, we will work in "natural unites" for which the speed of light and the Planck 
constant are one, i.e. c = /i = 1. Greek indices will be used for space-time coordinates (/x = 0, . . . , 3) 

and Latin indices for space coordinates {i = 1,2,3). The metric "sign" convention will be the one 
used by Landau & Lifchitz in The classical theory of field, for which the space-time interval is given 
by 

ds^ = di^ - d/2 + . . . . 

Therefore, the metric of Minkowski describing flat space-time is given in Cartesian coordinates by 

7?^,, = diag. (1,-1,-1,-1). 
Finally, the following abbreviations and conventions will be used throughout this thesis : 



CMB 

EoM 

FLRW 

FP 

GR 



G = 

9 

k = 
Ac 

n 



Cosmic Microwave Background 
Equation(s) of motion 
Friedmann-Lemaitre-Robertson- Walker 
Fierz-Pauli (theory of) 
General Relativity 
Newton's potential 

Newton's constant 
Curved metric 
Determinant of the metric 
Einstein tensor 

Three-dimensional momentum 

Cosmological constant 

Planck Mass 
Ricci scalar 
Ricci tensor 
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Motivations & outline 



When General Relativity (GR) was formulated in the beginning of the twentieth century, the anoma- 
lous perihelion advance of Mercury was an unsolved problem. When it was clear that Einstein's new 
theory predicted the same perihelion advance as the one observed, GR became the standard theory for 
the gravitational interaction. Its popularity, firstly motivated by its mathematical elegance and consis- 
tence, was reinforced by its prediction of the bending of light by a massive source, which was first 
observed in 1919. Since then, GR has been tested with an ever-growing precision [1, 2, 3, 4, 5, 6, 7]. 
However, gravity remains little-known mainly because of its particularly low coupling to matter'. 

GR predicts [8, 9] the existence of gravitational waves, in analogy with Maxwell's theory of 
electrodynamics. There are many similarities between gravity and electromagnetism, but also striking 
differences. For instance, while the existence of photon is well established, gravitational waves have 
still not been directly observed despite the huge efforts invested recently in developing gravitational 
waves detection experiments [10, 11]. At present, the best evidence one has for the existence of 
such waves comes from timing observations of binary pulsars [12, 13]; the secular decrease of orbital 
period of the system PSR 1913-1-16 has been shown to be compatible with the emission of gravitational 
waves predicted by GR [14, 15]. Hence, gravitational waves are in a situation similar to neutrinos 
after Pauli's proposal in 1930, when their emission could be inferred from energy and momentum 
conservation in beta decays without being observed in absorption experiments. 

According to GR, gravitational waves have zero rest mass, just as electrodynamic waves. If gravi- 
tons were to have a small, non-zero mass m, one would expect that all effects arising because of their 
mass would be small corrections to Einstein's theory of gravity, completely negligible at distances 
smaller than the graviton's Compton wavelength A ~ mT^ . Moreover, it is well known in fields 
theory that the very concept of a field rest mass is closely related to the range of the interaction it 
mediates. If the gravitons were to have a non-zero mass, one would expect the potential of a static 
source to have a Yukawa shape e^™''/r effectively cutting off the gravitational interaction at distances 
larger than mT^ . For this reason, much attention has been paid to the apparent absence of such cutoff 
[16, 17, 18]. The dynamic of gravitational waves allows for independent constraints on the gravi- 
ton rest mass through the modification of their dispersion relation, which implies that these waves 
propagate at velocity 



This relation implies that massive gravitational waves of different frequence should propagate at dif- 
ferent speed. Direct observations of gravitational waves could then be used to set a bound on the 
graviton mass, for instance from observations of compact binary systems which slowly spiral together 

'The gravitational interaction between two protons is 38 orders of magnitude smaller than the Coulomb interaction 
between them. 
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due to the emission of gravitational waves [19, 20, 21, 22, 23, 24]. Other bounds on the graviton rest 
mass could be set through the delay between electromagnetic and gravitational waves [25] or through 
the relativistic time delay [26] (for a recent review on graviton mass limits, see [27]). 

Although observations made so far seem compatible with massless gravitational waves, a very 
small graviton rest mass cannot be excluded. Hence, one may wonder whether it is possible to modify 
Einstein's theory of gravity as to describe massive gravitons. Giving a mass to gravitons is a chal- 
lenging problem [28], which has motivated many studies since the original work of W. Pauli and M. 
Fierz [29]. Because of the intricate structure of GR, it is much more complicated to give a mass to 
gravitons than imitating Proca's work on massive photons, with consequence that no clear consensus 
has emerged on the subject. 

Beside the theoretical interest of such question, theories with massive gravitons also have more 
pragmatic motivations. Surprisingly, one of them comes from Quantum Chromodynamics (QCD) 
[30]. Indeed, it is believed that QCD can be approximated by some sort of string theory. However, 
string theories often have massless spin-2 modes, while QCD does not. This apparent contradiction 
may be an illusion if it turns out that these modes can be removed from the massless sector of string 
theory by giving them a mass, a task which may seem very similar to giving a mass to gravitons. 
Yet, there are other original motivations for theories with massive gravitons. For instance, massless 
scalar fields only interacting with gravity could exist in principle. If so, they could be responsible for a 
non-zero graviton rest mass, just because of their interaction with gravity. The class of models which 
are the subject of this thesis makes use of this possibility. 

All these considerations legitimate the efforts invested in the attempt of building a consistent the- 
ory with massive gravitons. Yet, the main motivation for such theories is elsewhere. Generalizations 
of Einstein's theory of gravity with a small, non-zero, graviton mass may give rise to large scale (or 
infrared) modifications of GR. Despite that Einstein's theory has been directly tested from scales of 
a fraction of millimeters up to Solar System scales, there is no direct confirmation of it on larger 
distances, with consequence that infrared deviations from GR cannot be excluded. The recent de- 
velopments in theories describing massive gravitons were parts of this revival of interest for infrared 
modification of GR. 

Theories that predict infrared modifications of the gravitational interaction have two main moti- 
vations. The first one is to be found in theories with brane-worlds and extra-dimensions of large or 
infinite size (see [31] for a review). In these models, the usual matter is supposed to reside on a 3- 
brane embedded in higher-dimensional space, while gravity may propagate in all dimensions [32, 33]. 
Hence, they generally predict that gravity is modified at very large scales while being conventional 
at smaller ones. A notable example is given by the Dvali-Gabadadze-Porrati (DGP) model [32], for 
which gravitons may appear as a massive resonance [34, 35] to observers trapped in the 3-brane. The 
second motivation for large scale modifications of GR is to be found in cosmological observations. 

The dark sides of the Universe 

Because of the non-renormalizability of Einstein's theory of gravity, it is commonly believed that 
GR gives way to a quantum theory of gravity [36, 37, 38, 39] at high enough energies. For dimen- 
sional reasons, it seems natural to identify this energy scale with the Planck mass Mpj = (SvrG)"^, 
where G is Newton's constant. The standard assumption is that GR is valid as an effective field the- 
ory- for energies below Mpi, and therefore for length scales larger than Api = Mpj^. Nevertheless, 
at very low energies GR also faces some troubles in explaining observations made from galactic up 

^For reviews on effective field theory, see [36, 37, 38, 39]. 



3 



to cosmological scales; the agreement is only achieved after the introduction of the otherwise unde- 
tected dark matter and dark energy. However, with some assumptions about these dark components, 
a picture emerges which describes with an impressive precision the whole bulk of cosmological data 
[40, 41, 42, 43]. The question arises to what extent one should consider this agreement as a confirma- 
tion of GR itself. It is not inconceivable that large scale modifications of GR, such as theories with 
a small graviton mass, may account for these observations without need of dark hypothesis. Over 
the last years, there have been several attempts [44, 45, 46, 47, 48, 49] in constructing self-consistent 
theories which deviate from GR at large distances and time scale\ most of them being motivated by 
these "dark paradigms". 

The existence of dark matter in the form of an undetected planet was already assumed in the 
19th century to reconcile astronomical observations and Newton's theory of gravity; the astronomers 
Jean Joseph le Verrier and John Couch Adams both predicted the existence of a new, unseen planet by 
studying the anomalies in the orbit of Uranus. Although Neptune was discovered in 1 846 at the precise 
location predicted by these astronomers, this "dark matter"problem also had an alternative explanation 
as modification of gravity. Soon after, le Verrier started a complete study of orbital motions of the other 
planets. Motivated by his previous achievement, he postulated the existence of another undetected 
planet, Vulcan, between Mercury and the Sun to explain the periheUon advance of Mercury. However, 
this time the solution was not a missing dark matter, and the problems in the motion of Mercury were 
solved by Einstein's new theory of gravity. Nowadays, both dark matter and dark energy are needed 
to reconcile theory and observations. 

Modem dark matter was first advocated by Jan Henrik Oort [50] and Fritz Zwicky [51] in the 
early 1930s to explain vertical motions of stars in the Milky Way and radial velocities of galaxies 
in the Coma cluster, respectively. Since then, the number of observational evidences [52] for non- 
baryonic dark matter has burst out (see [53, 54, 55] for reviews): dark matter is needed from galaxy 
rotation curves [56, 57] up to structure formation [58]. Moreover, there were recently claims for direct 
evidences for the existence of dark matter from two colliding clusters of galaxies known as the Bullet 
Cluster [59]. 

The problem of modern dark matter is very similar to the old problem of unseen planets. One may 
either resolve it by assuming the existence of a large amount of unseen, dark, matter existing on system 
with size ranging from galactic to cosmological scale, or by assuming deviations from the standard 
laws of gravity at large distances and time scales [60]. Both approaches are popular. An example 
of the first approach is given by weakly interacting massive particles (or WIMP's) [54]. The idea is 
to suppose the existence of a stable, weakly interacting, massive particles which has been produced 
in a huge amount in the early Universe and which could explain the discrepancy between theory and 
observations through its interaction with gravity. The second approach, often illustrated by Modified 
Newtonian Dynamics (MOND) [61] and theories attempting to incorporate it in a full relativistic 
description [62], consist in searching for possible alternatives to GR which are often designed to 
explain rotation curves of galaxies without dark matter. 

Modern physicists face more serious troubles when trying to understand the nature of dark energy 
(for reviews on the subject, see [63, 64]). First evidence for dark energy dates back to the early 
1990s. At that time, it become clear from the analyze of the Cosmic Microwave Background (CMB) 
power spectrum that the Universe is nearly flat [65, 66]. Yet, the fractional energy density of all 
non-relativistic matter (baryons and dark matter) was bounded by r^Matter < 0.3 [67, 68] so that 
something was missing in order for the Universe to be flat. Nowadays, the existence of dark energy 

^The general philosophy of those works is to view these models as the low energy limits of an unknown fundamental 
theory, without worrying about issues like renormalizability or embedding into an ultraviolet-complete theory. 
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is inferred from many cosmological observations [69, 70] all indicating that the Universe is currently 
accelerating. In the framework of GR, this is only possible if the Universe is dominated by a fluid 
with negative effective pressure 

VK _ ^1 

— = U-A < --, 

PA 3 

called dark energy and which is some kind of effective version of Einstein's cosmological constant. 
Dark energy is assumed to be responsible for roughly 70% of the total energy content of the Universe, 
= 0.742 ± 0.030 [Al>f, with an energy density given by 

PA = ?,HlMl^^ ~ 4 . 10-6 ^ ~ 10-29^. 

This energy scale is much smaller than the energy scale one would associate, on dimensional grounds, 
to the three known interactions 

PA ~ lO-^^PQCD ~ lO-^'^pEW ~ lO^^^'^Pgravity 

Its seems then really difficult to understand the origin of dark energy. The unnatural smallness of p\ 
is a huge problem. It implies that the vacuum energy density-^ is canceled in the gravitational field 
equations to a very small value by a mechanism which does not seem to be related to any of the three 
known interactions. In fact, this problem could be divided into two parts. The first part consists in 
understanding why is dark energy not as high as it should be to agree with estimations of the vacuum 
energy density [63]; in other words, why is pA essentially zero. The second part of the problem consist 
in understanding why pA is in fact different from zero. 

Although it is not inconceivable that the first part of this problem may be solved by a mechanism 
which drives the cosmological constant to zero in the early Universe [71], its second part requires to 
understand what is hidden behind p\. There are various hypotheses as for the nature of dark energy 
[72, 73]: 

• In order to fit the cosmological observations, assumptions are made about our Universe. In 
particular, the Universe is assumed to be homogeneous and isotropic, although this is clearly 
not the case below scales of ~ 100 — 200 Mpc. The acceleration of the Universe's expansion 
could be an illusion if one actually hves in a kind of void, where the matter density is smaller 
than the spatial average density [74, 75, 76] (see [77] for a review). 

• Models with weakly interacting, time-varying, negative pressure fluids such as quintessence 
[78, 79, 80] are perhaps the most popular candidates for solving the dark energy problem (see 
[81] for an introduction to quintessence) since such fields naturally imply an acceleration of the 
Universe's expansion at late time. 

• The weak anthropic principle [82] could also somehow explain the observed value of the vac- 
uum energy density [63, 83]. Indeed, if the Universe is much larger than its visible part, and if 
the cosmological constant take different values in different regions of cosmological size, then 
we happen to live in a region where the effective vacuum energy density is small enough to 
allow galaxies and stars to form. In other regions where the effective cosmological constant is 
not that small, there is simply nobody to worry about this issue. 



Standard baryonic matter and cold dark matter amount only to Slflaryons ~ 0.04 and Hcdm = f^Maitcr — f^Baryons ~ 0.22, 
respectively. 

^Gravity couples universally to all forms of energy, including the vacuum energy which acts like a cosmological constant 
in Einstein's equations. 



5 



This list is not exhaustive. Independently of the relevance of all these hypothesis, large scale deviation 
from GR may also account for the unnatural smallness of pA [84, 85]. For instance, theories based on 
an arbitrary function of the Ricci scalar [86] are designed to fit the cosmological observations without 
need of dark energy. But what is even more interesting is that theories with massive gravitons may 
also contribute to the cosmological constant and therefore shed light on the dark energy paradigm (see 
chapter 6). 

Massive gravitons 

It is an open question whether all current observations could naturally fit into a unique theory of 
gravity without need of dark matter and, or, dark energy. In order to address this question an alternative 
model is needed whose predictions can be compared to those of GR. Theories with massive gravitons 
are candidates for such alternatives. They are interesting from both theoretical and phenomenological 
point of view, independently of the answer to this question. 

As already mentioned, the first attempts trying to describe massive gravitons were based on the 
original work of W. Pauli and M. Fierz [29]. These attempts were concerning a Lorentz-invariant 
model, in the linear regime, which was supposed to generalize Einstein theory of gravity with massive 
gravitons. This model, known as the Fierz-Pauli theory, is the only one which is both Lorentz-invariant 
and free of instabilities in Minkowski space-time. Unfortunately, it turned out that it could hardly be 
considered as a coherent candidate for a model of massive gravity (i.e., with massive gravitons). 
Although some of its problems may be solved by a non-linear completion, any theory with Lorentz- 
invariant massive gravitons of the Fierz-Pauli type suffers from strong coupling at very low energy 
scales [87]. Therefore, such theory requires an UV completion to describe the gravitational interaction 
at scales where GR has already been tested. 

It is not clear how to build consistently a theory of massive gravity. Yet, the example given by the 
Fierz-Pauli model enables one to understand that it is essential to have a full non-linear formulation 
for both theoretical and observational considerations. Because of the success of Einstein's theory of 
gravity, it seems natural to add to the action of GR a term which will, in the linearized approximation, 
give a mass to gravitons without modifying the kinetic terms coming from GR. However, no scalar 
term can be build out of the metric without derivative. Therefore, either one gives up the requirement 
for the invariance under the whole group of diffeomorphisms, or one adds dynamical fields to the the- 
ory. In fact, both approaches may appear equivalent through the Stuckelberg formalism [88]. Indeed, 
any action containing the Einstein-Hilbert term plus a non-gauge-invariant function of the metric can 
be thought of as the unitary gauge description of some gauge-invariant theory containing one scalar 
field for each broken diffeomorphism. Theories of massive gravity discussed in this thesis are of this 
kind; they are based on a gauge-invariant action containing a function of the metric derivatively cou- 
pled to four scalar fields [89], along with the usual Einstein-Hilbert term. Gravitons acquire a mass 
due to a mechanism which may be thought of as the low energy limit of the Brout-Englert-Higgs 
mechanism for gravity [30]: the vacuum expectation value of each scalar field may break one coor- 
dinate reparametrization invariance, and gravitons acquire a mass. One should note here that adding 
dynamical scalar fields is not the only possibility to restore the gauge-invariance of the action. For 
instance, bigravity theories [47, 90, 91] assume two dynamical metrics: their action is made of two 
Einstein-Hilbert term along with a gauge-invariant coupling between the metric. 

Instead of breaking the whole diffeomorphism group of GR, one may require the vacuum to break 
only the invariance under the boosts of the Lorentz group. Indeed, one already knows that if the 
Lorentz-invariance is not broken, either the vacuum is unstable (for non-Fierz-Pauly mass terms) or 
there is a strong coupling problem (for Fierz-Pauly mass terms). It has recently been argued that all the 
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standard problems of massive gravity could be avoided just by renouncing the Lorentz-invariance of 
the background. Consequently, theories of massive gravity with spontaneous breaking of the Lorentz 
symmetry have been considered intensively. Their study started a few years ago (for a review, see [28]) 
with discussions about their stability in Minkowski space-time [89, 92]. Because they are free of the 
usual problems encounter in the study of Lorentz-invariant models, they are candidates for healthy 
theories with massive gravitons. They may be thought of as generalizations of a model known as the 
Ghost Condensate [48] (see chapter 2) with which they share several properties, although gravitons 
are massless in the latter. 

For the sake of completeness, one should mention that there are other directions of search. One 
of them is based on topological massive gravity in 2 + 1 dimensions [93, 94], whose main motivation 
comes from quantum gravity considerations (one expects less severe short-distance behavior in 2 + 1 
than in 3 + 1 dimensions). Theories with fourth-order derivatives have also been studied [95, 84]. In 
these models, there is a massive spin two graviton along with the massless one. Unfortunately, such 
models are plagued by ghost instabilities. Another possibility is provided by models where the four 
dimensional graviton is a metastable resonance with a finite lifetime [34, 35, 96]. Such particular 
graviton state are found in brane theories with extra dimensions, as for example in the DGP model 
[32]. 

Outline of this tliesis 

Up to now, several aspects of theories of massive gravity with spontaneous breaking of the Lorentz 
symmetry have been studied, most of them concerning a minimal class of massive gravity models 
[97] for which there is no modifications of Newton's potential. In this thesis, we discuss some other 
theoretical and phenomenological aspects of this minimal class of models. The first part of this thesis 
is merely introductory. It contains two chapters in addition to the present one: 

• Chapter 2 consists of two parts. The first part overviews the ins and outs of Lorentz-invariant 
theories of massive gravity, with the aim of illustrating the problems which arise when trying 
to give a mass to gravitons. Although this discussion is not essential for the understanding of 
the original part of this thesis, it enables one to comprehend why people have started studying 
models with Lorentz symmetry breaking. In the second part of chapter 2, we introduce models 
of massive gravity with spontaneous breaking of Lorentz symmetry. This second discussion 
gives some details about the stability issues in these models, which lead to a minimal class 
of massive gravity theories possessing only two massive modes with helecities ±2. It is this 
minimal class of models that will be considered in the rest of the thesis. 

• Chapter 3 is devoted to the weak field limit of the massive gravitational field. This discussion 
has two aims. The first one is to introduce a sub-class of minimal models for which there is no 
modifications of Newton's potential. The second one is to provides the equations which will be 
the starting point of chapter 4. 

The second part of this thesis contains all the original contributions to the study of massive gravity 
theories with spontaneous breaking of the Lorentz-invariance. It contains four chapters: 

• In chapter 4, we discuss the first original contribution of this thesis, consisting in the study of a 
physical instantaneous interaction present in massive gravity theories. In GR, the gravitational 
potentials are all instantaneous potentials, as in classical electrodynamics. However, these in- 
stantaneous contributions cancel each other in observables, leaving the theory free of physical 
instantaneous interactions. For the models considered in this thesis, these subtle cancellations 
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are spoiled by the presence of four scalar fields which break Lorentz-invariance. One then 
expects these models to have physical instantaneous interactions. The existence of such inter- 
actions is related to the presence of a mode with the dispersion relation = which could be 
interpreted as a mode with an infinite propagation velocity. 

Beside these theoretical expectations, the presence of a superluminal interaction is required by 
a crucial feature of black holes: it has been shown that black holes do possess hair in massive 
gravity models with Lorentz-symmetry breaking [98]. This violation of the "no-hair"theorem 
of GR is supposed to be the consequence of the presence of physical instantaneous interactions 
who could carry informations outside of the black hole horizon. 

The aim of this chapter is to demonstrate the existence of a physical instantaneous interaction in 
massive gravity theories with spontaneous breaking of Lorentz-invariance. A concrete example 
will be studied, consisting in an instantaneous frequency shift of a light beam by a distant 
gravitational source. 

• In chapter 5, we obtain the exact static vacuum spherically symmetric solutions of massive grav- 
ity equations. In GR, this solution plays a crucial role. First, this solution describes the metric 
outside of spherical non-rotating bodies and gives rise, in the weak field limit, to the Newtonian 
gravity. It provides therefore a useful approximation in many astrophysical situations. Second, 
the Schwarzschild solution describes the result of a gravitational collapse, the black hole. 

Black holes are perhaps the most interesting objects to constrain alternative models of gravity, 
since by reconstructing the metric around an astrophysical black hole one should be able to test 
whether it has the Schwarzschild or Ken- form, and therefore probe GR in the strong field limit. 
In massive gravity theories with breaking of the Lorentz-invariance, the properties of black 
holes are expected to be different. In particular, rotating black holes are certainly modified, and, 
more generally, black holes are expected to have hair [98]. The possible existence of black hole 
hair in massive gravity models suggests that there might exist spherically symmetric solutions 
other than the Schwarzschild one. 

The aim of this chapter is to present a new class of vacuum spherically-symmetric solutions 
in massive gravity theories. These solutions constitute the second original contribution to this 
thesis. One will see that, in addition to the Schwarzschild radius, these solutions depend on 
one more parameter, called the scalar charge S. At zero value of this parameter the standard 
Schwarzschild solution is recovered, while at non-zero values of S the Schwarzschild metric 
gets modified. The modified solutions are non-linear at all distances; they cannot be obtained 
in the linear approximation. These new solutions may have event horizons and are, therefore, 
candidates for modified black holes. Both analytical and numerical examples of such modified 
black holes are discussed in this chapter. 

• Chapter 6 contains the third and last original contribution of this thesis. Given that massive 
gravity models with spontaneous breaking of Lorentz-invariance passes the most obvious con- 
straints, one may wonder if it reproduces correctly more subtle parts of modern cosmology, in 
particular, the theory of structure formation. This is the question which is addressed in this 
chapter. The answer is not obvious a priori since the vacuum in these models contains the 
condensates of four scalar fields whose perturbations mix with the matter density perturbations. 

Cosmological perturbations are studied in the minimal class of massive gravity models. They 
consist of two parts. The first part behaves identically to the perturbations in GR while the sec- 
ond, "anomalous"part is proportional to an unknown function ^'o(a:*) of the space coordinates 
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which arises as an integration constant. The growth of the "anomalous"perturbations depends 
on the value of a single parameter 7. For — 1 < 7 < they grow slower than the standard ones, 
so that the latter dominate, while at 7 = 1 they cancel out. Thus, at least in these two cases 
the perturbations behave in a standard way and massive gravity models are consistent with the 
formation of structures. This is the main result of this chapter. 

• Finally, chapter 7 contains the general conclusion of this thesis. 
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2 

9^assvvt gravity 



One may wonder whether it is possible to construct a well-defined theory describing massive gravi- 
tons. If such theory exists and if the mass m of the gravitons is small enough, or equivalently if its 
Compton wavelength A ~ vrT^ is large enough, one may not be able to distinguish this theory from 
GR in the cunent observations of the gravitational field. All effects arising because of the mass m 
would be small corrections to Einstein's theory. Unfortunately, things are much more complicated. 

In this chapter, one review the problems and issues of theories with massive gravitons. In section 
2.1, one discuses Lorentz-invariant models in their linear regime with the aim of demonstrating that 
they cannot lead to satisfying theory of massive gravity. Although this discussion is not essential 
for the understanding of the original part of this thesis, it enables one to comprehend why people 
have started studying models where the Lorentz symmetry is dynamically broken by the vacuum 
expectation value of four scalar fields, models who are introduce in section 2.2. This second discussion 
gives the details of the models which will be considered in the rest of the thesis. The original part of 
this thesis concern a minimal class of massive gravity models who will be introduced at the end of 
this chapter. 

2.1 Lorentz-invariant massive gravity 

Theories of gravity are in general non-linear as GR, since the gravitational field is a source of gravity 
itself. Such theories are complicated and it is difficult to find analytical solutions to them (we only 
have a few analytical solutions to GR equations). Therefore, when a solution is known, it is common 
to analyze small perturbations around this solution. For example, the Solar System is well described in 
GR by considering small linear perturbations around Minkowski background (this will be discussed in 
details in section 2.1.5). For this reason and because linear approximation reveals the particle content 
of a theory, massive gravity where first considered in the linear approximation, just by adding mass 
terms for the perturbations around Minkowski vacuum. 

First attempts in constructing a model describing massive graviton were related to a Lorentz- 
invariant generalization of Einstein theory known as the Fierz-Pauli (FP) theory of massive gravity 
[29]. But before discussing the FP theory, let us first have a look at the most general theory of Lorentz- 
invariant massive gravity in order to prove that, at the linearized level, the FP theory is the only healthy 
Lorentz-invariant theory describing massive gravitational waves. For this purpose, consider small 
fluctuations about Minkowski space-time parameterized as follow 

At the quadratic level in these small fluctuations, the action for Lorentz-invariant massive gravitational 
waves is simply given by the Einstein-Hilbert action plus the most general Lorentz-invariant mass 
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terms 



where a is a dimensionless constant and h = h^. Indices are raised and lowered with the Minkowski 
metric. In what follows the Lagrangian densities will be defined up to an overall Mpj factor. £eh 
includes the gauge invariant kinetic terms for the metric fluctuations (see appendix A.3.1 for details 
about the linearization procedure) 

Since the gravitational field is assumed to be zero on the border of the integration domain, terms with 
total derivatives are neglected. The mass terms break the gauge invariance of the theory. Hence, it is 
tempting to enlarge the field content of the previous action so as to restore the gauge invariance. This 
can be done through the Stiickelberg formalism [88] (for a recent review on the subject, see [99]). Let 
us introduce four new scalar fields and replace h^y in the Lagrangian by the combination 

With this definition, the Lagrangian for the massive gravitational waves becomes 

2 

2 

~ [d^ed.iu + (1 + 2a) {d^i.f) . (2.1) 
Now, it is straightforward to see that the theory is invariant under the following gauge transformations 

h^j,v hf_iu + d^(iu + du'ifj., Cfj. ^ CiJ. + 'ifi, (2-2) 

which correspond to an infinitesimal coordinate transformations x'^ ^ — <;^. This gauge-invariant 
formulation is completely equivalent to the previous one, since it is always possible to go to the unitary 
gauge = —^^ where there is no scalar fields = 0. Hence, despite the presence of the scalar 
fields, this Stiickelberg trick do not introduce new degrees of freedom. These scalar fields are known 
in the literature as "Goldstone fields". 

Quadratic Lagrangian reveals the particle content of a theory. In GR, there are only two propa- 
gating modes because of the gauge-invariance. These degrees of freedom correspond to the massless 
spin-2 graviton, with the helicity states ±2. Since mass terms break the gauge invariance in the metric 
sector of the massive theory (2.1), extra propagating degrees of freedom are expected to emerge. The 
massless graviton of GR should turn into a massive spin-2 field which has five helicity states: two 
tensor modes with helicities ±2, two vector modes with helicities ±1 and one scalar mode with helic- 
ity 0, according to the irreducible representations of the Poincare symmetry group of massive spin-2 
fields. Therefore, besides the usual tensor modes describing helicities ±2 states, the massive theory 
should posses at least three more degrees of freedom corresponding to the helicities 0, ±1. However, 
the situation is more complicated because of the scalar sector of the theory. 

One of the greatest problem when dealing with a theory of gravity which has more dynamical 
modes than GR is to guarantee that those modes do not imply instabilities, i.e. the vacuum should 
be stable against small perturbations. There are basically two kinds of highly dangerous instabilities. 
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The first one is the classical instability characterized by cj^ = — /c^ + m^. Modes with such dispersion 
relation are stable at low momenta < A;^ < m^. However, they grow exponentially at arbitrarily 
high momenta < k"^ and are therefore unacceptable. The second type of pathological behavior 
is known as ghost. Ghosts have healthy dispersion relations (w^ > 0), but they lead to an instability 
because they have the wrong sign in front of their kinetic terms, and hence negative energy. Therefore, 
it is possible to endlessly create ghost fields along with "healthy"fields without contribution of external 
energy [100, 101]. In order to clarify the particle content of the quadratic Lagrangian, let us introduce 
the transverse and longitudinal decomposition. 

2.1.1 Transverse and longitudinal decomposition 

A convenient way to study the fluctuations about Minkowski space-time in metric theories of grav- 
ity consist in decomposing the metric perturbations into transverse and longitudinal fields [102]. In 
massive gravity, this decomposition is generalized to include the four Goldstone fields 

/loo = 2(/?, 
hoi = Si + diB, 

hij = 2^>5ij - 2didjE - diFj - djF^ + H,j, 

ip, B, ^, E, ^0 and ^ are scalar fields under tree-dimensional rotations ' . Si, Fi and are transverse 
vector fields satisfying 

dA = diF, = diCT = 0. 

Summation over spatial indices is performed with the Euclidean metric. Because of the previous 
relations, each of these three vectors possess only two independent components. Hij is a transverse 
and traceless tensor field 

diHij = Hii = 0, (2.4) 

which possesses only two independent components. Therefore, in total there is one tensor perturbation 
Hij consisting of two components, tree vectors Si, Fi, ^f, consisting of two components each, and 6 
scalars ip, B, E, ^. Since gauge-invariant quantities play important role in gravity, it will be 
useful to introduce gauge-invariant fields. One vector and two scalar perturbations are gauge degrees 
of freedom. As a consequence, there are only two gauge-invariant vector fields 

Wi = Si + Fi, Ui = - Fi, 

and four scalar gauge-invariant fields 

^ = ip- E - B, E = C-E, E'^ = f-E-B, 

and The tensor perturbation Hij is also gauge-invariant. Over-dot denotes the time derivative. 

This decomposition is a powerful tool to explore any metric theory of gravity. Indeed, let ef^ 
and ef ^ be two tree-dimensional, orthogonal, unit-vectors perpendicular to the momentum ki of the 
gravitational field. Then, any transverse and traceless tensor Hij can be decomposed as 

H,, = h^e^ + h^e^, 

'The four scalar fields are scalars under the rotations of the three dimensional space. By rotations here, we mean 
the generalized rotations which include the rotations of the internal space along with the usual rotations of the three 
dimensional space. Such rotations leave Lagrangian (2. 1) invariant. 



--Co, 



(2.3) 
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with the two following polarization tensors 

It is straightforward to see that this decomposition verifies the four relations (2.4), and that these 
polarization tensors satisfy e^^'e^^' = S^"'''\ where summation over (a) indices is performed with the 
Euclidean, two-dimensional metric. Moreover, any transverse vector field vi could also be written as 

diVi = ^ Vi = V(a)ef'' ViVi = Vt^a)V(a)- 

The notations introduced here will be of constant use in this thesis to describe perturbations in 
Lorentz-invariant models of massive gravity as well as in Lorentz- violating theories. They will be 
generalized as to describe perturbations about a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) 
background in chapter 6. 



2.1.2 Getting around ghosts 

The previous notations enable one to simplify the analysis of the field content of the theory (2.1) 
because they preserve the rotational invariance of the Euclidean space. Therefore, if h^j,^ and are 
decomposed according to the notations introduced previously (2.3), the Lagrangian of the massive 
gravitational field decomposes into scalar, vector and tensor parts 

— -^tensor ~l~ -^vector ~l~ -^scalar • 

These tree sectors could then be consider separately. 

Let us start with the tensor sector. This sector describes the usual transverse and traceless graviton 
modes already present in GR (although they are massless in GR). They correspond to degrees of 
freedom with helecities ±2, and have a mass equal to m 

-^tensor = ('-' + ^'^) ^ij ^ 

where □ = d'^d^. This sector has healthy kinetic term and is free of tachyons provided that > 
0. Indeed, by making use of the decomposition of Hij into its two polarization states, the tensor 
Lagrangian reads 

Aensor = "^^x + m'^) hx - (□ + m^) 

Hence, we will assume that > in what follows to avoid tachyon instabilities. 
The vector part of the theory is described by the following Lagrangian 

1 



•^vectOi 2 



As can be seen from this relation, Wi is a non dynamical field since it appears without time derivatives 

7(a) 



in the Lagrangian. This field can therefore be integrated out. Indeed, recalling that wf = wf. and 



that WiCFi = W(^a)(^{a), the variation of the action with respect to zU(^a) gives a vectorial equation- 

= [dl — m^) uji — m^&i, 



^The variation of the quadratic Lagrangian with respect to ■cu^a) gives (9^ — m^) ci7(a) — m^(j(„) — which is a vectorial 
equation with indices (a) — 1, 2. Since the vector fields are transverse, -cui = tU(a)e*°' and ai — o-(„)eJ"', this equation 



implies (9^ — m^) -cui — m^&i = which is a vectorial equation with indices i = 1,2,3. 
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which can be used to express Wi in terms of ai in the Lagrangian. Note that in GR, this equation 
implies that Wi = 0, and there are then no propagating modes of helicity ±1. For the massive case, 
the vector part of the theory becomes 



In the momentum space 



dl 



(Jidlai 



dl 



>0, 



where k = ykf is the three-dimensional momentum. As a consequence, the term with time deriva- 
tives has the correct sign in the previous Lagrangian provided that > 0. Hence this conditions 
does not only guarantee that the tensor degrees of freedom are not tachyons, but also that there are 
no ghosts (fields with wrong sign of the kinetic terms) in the vector part of the theory. The following 
redefinition field 



m . 



291 



dl 



enables one to write the action of the vector sector in the canonical form 

1 



(2.5) 



The vector sector becomes dynamical since the Lorentz - invariant mass terms introduced in the 
quadratic Lagrangian (2. 1) contain kinetic term for the transverse vector fields. Hence, the vector 
Lagrangian describes two propagating modes of mass m and helicities ±1. 

Finally, the scalar sector describes the helicity degree of freedom. For a general value of a, there 
are two dynamical modes in the scalar sector whose Lagrangian is given by 



scalar 



m 



-rn^a ( <I> 



3^ - dll 



In GR, m = and <I> is a non-dynamical field. Still, this field appears hnearly in the Einstein-Hilbert 
Lagrangian with consequence that it can not be just integrated out as it was the case for the vector 
field Wi. Indeed, varying the action of GR with respect to $ gives an equation which is actually a 
constraint on ^' 



0, 



and which could therefore not be used to integrate <^ out of the Lagrangian. This situation occurs in 
GR because ^ is in fact a Lagrange multiplier imposing a constraint on ^. Using this constraint, the 
equation obtained by varying the Lagrangian with respect to ^' becomes = 0. Hence, there is no 
propagating modes in the scalar sector of GR. 
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The situation is different for the massive gravitational field. The Lagrangian of the scalar sector 
can be written as follow 



/:scaiar = 6^^ + 2^5/ {2A + 213 + 2E - ^ ] - m 



2 



a{A-?,-^ -dlE)\ (2.6) 



where A = ^ — E^ and B = E^ — E. In the massive case, the role of $ is played by the field A, 
which is not a dynamical fields since it enters the action without time derivatives. Let us first discuss 
the class of models for which a ^ —1. As can be seen from (2.6), A is not a Lagrange multiplier; it 
appears quadratically in the Lagrangian. Therefore, provided that a 7^ —1, the equation obtained by 
varying the Lagrangian with respect to A gives a relation which can be used to integrate A out. Since 
there is no terms with B and two time derivatives, this field is not dynamical either implying that it 
can also be integrated out. The equations obtained by varying the action with respect to A and B are 

= 2dl^ -rr?{a + l)A + rr?a (3^' + dlE) , = 4^' + m'^B. 

By making use of these two equations, the Lagrangian of the scalar sector becomes 

^scalar = 6-i'<f + 2^df I ^ - ^ + 2H - + (3^1/ + a^S) | 

-m" {^6ij + didjEf - m^^^ (3^ + S^H)' . 

This Lagrangian looks quite complicated. But at this point it is clear that both ^' and E are dynamical 
fields. One of them is a ghost, i.e. a field with the wrong sign of the kinetic terms. This can be seen 
from the terms with time derivatives. These terms can schematically be written as 



62 



where a and b are some constants. Depending of the sign of a, either or H is a ghost, since these 

a 

two fields have an opposite sign in front of their kinetic term. Therefore, Lorentz-invariant massive 
gravity models with a 7^ —1 have pathological behavior at arbitrarily high spatial momentum. 

The Lorentz-invariant massive gravity model characterized by a = — 1 is known as the Fierz-Pauli 
(FP) theory [29]. In this particular model, there is no quadratic terms in A in the Lagrangian (2.6). 
This field is therefore a Lagrange multiplier as in GR, imposing the following relation between ^ and 

= 29^* - (3^ + dlE) . 

This constraint kills one dynamical field out of two, leaving only one helicity-0 mode in the La- 
grangian. Using this constraint, and after integrating out B as previously, the Lagrangian of the scalar 
sector becomes 

Acalar = "6^ (□ + m^) I-. (2.7) 

This Lagrangian is healthy and describes one scalar propagating mode corresponding to the helicity 
mode of the massive graviton. Therefore, the Lorentz-invariant theory (2.1) with > and a = —1 
describes five dynamical modes corresponding to the five helicity states of a massive spin-2 graviton. 
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2.1.3 The Goldstone sector 



The previous conclusion concerning the vector and scalar sector could have been found directly from 
the study of the Goldstone sector of the theory. The four Goldstone fields have been introduced to 
restore the gauge invariance of the massive gravitational field. Yet, there is another motivation for 
the introduction of these fields: these scalar fields give a convenient way to single out and study 
the dangerous degrees of freedom of the theory [87, 89]. Since the Goldstone fields only appear in 
the vector and scalar sectors, let us concentrate on these two sectors, described schematically by the 
following Lagrangian 

/:vector + /:scalar = Ku (^'^""^ + Mf h^p + 2Mi^^h^,i^ + ^^.D^'^iu, 

where V is the usual two-derivative graviton kinetic operator coming from the Einstein-Hilbert La- 
grangian, Mh accounts for the metric mass terms. Mm is a one-derivative operator describing the 
mixing between the metric and Goldstone fields, and finally Dq is a two-derivative operator describ- 
ing the kinetic terms for the Goldstone fields. The field equations read 




0. 



Here, one assume that the gauge symmetry is fixed in the gravity sector, by choosing for example 
the de Bonder gauge^. Then, the Einstein-Hilbert operator V is non-degenerate and the first set of 
equations may be used to express the metric perturbations in terms of the Goldstone fields 

= ^ — I — ^r^Ca- 

This therefore gives a closed set of equations for ^q. 



For momenta larger than the mass, k ^ m, the mass matrix could be neglected with respect to V. 
Then, the second term in the last equation is of order while the first term is of order (u^, k'^) x m^. 
Consequently, when working atk^m, the field equations reduce to 

D^^p = 0, 

equations which are given by the Goldstone sector of the theory 

This short calculation shows that neglecting the metric perturbations is legitimate as long as modes 
with k"^ » m? are considered (except for the tensor modes which are completely described by the 
metric perturbations). This result is compatible with the previous analysis of the vector and scalar 
sector of the Lorentz-invariant theory. Indeed, with the decomposition (2.3) into transverse and longi- 
tudinal fields, the Goldstone sector of the Lagrangian (2.1) reduces to 



2 

m 

-^Gold. = 



iJu^J + 2{l + a) Ul+ {dUY) - mo? -Ui] + 2 (1 + 2a) iod^ 



2 



^^In the de Dander (or harmonic) gauge, the perturbations of the metric satisfy d'^ — Tii_t,yh) = 0. 
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This relation makes it clear that the .^^^-terms are perfectly healthy provided that rn^ > 0, and describe 
two propagating modes of helicity ±1. This vector sector is exactly the high-momentum k"^ » 
limit of the vector sector (2.5) which has been discussed previously. 

In order to discuss the scalar part of the Goldstone sector, let us focus first on the theories for 
which a / — 1. In the four-momentum space, the scalar-Goldstone sector is given by 

r { hP-\Mi^^\ M / -2(l + a)w2 + A;2 i(l + 2a)a;A; \ 

/:Goid.sc. = — ( Co j, M=y _,(i + 2a)u;fe u;^ - 2 (1 + a) /fc^ j " 

The dispersion relations of these two scalar fields are obtained by requiring that det M = 

detM = -2(l + a) (a;2_A:2)2 = ^ u? = k^. 

This relation clearly shows that there are two dynamical modes in the scalar-Goldstone sector, charac- 
terized by = k^ . Noting that det M = A+A_, where \± are the eigenvalues of M, the dispersion 
relations of the scalar-Goldstone fields could also be determined by finding the solutions to the fol- 
lowing equations 

, , TrM ± V(TrM)2 -4detM 

A± [u?) = with A± = ^' . (2.8) 

Two conditions have to be met in order for the Goldstone scalar sector to be free of instability. The 
first one requires tj^ > to guarantee the absence of classical instability. This condition is satisfied 
here. The second condition consist in requiring that the eigenvalues of M be positive definite on their 
roots ijj± 



> 0. (2.9) 



If satisfied, this second condition implies the absence of ghosts. Therefore, to avoid ghosts it is 
necessary that each eigenvalue X± has one zero, otherwise one of the eigenvalue will have two zeros 
and one of them will necessary violate the condition (2.9). This last observation enable to understand 
whether there is a ghost in the Goldstone scalar sector or not. Indeed, let us relate the dispersion 
relations to the eigenvalues of M. 

• If a > —1/2, the two dynamical modes are both solutions of the equations A+ = 0, since the 
trace of M is negative for these modes. Therefore, A_ has no zero while A+ has two, implying 
that one of them is necessarily a ghost. 

• If a < — 1 /2, the trace of M is positive when lu'^ = k'^ implying that the two dynamical modes 
are solutions of A_ = 0. Hence, there is also a ghost for these value of a. 

• If a = —1/2, A+ and A_ both have one zero on lu'^ = k'^. 

This discussion makes clear that, depending of the value of a (and recalling that we have excluded a = 
—1 from the beginning), there are two groups of solutions. First, if a = —1/2, the two propagating 
modes are associated to different eigenvalues, and it is straightforward to see from the matrix M that 
in this case is a ghost-like field while ^ has a healthy kinetic behavior. For the second group of 
solution characterized by a ^ —1/2, the two propagating modes are associated to the same eigenvalue 
implying that one of the modes is a ghost^. One should therefore conclude that there is always a ghost 



"^This situation is similar to what happens in theories with four derivatives acting on a single scalar field. It is well known 
that such models posses ghost [103, 104]. 
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like mode in Lorentz-invariant massive gravity. Yet, there is one special case for whicii this conclusion 
does not hold. Indeed, the only way to get rid of this pathological behavior is to choose a = — 1. In this 
particular model corresponding to the FP theory, there is no dynamical scalar mode in the Goldstone 
Lagrangian. 

The previous discussion shows that the study of the Goldstone sector is a powerful tool to analyze 
the field content of a theory. Yet, we have just shown that for a = —1 there are no helicity mode in 
the Goldstone sector while we have explicitly found such a mode in the previous section (2.7). This 
apparent contradiction is a particularity of the FP theory, in which the kinetic terms for the scalar 
mode are originated in the mixing between the metric perturbations and the Goldstone fields (recall 
that ^ is a Lagrange multiplier which impose a constraint between these terms). Therefore, the only 
scalar terms present in the Goldstone sector of the FP theory are mass terms which are negligible with 
respect to the kinetic terms coming from the mixing, and there is no contradiction. This particularity 
of the FP theory is directly related to the vDVZ discontinuity of this theory, which will be discussed 
in the next section. 

Before concluding the discussion of the Goldstone sector, it is important to remember that this 
sector is the high energy limit of the massive theory, corresponding to the k'^ ^> m? limit. There- 
fore, by studying the Goldstone sector we could miss some other modes with more slowly oscillating 
behavior. This remark should be kept in mind for the coming study of Lorentz- violating models. 



2.1.4 Interaction between two sources 

To agree with Solar System observations, the FP theory of gravity should reproduce GR predictions 
in the limit of vanishing graviton mass, since GR is quite successful in describing the whole Solar 
System. The predictions of the massive theory can be compared to those of GR by studying the 
interaction between two massive bodies, or between one massive body and light. For this purpose, 
let T and T' be two small, conserved energy-momentum tensors, corresponding for example to two 
massive bodies. These two energy-momentum tensors both source gravity through the linearized 
gravitational equation. For both GR and the FP theory in the unitary gauge (in which there is no scalar 
fields), the gravitational field produced by T is given by 

where "P^J^^ is the graviton propagator, obtained by inverting the equations of motion (EoM)''. The 
interaction between this gravitational field and T' is given by the action of this gravitational source 
[107] which contains all interaction terms between the metric and T' 

From this relation, it is straightforward to understand that it is sufficient to compare the graviton 
propagator in both theories to see if GR and the FP theory have similar predictions. After a four- 
dimensional Fourier transformation, the graviton propagator reads 

kxk^ ^^^"^"f^ ^M/3^i^« V^lur]af3 + ^-dependent terms] , 



FP : V-}, ^ 



^^"^ M^i {kxk^ - m2) 



2 

1^^la■r]u|3 + Tlfj.p'ni'a " -^VtivVaiB + fc-dependcnt terms 



""This could be done with the spin formalism of Barnes and Rivers [105, 106], but also with the transverse and longitu- 
dinal decomposition introduced previously. Indeed, the propagator could easily be calculated by adding sources in each of 
the three sectors of the theory, and then inverting the EoM. 
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The A;-dependent terms are gauge-dependent terms. They give a null contribution to 5int. since they are 
contracted with conserved energy-momentum tensors. Apart from these terms and from the presence 
of a non-zero mass, we see that the propagator in the FP theory differs from Einstein's one by a factor 
of 2/3 in front of the third term. This difference, which is mass-independent and therefore remains 
when the limit m ^ is considered, is known as the van Dam-Veltman-Zakharov discontinuity 
[108, 109, 1 10], or vDVZ discontinuity. This discontinuity implies that the addition of FP mass terms 
in the action of the gravitational field changes drastically the interaction between two arbitrary sources, 
even for arbitrarily small masses. 

The origin of this discontinuity can be traced back to the scalar sector of the FP model. Indeed, 
as discussed previously, the FP theory possesses five propagating modes corresponding to the five 
helicity states of the massive graviton: two tensor modes already present in GR, two vector modes 
and one scalar. It follows from the definition of the canonical vector field (j| that the vector Lagrangian 
(2.5) reduces to zero in the limit of vanishing graviton mass; i.e. the vector modes get kinetic terms 
from the Goldstone sector and decouple in the limit m ^ 0. Therefore only the standard contribution 
coming from GR survives in the limit of vanishing mass. The situation is completely different in 
the scalar sector. The kinetic terms for the scalar mode come from the mixing between h^y and the 
Goldstones. For this reason, the scalar mode does not decouple in the limit of vanishing mass, and 
the Lagrangian of the scalar sector does not reduce to zero when m — > 0. Therefore, the FP theory 
reduces to a scalar-tensor theory of the Jordan-Brans-Dicke type [111]. Indeed, in the presence of a 
conserved source T, the Lagrangian (2.7) of the scalar sector generalizes to 



where T = 7^. We conclude from this Lagrangian that, in the limit m — > 0, the ^ field becomes 



where and (□ + m?) ^ have to be understood as the Green functions of the d'Alembert and 
Klein-Gordon equations, respectively. The first term in this relation is the usual GR contribution to 
^f, while the second, proportional to the trace of the energy-momentum tensor, is a new contribution. 
Along with the contribution coming from $ and which is not independent, this new term is responsible 
for a l/3rjfj_ur]ai3 contribution in the propagator, which combined to the contributions of the tensor 
modes gives the abnormal 2/3 factor. 

The third term in the propagator is responsible for the coupling between the traces of the two 
energy-momentum tensors. Therefore, this extra term will contribute in the interaction between two 
massive bodies, but not in the interaction of a massive body and a light beam since 7^' is identically 
zero for light''. To illustrate the prediction differences between the Einstein theory and the FP one, let 
us assume that T and T' correspond to two static massive bodies. For such bodies, the only non-zero 
component of their energy-momentum tensor is the zero-zero component. Therefore, the gravitational 
interaction between two static massive bodies is given by 

T T' 



t^uap 3 ^2^ (^^^A _ ^2) 



*The energy-momentum tensor of light is given by T = diag. (p, — p, — p, — p). This tensor is trace-less. Indeed, p and 
p satisfy the following equation of state 3p = —p. 
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This implies that, in the most favorable situation corresponding to the m ^ limit, the FP theory 
predict an interaction which is 4/3 stronger than GR, unless the constant Mpj in the FP theory is 4/3 
bigger than its value in GR. This will re-conciliate the FP model with the Newtonian potential, and 
therefore with the observations of gravitational interaction between massive bodies, but will give a 
wrong angle deflection for light beam passing close to the Sun. Indeed, the energy-momentum tensor 
of light is trace-less. Therefore, the interaction between a massive static body and a beam of light is 
given by 



2Too7^o 



Then, if Mpj is 4/3 bigger in the FP theory than its value in GR, these last relations implies that the 
interaction between a massive body and a light beam is 3/4 the value predicted by GR and confirmed 
by observations. 



2.1.5 Breakdown of linearity 

The vDVZ discontinuity seems to indicate that the FP theory can be ruled out by Solar System obser- 
vations, since this theory predicts a deflection of light by the Sun which is 25% smaller than the value 
predicted by GR and confirmed by observations. In fact, this conclusion is not trustworthy [112]. 
Indeed, the appearance of the vDVZ discontinuity reUes on a linear approximation of a non linear 
theory. In GR, this approximation is valid in the entire Solar-System, while it is possible to show that 
the linearized FP theory is not valid in the Solar System for arbitrarily small masses. 

In both theories, GR and FP, the metric around the Sun can be linearized provided that there is a 
small parameter e in the metric to support the Taylor expansion 

ds^ = [g^,u + h^u + (e^)] dx^dx'^, | VI ~ ^■ 

Then, if e ^ 1 the second order terms could be neglected as compared to the first order terms, and 
the gravitational interaction is correctly described by the linearized approximation. In GR, the metric 
around the Sun is given to a very good approximation by the Schwarzschild solution, and the linear 
approximation is valid at distances 

r > 2GMq ~ 3km. 

Therefore, in GR the gravitational field around the Sun is correctly described by the linearized ap- 
proximation for distances much larger than 3 km, which corresponds to the Schwarzschild radius of 
the Sun and which is actually much smaller than the radius of the Sun. 

In the FP theory, the small parameter controlling the validity of the linearized approximation is 
ry/r, where is the Vainshtein radius [112]. As in GR, the linearized approximation is valid for 
distances much bigger than r„ 

For a graviton mass as big as m ~ (lO^^cm) ^ (which corresponds more or less to the size of the 
Solar System), the Vainshtein radius ~ 10*^ km is of the same order as the Earth orbit. For a much 
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smaller mass corresponding to the Hubble scale, m ~ Hq , the Vainshtein radius is ~ 100 kpc 
which is bigger than the actual size of the Milky way. 

The Vainshtein radius is singular in the graviton mass, and therefore the higher order corrections 
to the linearized theory are also singular in the graviton mass. This has for consequence that the linear 
theory breaks at a distance much larger than the Schwarzschild radius of GR. One has to conclude 
from this that the linearized FP theory can not describe gravity within the Solar System if the mass 
of the graviton is too small, since then the distance below which non-linearities become important is 
too big, and the light passing not far from the Sun will feel the non-linearity of the FP theory. This 
could solve the problem arising from the vDVZ discontinuity. Indeed, with the account of the full 
non-linear theory, one may hope that the non-linear completion of the model will modify the theory 
in such a way as to smooth the massless limit, and reconcile the model with observations. One can 
for example imagine that the full non-linear solution is continuous in the graviton mass, and still be 
the summation of singular terms. This is exactly what happens in the DGP model [32], in which 
the higher order corrections affect the scalar sector and decouple it from other modes in the small 
mass limit [113, 114]. It is currently debated whether non- linear completion of the FP theory could 
solve the problems raised by the vDVZ discontinuity, since the full non-linear FP theory is not defined 
unambiguously. Recently, it has been argued by studying the Goldstone sector of FP theory that for 
some non-linear models there is a Vainshtein-like recovery of the solutions of GR [115]. Still, ghost- 
like mode seems unavoidable when accounting for higher order terms [116]. This ghost can be seen 
as responsible for the cancellation of the attraction exerted by the non-decoupling scalar mode [117]. 
But it also indicates that the theory is unstable inside the Vainshtein radius. 

Non-linearities could also solve the vDVZ problem by selecting a non-asymptotically flat vacuum 
[47, 118]. Indeed, there is no mass discontinuity in de Sitter or Anti - de Sitter space-times [1 19, 120, 
121, 122]. 

2.1.6 Strong coupling 

Both GR and the FP theory should be treated as low energy effective theories valid at energies below 
a certain energy scale A, which corresponds to the energy at which an ultra-violet completion of the 
theory is needed. Below A, all quantum effect are negligible. Above this scale, terms with higher 
order derivatives become important and quantum effect could not been neglected. In Einstein theory, 
it is natural for dimensional reasons to identify this UV energy scale with the Planck mass Mpi [37]. 
One would therefore trust GR below this energy scale and for distances bigger than the Planck length 
Mpj^. Above this energy scale, it is supposed that GR has to be extended in a quantum theory of 
gravity. 

In the FP theory, the energy scale at which quantum corrections become important is much lower. 
This is related to the non-linearity problems of the FP theory, which at the quantum level implies that 
the theory becomes strongly coupled at an energy scale [87] 

A=KMp,)'/'. 

Therefore, the FP theory is trustable at distances larger than A^^. For a mass m ~ (lO^^cm) ^ of 
the order of the Solar System size, the effective theory is valid for distances larger than 4.5 km. For a 
much smaller mass corresponding to the Hubble scale, quantum correction are important for distances 
smaller than 10^^ cm, which encompass the entire Solar System. Since Newton's law has been tested 
to sub-millimeter scales [3, 4, 5, 18], this strong coupling effect is one of the biggest problems of the 
FP theory. 
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The energy scale A is related to higher order operators (terms with more than two derivatives) 
which should be present in the UV completion of the theory. Yet, there is another source for such 
terms. Indeed, the Goldstone fields have been introduced to restore the gauge invariance of the theory. 
But actually the transformation law of the metric perturbations is not strictly linear^ 

h^y hf^u + dfj,qy + dy(i^ + d^^adv'^'^ + ... . 

Therefore, the definition of h^y should contains terms quadratic in the Goldstone fields, otherwise the 
theory would not be gauge invariant at the quadratic level. The previous definition of h^y should then 
be completed 

h^,y = h^u - df,iy - dy^f, - d^iadvi"^ + ... . (2.1 1) 

Consequently, the Lagrangian of the FP theory actually contains higher order terms which we have 
neglected up to now. These terms are responsible for the particularly high UV cutoff of the FP theory. 
Let us introduce one vector and one scalar fields through the following relations 

Cm = ii^ + dpi(t)- 

These notations replace four variable by five. Therefore, this introduces a fake U{1) symmetry 

4 ^ + ^M-^' (/> ^ </> - A. 

Then, the kinetic terms of the scalar sector can be diagonalized by the following conformal transfor- 
mation 

which eliminates the kinetic mixing between h^y and 4>. If h^j^y and are decomposed into transverse 
and longitudinal fields according to (2.3), the scalar Lagrangian reads^ 

The U{1) symmetry enable to cancel either ^' or 0, since these fields transform as follows: 

2 

TYl 

— ^A, (j)^4>-\. 

Let A be such that ^' = 0. The higher-order interaction terms for the Goldstone fields are all 
invariant. Therefore, with the addition of these higher-order operators, and recalling that there is an 
overall Mpj factor in the definition of the Lagrangian, the scalar Lagrangian is 

M^i£,,,i,, = -^-<t,^ (□ + m") 4>^ - ^Ucj,- {d^dyck'd^'d^<i^' - {Ucl^-f^ + . . . , (2.12) 

where the following rescaling (j)^ = m^Mpi(/> has been made in order to get the canonically normalized 
scalar Lagrangian. It is obvious from this relation that the cubic interaction terms for (f), coming from 
the definition of h^^, are suppressed by a parameter which is precisely = m^Mpi. Since other 

'This assumes that the gauge invariance of the full non-linear theory is the same as in GR. 
is a gauge invariant field. 
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higher order operators, coming from the non-linear and from the UV completions of the linearized 
theory, are less coupled than the cubic ones [87], A correspond to the strong coupling scale. This 
value of A has been confirmed by direct calculations of the four-graviton scattering amplitude [123]. 

One can imagine that the UV completion of the theory will bring terms which will exactly cancel 
these cubic operators. If such a UV completion exist, the strong coupling scale will be lowered to 

1 /3 

A = (m^Mpi) which is the lowest that can be achieved [87]. For a graviton mass corresponding to 

the Hubble scale, this new cutoff scale is still too low, A = (lO^ km) ^. In principle, this low cutoff 
scale could be raised in a non-flat background due to curvature terms. This occurs in the DGP model 
[124], but apparently such effects are absent in the FP theory because of the presence of higher order, 
ghost-Uke propagating mode [116, 1 17] 

It is worth noting that Lagrangian (2.12) enables one to understand the origin of the Vainshtein 
radius and its link to the strong coupling problem. Indeed, let us consider a static massive body for 
which T^jy = 6^o6uoM6^ (x). Then, if we neglect the higher order operators, it is not difficult to 
convince ourselves from relation (2.10) that (/> is given by 

where the homogeneous solution has been neglected since it does not play any role in the present 
discussion. For r <^ m^^, this solution can be approximated by 4''^ ~ Mr~^. Hence, the cubic 
operators could be neglected as compared to the second order terms if and only if 
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Apart from the coefficient of order one, this is exactly the definition of the Vainshtein radius. 



2.1.7 Instabilities in curved backgrounds 

In curved background, the relation a = —1 is violated by curvature terms. This relation, which 
defines the FP theory in Minkowski space-time, guarantees that the theory is free of pathological 
behavior. Its violation implies that ghost instabilities reappear in the spectrum. Such instabilities in 
curved background are known as Boulware-Deser instabilities [107]. They exist for arbitrarily high 
spatial momenta and are therefore unacceptable. 

The origin of the Boulware-Deser instability is easy to understand from the study of the field cf). 
Let us consider a nearly-flat background'''. If the curvature is small enough, the background could 
be described to a good approximation by the first order perturbation theory about Minkowski solu- 
tion. Then, the perturbations about this nearly flat background will be given by expanding the 
Lagrangian for at the quadratic order in the field X = 4''^ ~ 4b, where (pb is the background. The 
Lagrangian for x will contain four time derivatives. Indeed, cubic terms {d'^(j)Y Lagrangian (2.12) 
give d'^4>b{d'^x)'^ ■ Independently of the way that the FP theory is generalized in curved space-time, 
the quadratic Lagrangian for x reads schematically [116] 

~ i^x) -m X - -j^ [d X) ■ 

'in Ref. [28], these Boulware-Deser instabilities are explicitely found and discussed for a Friedmann-Lemaitre- 
Robertson- Walker (FLRW) background. In Ref. [125], this instability is discussed in a fully non-linear generalization 
of the FP theory through the analyze of the Goldstone sector. 
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This Lagrangian has four derivatives acting on a single scalar fields. It is well known that such model 
possesses a ghost [103, 104] This Boulware-Deser ghost as a coordinate dependent mass given by 
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'ghost 



Since the FP theory is only vahd for energies below A, this ghost instability has no physical meaning 
as long as rrighost > A. Moreover, when approaching Minkowski background, the mass of the ghost 
increases and therefore the distance below which this instability shows up decreases. But this does 
not save the FP theory. Indeed, at distances r S> from a static source, the background solution (f)^ 
is given by (2.13). Then, the ghost mass is of order 



where rghost is the distance at which mghost drops below A. This interval is not empty with consequence 
that FP theory has ghost instability in curved background. 

To summarize, the FP theory is the only Lorentz-invariant model describing massive gravitational 
waves which is stable in Minkowski background. However, its suffers from several problems. Some 
of them might be cured by higher-order terms, or by an UV completion of the theory. Recently, it 
has even been argued that the helicity ±1 modes could be superluminal in slightly curved background 
[126]. Although the FP theory might be less pathological if the vacuum is different from Minkowski 
space-time [119, 120, 121, 122], it seems that the FP theory has too many pathologies to be considered 
as a serious candidate for a theory of gravity. This justifies the study of Lorentz- violating massive 
gravity, since breaking this symmetry is a way to go around all the usual problems of massive gravity 
discussed in this section. 

2.2 Lorentz-violating massive gravity 

Recently, it has been argue that the standard problems of massive gravity could be avoided just by 
renouncing the Lorentz-invariance of the background [89, 92]. Those works were part of a more 
general tendency based on the idea that, within the current observational limits, the Lorentz symmetry 
could be spontaneously broken in the gravity sector' ' . This would appear audacious, as the Lorentz- 
invariance is perhaps the most important symmetry in physics responsible for the giving up of absolute 
time and the emergence of the space-time continuum. Likewise, Lorentz-invariance could be the limit 
of a more general unknown symmetry, just as the Galilean group is the small velocity limit of the 
Lorentz group. 

Up to now, Lorentz-invariance has been tested in many sectors of the Standard Model of Elec- 
troweak Interactions [129, 130, 131, 132]. Yet, the idea that Lorentz-invariance could be broken by 

'"in general, the presence of ghost in theories with fourth-order operators acting on a single scalar fields can be demon- 
strated by reformulating the theory as a standard two-derivative Lagrangian for two independent scalar fields through the 
definition 9 = X- 

"There are some hints that quantum gravity and string theory may both suggest that Lorentz-invariance is broken at high 
energies [127, 128]. 
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and the ghost instability occurs for distances 
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the vacuum expectation value of some field is not conceptually different from the way that any time- 
dependent cosmological fluid actually breaks this symmetry. However, those cosmological fluids are 
not the vacuum state of the theory since they carry energy and dilute away with the Universe's ex- 
pansion. An example of such a cosmological fluid is given by the Cosmic Microwave Background 
(CMB). Indeed, any observer at rest with respect to the CMB would conclude that the CMB is highly 
homogenous and isotropic, while other observers working in other reference frames would observe a 
dipole because of theirs own motion with respect to the CMB's photons. 

Different models have emerged with spontaneous breaking of the Lorentz-invariance [127, 133, 
134, 135, 136, 137, 138, 139, 140, 141]. One common feature of all these models is that they "dy- 
namically"break the Lorentz-invariance. Indeed, if the preferred frame were to be fixed by an external 
mechanism, it will violate the general diffeomorphisms invariance of the theory, and the divergence of 
the matter energy-momentum tensor will not vanish'-. Therefore, the Lorentz symmetry has to be bro- 
ken dynamically by fields, whose EoM will imply the conservation of the matter energy-momentum 
tensor. 

The Ghost Condensate model [48] is an example of a theory in which the Lorentz-invariance 
is broken by the time-dependent vacuum expectation value of a scalar field. The condensation of 
this field leads to infrared modifications of the gravitational interaction, without involving any of the 
pathologies discussed previously. Although gravitons are massless in this model, the Ghost Con- 
densate shares several properties with the Lorentz-violating massive gravity theory discussed in this 
thesis. Another interesting class of model is given by bigravity theories [90, 91]. There are two dy- 
namical metrics in those models, one being not Lorentz-invariant. Therefore, bigravity theories posses 
two graviton fields associated to each metrics, one being massive while the other massless. 

In those two last examples, as in Lorentz-violating massive gravity, the Lorentz symmetry is 
broken by a mechanism which could be seen as an analog of the Brout-Englert-Higgs Mechanism for 
gravity. Lorentz-violating models of massive gravity will be introduced in the rest of this chapter. The 
purpose of this discussion is the motivate the particular class of models which is the subject of this 
thesis and which can be thought of as "minimal models"of massive gravity. 

2.2.1 Brout-Englert-Higgs Mechanism for gravity 

Gauge interactions are mediated by massless particles. Photons are one example of such gauge par- 
ticles corresponding to the unbroken U{1) gauge symmetry of electromagnetic interactions. Upon 
spontaneous breaking of a gauge symmetry, the gauge particle acquire a mass via the Brout-Englert- 
Higgs Mechanism. An example of this mechanism can be found in the Standard Model of Electroweak 
Interactions, in which some of the bosons associated to the electroweak SU{2) x U{1) gauge interac- 
tion become massive through their coupling to the Brout-Englert-Higgs field. 

The gauge group of GR is the general reparametrization invariance associated to a massless spin 
2 particle, the graviton. If the diffeomorphism invariance is spontaneously broken, the graviton is ex- 
pected to acquire a mass. In four dimensional space-time, there are four reparametrization invariances 
associated to the four coordinates. Each of these symmetries could then be individually broken, for 
example by a scalar field with vacuum expectation value depending on one particular coordinate. Such 
scenario could be thought at as the low energy limit of some kind of Brout-Englert-Higgs Mechanism 
for gravity ' \ As an example, let us consider the following Lagrangian which depends on four scalar 

'^In GR, the gauge invariance implies tiiat the divergence of the Einstein tensor is identically zero with consequence that 
the matter energy-momentum tensor is also divergence-less. 

'^The analog of the Brout-Englert-Higgs Mechanism for gravity will include the UV completion of the theory. 
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fields (fP , (j)^, (jP' and cf)^ through a function V{Q) 

where qafs is a constant metric for the scalar fields''^. The parameter A is a UV cutoff. The gravita- 
tional field equations for this theory read 

1 / 1 

The Minkowski vacuum is then solution to the EoM if the scalar fields are given by 

(fp = at, (f)^ = bx^, (jP' = crc^, (fP = dx^, 

where a, b, c and d are four constants which values are determined by the requirement that the energy- 
momentum tensor associated with these four scalar fields vanishes in the Minkowski background 

2 ^ V ,2 ^ 'P 2 ^ 'P Jl ^ P 

These four equations are satisfied for example if V{Q) = 'P'(Q) = 0. For arbitrary value of the 
metric qap, each of these four Goldstone scalar fields breaks one of the four coordinate invariances of 
the Lagrangian by selecting one coordinate parametrization among all possibilities. Therefore, these 
scalars are selecting one reference frame. 

Rather than breaking the whole diffeomorphism invariance of the theory, this mechanism could be 
used to break some specific symmetries while keeping the residual reparametrization invariance. For 
example, the background is invariant under the group 50(3) of three-dimensional rotations if goi = 
and Qij oc 5ij, because then b = c = d and the Lagrangian is invariant under a 50(3) rotation of the 
fields ^^j'^- Likewise, the invariance under boosts will not be broken if q^y oc rj^y. In these 

two last examples, the invariance under three-dimensional rotations is not broken and the transverse 
and traceless graviton modes with helicities ±2 satisfy the following equation 

0= [n-m^r)Hij, 

which shows that the gravitons are massive provided that V{Q) < in Minkowski background. Yet, 
this example does not represent the most general case. For arbitrary function V, the four constants 
a, b, c and d have to be such that V{Q) = in order for Minkowski space-time to be solution of the 
EoM. In general, the mass of the gravitational waves has to include the second derivative of V in order 
to be non-zero. 

Before introducing the action for massive gravity, let us remark that bigravity theories with spon- 
taneous breaking of the Lorentz-invariance are based on such symmetry breaking as well. Indeed, 
these models includes two Einstein-Hilbert terms along with a coupling function V{g'^"fau) [90, 91]. 
Then, Lorentz-invariance is broken by assuming light propagates at different speed in the two space- 
times. Those models are somehow related to the mechanism presented here, since d^cjf^dycji^qajs 
could be seen as a second space-time metric Then, the speed of light in this second space-time 
will be given by a? qoo / {b"^ qii) while being equal to one in the first metric. In the Ghost Condensate 

'■*For related works, see [30, 142]. 

'^Of course, bigravity theories are completely different from the present example because of their field content. 
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model [48], the metric for the scalar fields satisfy 5^/3 = SaoSpo. Therefore, the invariance under 
boosts is broken in this model. V and V' have to vanish in Minkowski space-time in order to satisfy 
the EoM for this background implying that the gravitons remain massless. Yet, this model has a very 
interesting phenomenology. 

2.2.2 The action for massive gravity 

This mechanism for the spontaneous breaking of the reparametrization invariance of the gravitational 
field can be used to construct a full non-linear theory of massive gravity, with a dynamical breaking 
of Lorentz-invariance by the space-time dependence of four Goldstone scalar fields. Indeed, the idea 
is to add terms to the action of GR in order to describe massive graviton without Lorentz symmetry. 
This could be done by adding an arbitrary non-Lorentz-invariant function of the metric to the usual 
Einstein-Hilbert action 

S = Sgr + with S^ = K^ j A^x ^gT (5°° , , d'') ■ 

Since this function involves only the metric without derivatives, one expects that such terms will 
give a mass to the gravitons. But in this action, the Lorentz symmetry is not dynamically broken 
with consequences that the action is not gauge invariant. Yet, it can be seen as the unitary gauge 
description of a gauge invariant action through the Stiickelberg trick. Indeed, the previous model can 
be generalized in a gauge invariant formulation by adding four scalar fields , 0^, (l? and (\? to the 
action, which then reads [89] 

64, = A* J d^x ^gT (X, V\ W'^) , (2.14) 

with 

X = g^"'^f,4>°^^(|)°, V' = g^"'^^(|P^u(t)\ W'^ = g'^'d^cp'd^cp^ — . (2-15) 

X 

Note that the fields cj)^ have dimensions of length so that X, and W^^ are dimensionless. This 
new action defines a class of models which are gauge and Lorentz-invariant. As previously, A is a 
UV cutoff with dimension of mass. Therefore, these theories have to be understood as low energy 
effective theories valid for energies below A. This low energy effective action will be our starting 
point in the study of Lorentz- violating massive gravity. But before studying these models in details, 
two remarks are in order: 

• For general functions T of X, and W^^ , we could replace W^^ in the action by a simpler 
variable Y'^ = (9^0* 5^ ^-^ without loss of generality, J" {X, V\ W'^) = P (X, V\ Y'^). The 
reasons why W^^ has been chosen instead of Y^^ will appear in the next sections. 

• Instead of breaking the whole Lorentz group, one requires that the background breaks only 
the invariance under boosts. Then, the function JT has to be understood as a function of the 
following scalar variables 

X, V^V, Tr(VF"), 

where the dots stand for scalar terms which mix and W"^^ (with all indices contracted) while 
n = 1,2,3. This last requirement guarantees the invariance of the action under the rotations 
_^ 0^j(^ of the 0' internal space, which is necessary for the background to be invariant 
under the 50(3) group of three-dimensional rotations. 



27 



Now, let us assume that the following field configuration is the vacuum solution of these theories 

gij.u = Vfj-u, = at, = bx\ (2.16) 

where a and b are two constants. This background configuration obviously breaks the Lorentz sym- 
metry. For this solution, one has that X = a^, = and W^^ = —6^^W/3 with W = 6^. 
Consequently, in this vacuum the first derivative of with respect to is identically zero while the 
first derivative of with respect to W^^ is proportional to 6ij 

The requirement that (2. 16) is a solution of the EoM implies that the energy -momentum tensor of the 
Goldstone fields has to vanish in Minkowski space-time. If J^x is the first derivative of with respect 
to X, this condition translates into two algebraic equations 

= J^-2XTx, = J' + 2WJ^w, (2.17) 

which determine the value of the two constants a and b. Therefore, the vacuum configuration (2.16) 
is a good ansatz and there is a solution for any generic function JT. Actually, once a function T is 
chosen, the two constants a and b could be set to one by a redefinition of the scalar fields. Therefore, 
we could assume that this redefinition has been done and only consider functions for which a = b = 1. 
This will be assumed in what follows, and therefore the vacuum solution will be 

9f.u = Vf.u, cj)'^ = t, cj)' = x\ (2.18) 

Before discussing these models in more details, it is worth noting that the action (2.14) is invariant 
under the usual translations of the coordinates and under shifts of the scalar fields + A"^ by a 

constant vector A". The vacuum configuration (2.18) breaks both invariances. However, for the class 
of models based on (2. 14), the usual invariance of the background under four dimensional translations 
is replaced by translations of the coordinates along with a shift of the scalar fields 

x"^x"-A", (/.°->0" + A". (2.19) 

As a result, there are different possible definitions of the energy for the Goldstone scalar fields [143]. 
For usual matter fields, the gravitational energy is the one which couples to gravity. This energy is 
given by the energy-momentum tensor, and the conservation of this tensor emerges naturally from 
the requirement that the matter fields action be invariant under four-dimensional diffeomorphisms [9]. 
The particle physics energy, or inertial energy, is the Noether charge of the time translation symmetry 
of the vacuum. For usual matter fields, these two definitions are equivalent in Minkowski space-time. 
But things are different for the Goldstone fields. Indeed, their vacuum expectation value breaks the 
time translational symmetry down to symmetry (2.19). Consequently, their gravitational energy is not 
equal to their inertial energy. Indeed, the four currents associated to the four symmetries (2.19) read 

= -2mX + J^, 

where ta is the energy-momentum tensor deduced from 5^, while Jq , Jj*, J2 and J3 are the currents 
associated to the shift symmetries of the four scalar fields. 
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2.2.3 The Goldstone sector 



It is not expected that all models described by an action such as (2.14) will be well-defined theories 
free of pathologies. The study of the Goldstone sector of this action enables to restrain the class of 
viable models by excluding those possessing instabilities [89]. Indeed, as discussed in section 2.1.3 
for Lorentz-invariant massive gravity, this sector describes the k"^ > m? limit of the vector and scalar 
sectors of the perturbation theory about Minkowski space-time. Therefore, any dangerous instability 
present in these two sectors will show up in the Goldstone sector. 
Let us consider perturbations of the vacuum solution (2.18), 

g^,u = V^.u + h^u. <tP = t + e, (t)' = x' + e- (2.20) 

Then, the Goldstone sector is given by all the terms of the Lagrangian quadratic in the Goldstone 
perturbations. After decomposing ^* into transverse and longitudinal fields through (2.3), and for an 
arbitrary function JT, the Goldstone sector of the massive gravity models described by (2.14) reads'"^ 



{di) +(5,eo)' +{ml-ml){dUf + {2ml-ml)i^dU, 

where (n = 0, . . . , 4) are five Lorentz- violating graviton mass parameters. They are combinations 
of the function and its first and second derivatives with an overall factor (see appendix A. 2 for a 
detailed expression of these mass parameters). It is worth noting that the FP theory is obtained when 
rn-Q = and all other masses are equal 

f) O Q O O 

FP : ttiq = 0, rrii = m2 = = m^. 

Here and below, m denotes a parameter of dimension of mass which gives the scale to the Lorentz- 
violating masses 

mj oc with w? = — ^ . 

Ml 

The Goldstone sector does not tell anything about the tensor modes, which will be explicitely dis- 
cussed in the next chapter. However, before going further in the study of the Goldstone sector, one 
should stress that the mass parameter rrig is special in the sense that it corresponds to the mass of the 
helicity ±2 modes. Therefore, one has a first constraint on the function which is that m\ has to be 
strictly positive. Hence, only functions satisfying > will be considered in the rest of this thesis. 

The vector part of the Goldstone sector is given by the quadratic terms in . For m\ > 0, there 
are two healthy propagating degrees of freedom with helicities ±1 in the Goldstone sector. If m\ is 
negative, the vector sector possesses two ghosts while for m\ = there are no helicity ±1 modes. 
Therefore, we will assume that the two following conditions are satisfied 

m2 > 0, > . 

As in the case of Lorentz-invariant massive gravity, the scalar sector is much more tricky. Indeed, 
in the four-momentum space, the Goldstone scalar Lagrangian reads 

>CooM. sc. = f ( kC ) M ( % V M=( ^ ' 



'^The Goldstone sector of Lorentz- violating massive gravity was first studied in Ref. [89]. In this section, we resume the 
arguments presented in that paper. 
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where 




2"! = 2 




/^3 = 




Therefore, the question we would hke to answer here is whether it is possible to find values of the 
three parameters /i, such that the Goldstone scalar Lagrangian is free of classical instabilities and 
ghosts. Let us first assume that 772^ ^ 0. The d.eteriiiina.nt and the trace of the matrix AI are given by 

det M = nioj^ + H2k^ + oj^k'^ (l + fiiH2 -/"!), Tr M = cj^ (/i^ + 1) + fc^ (/xg + 1) . 

The dispersion relations of this sector are given by the equation det M = which has two solutions, 
namely and w^. The product of these two solutions reads 



There are therefore three different classes of models: 

• For /ii / and fi2 7^ 0, there are two propagating modes in the Goldstone scalar sector. 

• For fii = 0, there is at least one non-propagating mode in this sector characterized by k^ = 0. 

• For fi2 = 0, there is at least one non-propagating mode in this sector characterized by = 0. 

For the sake of example, let us discuss the two first class of models. A complete discussion of the 
third class can be found in Ref. [89]. 

Phase /ii 7^ and /i2 7^ 

In this class of models, there are two propagating scalar modes. The absence of classical instabilities 
implies among others that 



As discussed in section 2.1.3, the requirement that the scalar sector is free of ghost implies that each 
eigenvalue of M (given by eq. (2.8)) has to have one zero. This condition is not satisfied if fii, fi2 > 0, 
since then the trace of M is positive-definite. This implies that the equation A+ = has no roots while 
the equation A_ = has two roots, one corresponding to a ghost. 

If /Ui, /i2 < 0, the Goldstone sector also possesses a ghost. Indeed, if we assume that has only 
one zero, this mode correspond to a ghost since det M < for large value of lo'^ implying that A_ is 
negative at large w^. Therefore, the condition (2.9) cannot be satisfied. 

Phase /Ui = 

For fii = 0, the solutions to the equation det M = reads 



22 P2 , 4 

',uj_ = — k . 





This class of solutions reduce to three subgroups. 
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• For /i2 7^ and ^| 7^ 1, there is one propagating mode in the Goldstone scalar sector: 

- If fi2 and /ig — 1 are both positive, the trace of M is positive definite, implying that this 
propagating mode is a zero of A_. Since det M < and therefore A_ < for large a;^, 
one has to conclude that this propagating mode is a ghost; it does not verify the condition 
(2.9). 

- If fi2 and — 1 are both negative, the trace and determinant of AI are both positive at 
large w^, implying that both eigenvalues are also positives in this limit. Therefore, the 
derivative of the eigenvalue which has a zero is positive in its zero, and the propagating 
mode is perfectly healthy. 

• For /U2 = and /ig 7^ 1, there are two non-propagating modes in the Goldstone scalar La- 
grangian, with dispersion relation = and lo'^ = 0. 

• For /U2 7^ and = 1, there are two non-propagating modes in the Goldstone scalar La- 
grangian with dispersion relation /c^ = 0. 

Finally, there is a possibility that both fi2 and /i| — 1 vanish. This will be equivalent to the situation 
of the FP theory, where det M is identically zero and nothing could be said from the analysis of the 
Goldstone scalar sector. 

Phase ml = 

The three previous class of models assume that ttiI 7^ 0. Therefore, there is a fourth class of models 
characterized by mf = 0. If this last relation holds, the determinant of the matrix M reduces to 

det M = 4 [ml [mj - mf) - mj] k^u?. 

Consequently, if ttiq (m| — m^) — m\ 7^ 0, there are two non-propagating modes with dispersion 
relations /c^ = and = 0. If this relation between the parameters is not satisfied, then det M is 
identically zero and the Goldstone sector does not tell anything about the scalar sector. 

To conclude this discussion of the Goldstone sector of Lorentz- violating massive gravity theories, 
let us stress that models with m\ = have only two propagating massive tensor modes, and can be 
thought of as minimal models of massive gravity. The exact linearized approximation of these models 
will be discussed in the next chapter. 

2.2.4 Symmetry as stabUizer 

Imposing certain relations between the five mass parameters enables to get rid of ghosts scalar degrees 
of freedom. But this is not sufficient to make the theory healthy. Indeed, as discussed previously 
for Lorentz-invariant massive gravity, the absence of ghost is guaranteed by the fine-tuning relation 
a = —1. However, instabilities may reappear in non-fiat background since curvature terms could 
be responsible for a violation of this fine-tuning relation. Likewise, the stability of Lorentz-violating 
massive gravity in Minkowski space-time is also protected by fine-tuning relations, implying that 
Boulware-Deser instabilities may also be present in those models. Moreover, instabilities may appear 
in the spectrum upon the addition of higher order operators. 



31 



To illustrate this problem, let us once again consider the Goldstone scalar sector. In a nearly flat 
arbitrary background, the Lagrangian of this sector reads 

>CGolds. sc. = {ml + ...)i^ + ^{mj + ...) (ditj \^{ml + ...) {di^of 

+ {ml-ml + ...) {dftf + {2ml - mj + . . .) ^odf^ + ..., 

where the dots stand for curvature terms. If the five mass parameters are non zero, it has been demon- 
strated in section 2.2.3 that such Lagrangian describes two propagating degrees of freedom, one of 
them being a ghost like mode. Yet, it is possible to get rid of this dangerous scalar mode by imposing 
for example that there are no quadratic terms for with two time-derivatives. This last requirement 
reads 

ml + ... = 0. 

This condition is easily satisfied in Minkowski space-time if the function is chosen such that ttiq = 
0. Still, this condition seems impossible to satisfy for arbitrary curved space-time. Indeed, let us take 
for example the following function 

Then, the vacuum configuration (2.18) with a = 6 = 1 is a solution of the EoM with rn^ = m\ = 0. 
According to the previous discussion of the Goldstone Lagrangian in flat space-time, this particular 
model possesses four massive modes with helicities ±2 and ±1, respectively. The two scalar modes 
are non-dynamical with a dispersion relation A;^ = implying that the vacuum is perfectly stable. 
However, this theory is unstable for the following background configuration'^ 

4>^ = t + il^{t), 4>' = x' + ii^, 

where ^bg(i) and ,^^g are solutions of the linearized EoM above Minkowski space-time. Indeed, for 
such a background one has the following terms in the Goldstone scalar sector 

(m2-12m29oeOgW9°° + ---)4VO. 

Therefore, there are quadratic terms for with two time-derivatives in the Goldstone Lagrangian and 
this model possesses a ghost-like mode in curved space-time. Fortunately, there is an elegant solution 
to such problems. 

In several cases, requiring that these fine-tuning relations are a consequence of unbroken gauge 
symmetry is sufficient to protect Lorentz-violating models of massive gravity from becoming patho- 
logical upon addition of curvature terms and/or higher order operators [89]. Indeed, if a symmetry 
forbids one particular term, this term will reappear neither in curvature terms nor in the UV com- 
pletion of the theory. Therefore, if the absence of this term is sufficient to guarantee the stability 
of the model under consideration, then this model will be UV stable and free of instabilities of the 
Boulware-Deser type. An example of theory whose stability is insensitive to its UV completion is 
given by the Ghost Condensate. This model is required to be invariant under space-time dependent 
spatial reparametrization 

x' ^ x' + cf {t,x') , (2.21) 

"The Friedmann-Lemaitre-Robertson- Walker (FLRW) background will correspond to this configuration of the scalar 
field in the limit of small curvature. 
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implying that all masses but mg have to vanish. There are other symmetries one can think of. For 
example, the symmetry 

t ^ t + ?° (t, x') , (2.22) 
guarantees that niQ = mi = m4 = while the symmetry 

^x' + (t) , (2.23) 

implies that mi has to be zero. 

It is worth stressing that precautions have to be taken when considering models where some of the 
fields are non-dynamical; they may become dynamically unstable after taking into account higher- 
order terms and/or curved background. An example of such behavior is given by models invariant 
under the symmetry (2.22), which become dynamically unstable after accounting for higher order 
operators [28]. Let us discussed the two other symmetries separately. 



2.2.5 The symmetry a;* — * a;* + (t, x*) and the Ghost Condensate 

The class of models invariant under the symmetry (2.21) are known as Ghost Condensate models 
[48]. Because of the relation between the four scalar vacuum expectation values and the four coordi- 
nates, the requirement that the theory is invariant under this symmetry independently of the gauge is 
equivalent to requiring that the full action is invariant under the following transformation of the scalar 
fields 

</.^^</.' + x'(</'°, (2.24) 

where are arbitrary functions of cf)^ and (p^. Therefore, the action of Ghost Condensate models 
depend only on the field (j)^ and reads 

S^ = A^ J d^x ^gT (X) . (2.25) 

The Minkowski vacuum is a solution of the EoM for = t\i the function T is chosen such that the 
following two relations hold 

T\x=\ = and Tx\x=\ = 0. 
At the quadratic level, the Lagrangian of the Ghost Condensate simply reads 

= + ^ (/ioo-2eo)'. 

Therefore, the new term will not affect the tensor and vector sectors of the gravitational field; i.e., the 
tensor modes are massless and there is no vector degrees of freedom. It is then straightforward to see 
that there is no scalar propagating mode either, as the dispersion relation deduced from the Goldstone 
sector is simply uP' = 0. Indeed, the field equations for the scalar sector are 

= 2dl^ + - Ho) , = ^' + - ^) , = (6 - Hq) , 

whose solutions read 

* = ^0 {x') , ci> = ^0 , Ho = (^1 + M/Q [x^) , 
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where ^'o (x*) is a time-independent integration constant. Note that if ^I'o = 0, this solution reduces 
to the solution of GR. Ghost Condensate models could therefore be seen as Einstein gravity plus an 
initial condition on the potentials. 

The UV stability and the absence of dangerous Boulware-Deser instabilities of Ghost Conden- 
sate models is easily understood by studying the Goldstone sector. The absence of gradient terms 
in this sector signals that the non-propagating modes will become dynamical upon addition of cur- 
vature terms or higher order operators. Let us first concentrate on the UV completion, which brings 
terms with more than one derivative acting on cfp. Because of the symmetry (2.24), the most general 
quadratic Goldstone Lagrangian for the perturbation is given, in the four-momentum space and up 
to fourth order operators, by 



2/2 k'^^'^ \ 

£ = mo ( w - co^ - ci^^ - + . . . ) eo^o 



where Cj are constants of order one. These higher order terms are suppressed by the cutoff scale = 
mMpi and are therefore small corrections to the low energy effective action (2.25). The dispersion 
relations deduced from this Lagrangian are 

"^ = --^f'^^-2X^+oy , '^^ = «p + 0(^). (2.26) 

For the low energy effective theory with ^ A^, the first relation is irrelevant. The second disper- 
sion relation is the one of interest. It shows how the non-propagating mode with = becomes a 
slowly varying dynamical mode upon addition of higher order terms. One could therefore conclude 
that imposing tuq > is sufficient to guarantee that the Ghost Condensate will be free of ghosts after 
considering higher order terms. Note that C2 has to be positive otherwise this modes will be a slowly 
varying classical instability. 

Therefore, Ghost Condensate models are truly healthy in the absence of gravity in contrast with 
the analogous sector of the DGP model which only becomes healthy trough its coupling to gravity 
[124]. When coupled to gravity, the scalar sector is classically unstable at low momenta. Indeed, 
after taking into account the higher order terms which modify the relation = 0, the full scalar 
Lagrangian reads in four-dimensional Fourier space 



:0* 



—Ik"^ niQ iniQU) 

y -imluj ml ^io'^ - J j 




The determinant of this matrix shows that two of the three scalar fields are non dynamical with the 
dispersion relation A;^ = while the third field satisfies the relation 

This last relation shows that when coupled to gravity, the ghost condensate is classically unstable for 
small momenta 2k^ <^ ml- Such low-momentum instabilities are nothing else than the analogs of the 
Jeans instability for ordinary matter fluid. 

Similarly, Ghost Condensate models do not exhibit dangerous Boulware-Deser instabilities, con- 
trary to Lorentz-invariant massive gravity. Indeed, the propagating mode will never be a ghost since 
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curvature terms cannot change the sign of the leading term proportional to up' . This is completely 
different from the situation encountered in the FP theory discussed in section 2.1.7. In that model it 
is a non-propagating ghost which becomes dynamical in slightly curved space-time. Still, in slightly 
curved background, the dispersion relation (2.27) could acquire additional terms with small coeffi- 
cient controlled by the curvature. Those terms could for example be gradient terms with a negative 
sign — fc^ implying a Jeans-like instability at low spatial momenta. However, such low momenta clas- 
sical instabilities are not dangerous and are even interesting from a phenomenological point of view. 
Indeed, in GR such instabilities are believed to be responsible for the formation of the structures of 
the Universe. 

Ghost Condensate models are interesting as toy models which enable one to understand the fea- 
tures of theories with spontaneous breaking of Lorentz-invariance. Moreover, some properties of 
those models leads to an interesting phenomenology. For example, the Newtonian potential is modi- 
fied at distances larger than ttiq ^. However, this modification only shows up after a time tc = A/ttiq 
[48], since it is a consequence of the fc^ term which is suppressed by the UV cutoff A (recall that, 
without the UV completion there is no modification of the scalar potentials). This modification of the 
Newtonian potential could be interpreted as resulting from an effective mass term produced by lumps 
of the scalar field excitation ^o- Such lumps could have a positive or negative effective mass leading 
to an attractive or repulsive contribution to the Newtonian potential. 

Since the retardation effects are very strong in the Goldstone sector, the time delay between the 
passage of an object at a given point an the appearance of its track in the potential at this point is of 
order t^. Consequently, the Ghost Condensate could be thought of as a bubble chamber in which all 
moving massive objects leave tracks [144, 145], and observer in motion will not have time to see these 
new effects in the potential while an observer at rest will see them after a time t^. 

Due to the breaking of Lorentz-invariance, the non-linear dynamics of Ghost Condensate models 
is quite rich. For example. Ghost Condensate tends to form caustics [143]. The Ghost Condensate 
seems also to be non-perturbatively unstable. This leads to the formation of microscopic bubbles 
with negative energy density in their centers [146]. However, those "holes"of negative energy do not 
expand and have little effect on the outer part of the bubbles. 

From a cosmological point of view. Ghost Condensate could be a dark matter candidate [143, 147]. 
Moreover, there is a term proportional to T in the EoM which characterizes a fluid with equation of 
state p = —p. This fluid could give rise to a de Sitter expansion phase of the Universe without need 
for a cosmological constant Ac [48]. It has also been shown that Ghost Condensate could drive an 
inflationary de Sitter phase of expansion if the shift symmetry is broken by a potential V[^^^ [148]. 

To conclude this discussion of the Ghost Condensate, let us mention another very interesting prop- 
erty of Ghost Condensate models concerning black holes physics. Due to the higher order terms, 
there is an accretion of the Ghost Condensate onto black holes with a modification of the usual 
Schwarzschild black hole solution [149, 150]. Moreover, black holes may have hair [98] as informa- 
tions may escape frome black holes through a possible violation of the usual laws of thermodynamics 



2.2.6 The symmetry ^ + it) and a minimal model of massive gravity 

Similarly to the situation of the Ghost Condensate model, massive gravity theories invariant under 
the transformation — > + C {t) are perfectly healthy [89, 125]. In the Goldstone picture, this 
symmetry reads 



[151]. 




(2.28) 
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where are arbitrary functions of (fp. This transformation is a symmetry of the action if it has the 
following form 

S = j d^x [-M2^7^ + k'^J^ [X, W'^)] . (2.29) 

Indeed, both X and W^^ are invariant under this transformation while is not. This last observation 
motivates the fact that up to now W^^ has been used as a variable of while it would have seem 
simpler to use Y^K At the level of the quadratic action, the requirement that does not depend on V" 
implies that m\ = 0. Therefore, according to the discussions of section 2.2.3, in Minkowski space- 
time there are only two modes with helicities ±2. But contrary to Ghost Condensate models, these 
modes are now massive with consequences that (2.29) leads to true modifications of the gravitational 
interaction. 

To address the issue of the UV stabihty of these models, let us once again concentrate on the Gold- 
stone sector. Under the assumption that higher order operators are also invariant under the residual 
gauge symmetry ^ + {t), the UV completion cannot contain (^q?*)^ and is therefore given 
up to fourth order derivatives by 

where the q with z = 0, . . . , 6 are constants of order one. Let us first have a look at the vector sector. 
The Goldstone vector sector in the Fourier space reads 

r 2 ^'^ 2 C4 2\ \ i-T\2 
^Golds. vect. = "i2C3-r^ "7^ k |Ci | ■ 



A2 V 2C3 C3 

One conclude from this Lagrangian that upon addition of higher order derivatives, the non-propagating 
vector modes already present at the classical level remain non-dynamical. Indeed, the dispersion 
relation fc^ = is not affected by the UV completion. Moreover, there is a second vector mode 
which falls out of the validity region of the low energy effective theory since it describes a mode with 
frequency of the order of A. Therefore, the vector sector is not sensitive to its UV completion. 

The Goldstone scalar sector up to fourth order derivatives is much more tricky. For simplicity, let 
us neglect the contributions proportional to cq and ci Then, the field equations for the scalar sector 
are given in the four momentum space by 




-iojkml ( 1 + C6 



A2 



k'^ [ ml - ml + ^ (c3a;2 + (cs - C4) k'^ ^ j j 



^0 
ki 



0. 



As already discussed, the dispersion relation are simply given by the equation det M = 0, where M 
is the matrix of the EoM. It is straightforward to see that these dispersion relations are 

I 1 + , , "^^ ^ I ^ C2^ + . . . , ^ = k^{u? - const. A^ - const, /c^ + . . .) . 

\ rriQ \m\ — m|j / A^ ^ 



'**Recall that we have seen in the discussion of the Ghost Condensate that these terms give rise to a mode (2.26) which is 
unimportant for the low energy effective theory. 
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The first of these two relations is the generalization of the relation (2.26) already met in the study 
of the Ghost Condensate. Indeed, this first relation precisely describes the Ghost Condensate mode. 
Therefore, the Ghost Condensate is somehow included in these massive gravity models. The condition 
rriQ > that guarantees that this mode is not a ghost-like degree of freedom in the Ghost Condensate 
models is replaced here by 

4 

ml + „ > 0. (2.30) 

TTig — mg 

It is sufficient to impose this condition for the lj^ term present in the classical two-derivative theory 
to have the correct sign. As for the Ghost Condensate models, this mode is a slowly varying classical 
instability if C2 < 0. 

Similarly to the situation in the vector sector, the second scalar dispersion relation fc^ = is in- 
sensitive to the UV completion. This mode is still a solution of det M = after taking into account 
higher-derivatives terms. This is a direct consequence of the residual symmetry which forbids terms 
quadratic in ^ without at least two spatial derivatives. Similarly to the vector sector, there is a new 
mode uj^ = const. +const. fc^ + . . . which falls out of the low energy effective theory. Consequently, 
we have just demonstrated that massive gravity theories based on (2.29) are perfectly healthy upon 
addition of higher order operators. 

Those models do not exhibit dangerous Boulware-Deser instabilities either, since curvature terms 
cannot change the sign of the term proportional to or introduce such terms for the ^ fields. Hence, 
there is no ghost propagating mode in curved space-time contrarily to the situation of the FP theory. 
Still, in slightly curved background, the dispersion relation o;^ = could acquire additional terms 
with small coefficient controlled by the curvature. Those terms could for example be gradient terms 
with a negative sign —k"^ implying a Jeans-like instability at low spatial momenta. 

At this stage, it is clear that these models are candidates for consistent massive gravity theories. 
Consequently, from now on we will only focus on these models which are very interesting from a 
phenomenological point of view even without taking into account their UV completions. 
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3 

WealijravitationaC interaction 



GR is a full non-linear theory, with consequence that only a few exact solution of the Einstein's field 
equations are known. Fortunately, the gravitational interaction is so weak that for most practical situ- 
ations it is sufficient to study the solution of the linearized EoM. For instance, in GR the gravitational 
field of the Sun is static to a very good approximation implying that the metric around the Sun is given 
by the Schwarzschild solution 

ds^ = (1 + 2^>) df - (1 + 2^>)"^ dr^ - [dO"^ + sin^ Odip^) , $ = 

Thus, for distances much larger than the Schwarzschild radius = 2GMq ~ 3km, the potential 
$ is a small quantity. Consequently, the Schwarzschild metric can be approximated by its Taylor 
expansion around rg = 0, and all terms of order {rg/r)^ with n > 2 are completely negligible as 
compared to the linear terms 

ds2 = (l - ^) dt2 - (l + ^) dr2 - r2 (d^^ + sin^ Od^^) + .... 

Therefore, the gravitational field of the Sun is correctly described by the linearized approximation for 
distances r » 3km. 

In this chapter, we discuss the perturbations above Minkowski space-time induced by a small- 
amplitude gravitational source ST^j^y in Lorentz-violating models of massive gravity described by ac- 
tion (2.29). The linearized theory which results from this corresponds to the weak field limit of the 
gravitational interaction produced by 5T^y as long as these perturbations remain small as compared 
to the metric of flat space-time. This chapter has two aims. The first one is to introduce a sub-class 
of minimal models for which there is no modifications of Newton's potential. The second one is to 
provides the equations which will be the starting point of chapter 4. 

3.1 Energy-momentum tensor 

Before studying the linearized approximation of the massive gravitational field, it is worth noting that 
there is a convenient way to parameterize the energy-momentum tensor of any gravitational source. 
Let with v^j^v^ = 1 be the velocity of an observer comoving with a fluid described by the energy- 
momentum tensor 7^j^. The two following tensors are projection tensors along and perpendicular to 
[152], respectively, 

satisfying 

J — KJ, — o, J jau — Ju- 
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These two projection operators enable one to introduce natural notions of time and space for the 
observer since the space-time interval can be written 

With these projectors, any energy-momentum tensor can be decomposed as [8, 152] 

T^u = pv^Vy + q^Vu + qyV^ + pj^^ + -k^u, (3. 1) 

where p = T^yV^^v^ and p = T^yj^'^ /"i are the energy and pressure densities measured by the co- 
moving observer respectively, q^^ = —Tyxv^p^ is the energy flux perpendicular to v^^ and vr'^'^ = 
'^x-fj^^P'^ —j^'^'^X'yP^ /3 is the anisotropic pressure tensor. The following relations are consequences 
of the previous definitions 

3.2 Weak field 

In order to determine the weak field approximation in Lorentz- violating massive gravity models, 
one has to consider small perturbations about the flat vacuum solution (2.18) for which the energy- 
momentum tensor of the usual matter is zero, T^y = 0. A small perturbation 1 6Tfj_y | ^ 1 of the matter 
energy-momentum tensor will produce metric 6g^y = h^y and Goldstone perturbations 5<j)'^ = 
about the vacuum solution. 

With the notations introduced previously, the energy-momentum tensor of a small-amplitude 
source of gravitational field above Minkowski space-time is parameterized as follows, 

^T^v = {5p + 5p) Vf^uSp + {v^5qy + Vy5qfj) + 6-K^y, 

where Sp and Sp are the matter density and pressure measured by the comoving observer, Sq^ is the 
energy flux perpendicular to {v^^q^j. = 0) and Sir^j^y is the anisotropic pressure tensor {v'^dir^y = 
^Trjl = 0). The velocity of the observer obeys the geodesic equation v^dyV^^ = 0, with f'^w^ = 1. 
This implies that the affine parameter of the observer can be chosen such that = (1, 0, 0, 0). As a 
consequence, Sq^ = diTQy = and 

5Tqq = 5p, 5%i = 6qi, 5Tij = 6ij5p + 5-nij. 

The energy flux 5qi and the anisotropic stress di^ij can be parameterized in the following way, 

5qi = Ci + diC, 

STTij = {3didj -6ijdl)7r + diTrj + djTri + 7rij, 

where the vector perturbations Q and vTj are transverse while the tensor perturbation vTjj is transverse 
and traceless ' . With these notations, the energy-momentum conservation reads 

S'p = dfC, C = Sp + 2df7r, Ci = d]7ri. (3.2) 

Hence, the ten independent components of the linearized energy-momentum tensor are expressed 
through four gauge-invariant scalars 6p, 6p, ( and vr, four vector degrees of freedom in the form 

'For illustration, an imperfect fluid with shear viscosity rj will have a transverse and traceless anisotropic stress given by 
^ij = —vHij [8]. 
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of two transverse gauge-invariant vectors C,i and vTj, and two tensor degrees of freedom through the 
transverse and traceless gauge-invariant tensor vTjj. 

As mention earlier, the presence of the source 5T^u implies a modification of the solution (2.18) 
for the metric and Goldstone fields. The metric and Goldstone perturbations are described through the 
decomposition (2.3) into transverse and longitudinal fields. With the account of all these notations, 
and because of the invariance of the vacuum solution under three-dimensional rotations, the linearized 
Einstein equations split into scalar, vector and tensor equations (appendix A.3.3). 



Mi 



liu) 



linearization 



tensor equations 



vector equations 



scalar equations 



Where Q^y is the Einstein tensor and where t^y is the energy-momentum tensor of the four Goldstone 
fields deduced from (see appendix A for a complete expression of this tensor). This splitting has 
already been observed in the previous chapter, where the quadratic Lagrangian of the gravitational 
field decomposes into scalar, vector and tensor sectors. The equations of these three sectors could 
then be consider separately. The discussion presented here was first made in Ref. [97]. 

3.2.1 The tensor and vector sectors 

The tensor linear equations states that any transverse and traceless anisotropic pressure tensor vrjj will 
generate massive gravitational waves 



= (□ + ml) Hij + 



(3.3) 



The solution to this relation can be formulated in the following way 



Hi 



M 



pi ° + "^2 



2 ' 



(3.4) 



where (□ + ml) ^ has to be understood as the Green function of the Klein-Gordon equation with 
mass m2. In this equation, we neglect the tensor modes which are solutions to the homogeneous 
equation (□ + m^) Hij = 0. These modes are present in the vacuum and do not play any role in the 
interaction between gravitational sources. 

There are two vector equations deduced from the linearized EoM 







djVJi 



2Q_ 







Wi 



m^ai 



M 



(3.5) 



pi 



One concludes from these two equations that there is no difference in the vector sector between GR 
and the model of massive gravity given by (2.29). Indeed, the first equation is exactly the equation 
deduced in GR while the second equation reduces to m^ai = 0. Then, using the energy-momentum 
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conservation, the equation for Wi can be written as = 2Mpj^ (tTj — C,i) implying that the solutions 
to the vector equations are given by 

= 2 " ' 2 ' = 0, (3.6) 

since we are interested in models with massive gravitons, i.e. > To conclude the discussion 
of the tensor and vector sectors, it is worth stressing that there is no vDVZ discontinuity in these two 
sectors; the solutions to the tensorial and vectorial equations both correspond to the solutions of GR 
in the limit of vanishing mass m? 0. 

3.2.2 The scalar sector 

As usual, the scalar sector is more complicated. There are four scalar equations which read 

= -2df^ + yi-^5p + ml - + (S^H + 3^) , (3.7) 

= -2^ + Mpi25^, (3.8) 

= ^> - ^' + m^H + SMpi^, (3.9) 

= 2^ + df{^--^)- Mpi^ {5p - dfn) + (3^ + dfE) - 

-ml($-S°j. (3.10) 

By making use of the conservation relations (3.2), eq. (3.8) implies that 

The first term in this relation is the usual GR contribution to the potential ^'o is a function of the 
space coordinates only which corresponds to a new contribution to this potential. Combining the three 
other equations one finds that Newton's potential ^ is given by 

^ = ^U^-3- +4:-^+ l---4 + ^^7#Wo(xO, (3.12) 



Ml,\2df J 2p Ml^df V f^rnl /i 5^ 
with 

ttiq — ml) + ml 

2 2 ' 

m2mQ 

The first contribution to the potential is the usual GR term while the two other contributions are new 
terms absent in GR. Since ^ is gauge-invariant, the previous relation is valid provided that <I> ^ 1. 
Therefore, the second contribution proportional to is responsible for a breakdown of perturbation 
theory at large distances. Indeed, the gravitational potential of a static source 6T^,y = 6o^6oi,M6^ (x) 
has the form 

^ ^.rf 1 2-3^ 2\ (. 2m? 3;U-2mi\ . 



"For models characterized by — 0, gravitational waves are massless and there are no constraints on the value of Ui. 
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where G = (^SvrMpj j is Newton's constant. Since the second contribution to this potential grows 
linearly, the perturbation theory breaks down when this term reaches one, that is, at distances 



2^ 



GMml2 



3fi 



< 



As an example, considering the Sun and assuming a graviton mass of the order m ~ (lO^^cm) 
(which correspond more or less to the size of the Solar System), this last expression tell us that the 
linearized approximation breaks at distances r > {GMm?^ ^ ~ lO^^cm which is beyond the actual 
tests of the Solar gravitational field. For a smaller mass of the order of the Hubble scale m ~ Hq^, 
the linearized approximation breaks at distances much larger than the actual Hubble scale. Therefore, 
this breakdown of perturbation theory is not problematic from a phenomenological point of view 
contrarily to the breakdown of the FP theory who occurs at relative small distances (section 2.1.5). 
Still, this growing term puts strong constrains on the graviton mass (2 — 3/u) m^j^ix <^ r^^ since 
Newton's potential is verified to a very good approximation in the Solar System [6]. 

With the previous solutions for the two gravitational potentials 4* and <1>, one can invert eqs. (3.9) 
and (3.7) to obtain the following expressions for E and H", 



2 — 3/i Sp 
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(3.13) 
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(3.14) 



These last two equations complete the analysis of the scalar sector of the massive gravitational field in 
the presence of a sources ST^y. Note that all these relations imply that while $ and H are static for 
a static source, rP is growing linearly with time. Therefore, for a static source, the solution found here 
is the approximation of a stationary solution for which there is an accretion of the Ghost Condensate 
fluid rP by 5T^,. 

The previous solution for the scalar sector contains the time-independant function ^'q. This func- 
tion corresponds to a mode with dispersion relation = (section 2.2.3) analogous to the Ghost 
Condensate mode which should become dynamical upon addition of higher-derivative terms. Yet, 
this function is fixed at the two-derivative level by the initial conditions on the gravitational system. 
It follows from eqs. (3.13) and (3.14) that a non-zero ^'o would mean the presence of a non-trivial 
condensation of the scalar fields (j)^. This function does not depend on STf^y and is therefore present 
in the solutions of the scalar sector even without matter fields. Hence, this function does not play 
any role in the interaction between sources of gravitational field and we will disregard its contribution 
when studying the interaction between gravitational sources in the rest of this thesis. 

Contrarily to the situation of FP theory, there is no vDVZ discontinuity in the scalar sector of mas- 
sive gravity models described by (2.29). Indeed, mf oc m^ with the consequence that the parameters 
p is of order one even in the limit m ^ and relation (3.13) implies then that ^'o oc m^H ^ in the 
limit of vanishing mass. Consequently, both \I' and $ tend to their values in GR. 
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3.3 Dilatation symmetry 



For particular models satisfying 3^ = 2, the term responsible for the breakdown of the perturbation 
theory disappears and Newton's potential reduces to its value in GR (plus the contribution of ^'o). 
Therefore, those models are very interesting from a phenomenological point of view. They describe 
massive gravitons without any modifications in the vector and scalar sectors. A way to achieve this is 
to impose the following dilatation symmetry as a symmetry of the Goldstone action [97] 



(3.15) 



This dilatation is a symmetry of massive gravity theories if their action depends on the variables X 
and W^^ through a new variable Z*-' = X'^W^^ 



(3.16) 



For the vacuum configuration (2.18) one has Z*-^ = —6"^^. For models based on a function of a 
single variable Z*-^, there are two relations between the mass parameters (see appendix A.2) 
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which imply that ?>n = 2. Consequently, for those models the scalar potentials read 



1 



while the Goldstone perturbations are given by 
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(3.17) 
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(3.18) 



Note however that in the absence of anisotropic pressure vr = 0, GR predicts $ = This relation 
is not satisfied in the massive gravity model because of the appearance of the function ^'q. Hence, if 
^'o 7^ the gravitational potential ^' created by a massive source is not the potential ^> responsible 
for the geodesic motion of an observer around this source'. If we neglect this arbitrary function there 
is no difference in the scalar sector between GR and the models given by (3.16) at the linear level. 

For such models, there are only a few constrains on the graviton mass coming from direct or 
indirect observations of gravitational waves. For instance, the secular decrease of orbital period of 
binary pulsar puts bounds on the graviton mass 1712; the slow-down of orbital period of binary pulsars 
has been shown to be compatible with the emission of gravitational waves as predicted by GR [15, 
153]^. This implies that the mass of the gravitational waves cannot be larger than the frequency of the 
waves emitted by these system [14] 



1712 

2^ 



< w ~ 3.6 10"^ Hz ~ 



flO^^cm) 



10"^^ GeV. 



This frequency corresponds to the period of the orbital motion which is of order of 8 hours. Apart from 
the relatively high graviton mass allowed in such models, there is another motivation for studying the 

^The Cassini spacecraft has put a bound on the differences between those two potentials '5/$ — 1 < 10~^ [2]. 
'^For a recent discussion on bounding the graviton mass with binary pulsars, see [154]. 
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sub-class of models (3.16). Indeed, we will see in chapter 6 that those models are attractors of the 
cosmological expansion [125]. 

To conclude, let us stress that the only differences between GR and massive gravity models de- 
scribed by (3.16) lei in the tensor sector which possesses two massive modes. This explain why these 
models have been studied intensively, and why most of the original contributions of this thesis will 
concern them. Before closing this chapter, let us note that the solutions to the linearized field equations 
which have been discussed here are the starting point of the coming study of physical instantaneous 
interactions. Moreover, these equations and their solutions will be generalized as to describe pertur- 
bations about a flat Friedmann-Lemaitre-Robertson-WaUcer space-time in chapter 6, in order to study 
the formation of structures in massive gravity theories. 
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Part II 

S^pects of massive gravity physics 
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4 

Instantamous interaction 



One peculiarity of massive gravity models based on the action (2.29) consist in the presence of a 
physical instantaneous interaction [155]. In the conventional GR, the gravitational potentials are all 
instantaneous potentials. Still, there is no physical instantaneous interaction since the instantaneous 
contributions cancel in the graviton propagator. An analogous situation exists in classical electro- 
dynamics, where the instantaneous contributions to the potentials Aq and Ai cancel each other in 
the photon propagator, leaving the theory free of instantaneous interactions [156]. In models where 
these subtle cancellations are spoiled by the presence of fields that break Lorentz-invariance, physical 
instantaneous interactions are possible. An example of such models is given by Lorentz-violating 
electrodynamics [157, 158]. 

Before considering the gravitational field, we review the Lorentz-violating electrodynamic model 
as a warming up. This model shared all the relevant features of massive gravity theories. The emer- 
gence of the electrodynamic instantaneous interaction in that model is analogous to the emergence 
of the instantaneous interaction in our class of massive gravity theories. Then, we will discuss the 
physical instantaneous interaction of massive gravity and give a concrete example consisting in a in- 
stantaneous frequency shift of light by a distant gravitational source. This chapter constitute the first 
original part of this thesis. 

4.1 Lorentz-violating massive electrodynamics 

As a toy model for physical instantaneous interaction, consider the Lorentz-violating electrodynamics 
of Ref. [157, 158], described by the following Lagrangian 

£ = -\f,,F^^ + \mlAl - \mlAl + A,J^, 

where is a conserved current 9^ J'^ = 0. The Lorentz-invariant Proca theory of massive electrody- 
namics is recovered by imposing mo = mi. Note that the gauge invariance can be restored by making 
use of the Stiickelberg trick. Hence, this action may be thought of as the unitary gauge description of 
Lorentz-violating massive electrodynamics. 

This Lagrangian is still invariant under spatial rotations implying that the following parametriza- 
tion of the spatial vectors Ai and Ji will simplify the equations of motion: 

Ai = Ai + diA , Ji = Ji -\- diJ , 

where and jf are transverse. It is important to notice that the transverse and longitudinal currents 
jf and diJ^ extend over the whole space, even if Jj is localized [159]. Indeed, for all localized 
currents Jj, the transverse current jf possesses an delocalized contribution which is the opposite of 
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the delocalized contribution of diJ^. As a concrete example one can think of a dipole which appears 
at t = on X* = 

Jo = -A^iiiQ {t) 9,53 (x) , Ji = -ATTfi,6^ (x) . 

Despite the fact that this current is conserved and perfectly locaUzed, its longitudinal and transverse 
components extend over the whole space 



3x^x^ 8tt5. 



6-'{x) 



3x^x^ 



6ij6^ (x) 



The transverse and longitudinal currents jj" and diJ^ associated to the localized current Jj can be 
determined using the transverse and longitudinal projection operator introduced in [156] (see appendix 
B. 1). Therefore, for any localized current appearing in a small region of space there are one transverse 
and one longitudinal currents which extend to the entire space. 

For an arbitrary current J^, the EoM deduced from the previous Lagrangian are 



{Ao - A^) 



- m^oAo - Jo 



do (A 



J' 



0, 

0, 
0. 



(□ + ml) Af + J, 

The potential Ao — A^ is invariant under the U{1) gauge transformations, and for mo = mi = 0, 
these equations are precisely those of classical electrodynamics. The first two equations imply the 
generalized Proca constraint mQ^o = m\dfA^. Let us concentrate on the model defined by mi > 
and mi ^ mo- Then, mo has to be zero to guarantee the absence of ghosts and rapid instabilities. 
Since we are interested in solutions which do not grow at infinity, the Proca constraint implies then 
that A^ = 0. Hence, the field equations are equivalent to 

Jo 

dV 



Ao 



[U + mf) Aj = - J: 



If Jq 7^ 0, it is obvious from these equations that Ao is an instantaneous potential, since Jq and Ao 
have the same time dependence. In fact, is also an instantaneous potential. Indeed, the transverse 
current ^ extends over all space and implies a non- vanishing in all space. As a consequence, for 
all localized currents Ji, Af is an instantaneous potential '. 

Using the current conservation, the equation for ^0 can be written OAo = J^ — Jo- Hence, the 
solutions to the field equations are 



A 



dVG+ {t,x,t',x') {t',x') - Jo {t',x) 
j dVG+ (t,x,t',x') Jf (t',x'), 



where {t,x,t',x') and {t,x,t',x') are the retarded Green functions of the d'Alembert and 
Klein-Gordon equations 

aG+ {t,x,t',x') = 6'^{x-x'), 
(□ + m^) G+ (t,x,t',x') = 6^{x-x'), 



'There is only one possibility for Aj and Ao not to be instantaneous potentials ; A[ will not be an instantaneous 
potential if the localized source is transverse Ji — jf while Aq will not be an instantaneous potential if the localized 
source is longitudinal Ji — diJ^ . 
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respectively. 



In classical electrodynamics, the instantaneous contributions to the potentials Aq and Aj cancel 
each other in observables such as the electric and magnetic fields [156]. There is then no instantaneous 
interaction. This is no more true in the Lorentz-violating model, because the mass mi modifies 
the dispersion relation of the transverse modes Af without affecting the dispersion relation of Aq. 
Therefore, the instantaneous contributions to the potentials do not cancel each other anymore, giving 
rise to a physical instantaneous interaction. This can be illustrated by looking at the electric field 
Ei = Foi, which reads 



where AG+ {t, x, t', x') = {t, x, t', x') — {t, x, t', x') is proportional to mi. The first integral 
in (4.1) is the retarded electric field of classical electrodynamics. This term vanishes outside of the 
light-cone of the source: it involves only localized current J^, whose effects are retarded because 
of the retarded Green function G+. The second integral is a contribution specific to the massive 
case, which does not vanish outside of the light-cone of the source since it involves the unlocalized 
transverse current jj. Therefore, this second integral extends over the whole space and represents an 
instantaneous contribution to the electric field. 

It is worth noting that the presence of this physical instantaneous interaction is related to the 
breaking of Lorentz-invariance. Without breaking of this symmetry, the modifications to the dis- 
persion relation of Aq and would be proportional and the instantaneous contributions of these 
potentials would cancel each other in the expression (4.1). The discussion of physical consequences 
of instantaneous interactions in massive electrodynamics can be found in [157]. 

4.2 Instantaneous gravitational potentials 

The emergence of a physical instantaneous interaction in the massive gravity model described by 
(2.29) is also a consequence of the non-cancellation of unlocalized contributions [160], whereas those 
contributions cancel in GR. Therefore, the same approach as the one used in the previous section will 
be useful to understand how a physical instantaneous interaction appears in this model. For simplicity 
and because of the particular interest of those models, we will focus in the current discussion on 
massive gravity theories based on a function of Z*-' (action 3.16). 

As already mentioned, a small perturbation ST^^y of the matter energy-momentum tensor will pro- 
duce metric and Goldstone perturbations 5g^y and Scf)^ of the vacuum solution, just as a non- vanishing 
current is responsible for the electromagnetic potentials It follows from the discussion of sec- 
tion 3.2 that the tensor potential Hij is instantaneous for the same reason that the potential Aj of 
electrodynamics. Indeed, in all physical situations 5Tij is localized in space. But its transverse and 
traceless part TTjj extends over the whole space, implying a non-vanishing tensor potential Hij ev- 
erywhere in space. Relation (3.5) makes it clear that the vector potential voi is also an instantaneous 
potential. By making use of the the Green functions of the d'Alembert and Klein-Gordon equations 




(4.1) 
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equation, the tensor and vector potentials (3.4) and (3.6) are given by 
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dV 
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G± (t,x,t',x') [iTi{t',x')-Ci{t',x')] 



(4.2) 
(4.3) 



It is also obvious from eqs. (3. 17) that the two gravitational potentials and $ are two instantaneous 
potentials, such as the Coulomb potential Aq—A^ of electrodynamics. With the account of the energy- 
momentum conservation, and by making use of the Green's function of the d' Alembert equation, these 
scalar potentials read 
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5p {t',x') + 35p{t',x') 



(4.4) 
(4.5) 



^0 has been neglected since this arbitrary function does not play any role in the interaction between 
gravitational sources. 

The modification of the dispersion relation of Hij induced by m2 is analogous to the modifications 
of the dispersion relation of the vector field Aj of electrodynamics induced by mi. Since there is no 
modifications of dispersion relations as compared to GR in the other sectors of the theory due to the 
breaking of the Lorentz-invariance, the instantaneous contributions to the potentials will not cancel 
themselves in the graviton propagator. Therefore, it is necessary to search for an observable involving 
Hij in order to see a physical instantaneous interaction. 



4.3 Frequency shift 

To illustrate the physical instantaneous interaction present in our class of massive gravity models, 
consider the measured energy of a light beam. Let be the four-velocity of an observer, with f ''f ^ = 
1, and the vector tangent to the light's geodesic, with u^^u^ = 0. The light's momentum is defined 
by p'^ = loqu'^, where ujq is a constant, and the frequency measured by the observer is given by 



to 



(4.6) 



In the vacuum T^j^ = and the solution (2.18) holds. Then, the affine parameter Tob of the observer 
can be chosen such that v^^ = (1, 0, 0, 0) while the affine parameters Tph of the light's geodesic can 
be chosen such that = (l, n*) with n? = 1. Therefore, the frequency measured by the observer is 

UJ = Wq. 

If there is a small extra source for gravity T^y = (57^jv, the vacuum solution is slightly modified. 
Because of the geodesic equations, the perturbations of the metric imply perturbations 5v^^ of the four- 
velocity of the observer and perturbations 5u^ of the vector tangent to the light's geodesic. Therefore, 
5v^ and 8u^ are given in terms of the metric perturbations and there is a shift of the measured fre- 
quency, 5uj = Lo {5v^u^ + v^5u^). At the linearized level, eq. (4.6) gives 



5lv 



dr, 



ph 



^ - n'di^ + n'n^ ( I Hij - diW^ 



(4.7) 
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Here, we only present the results; the details of the calculations can be found in appendix B. In a 
FLRW background, this relation is known as the Sachs-Wolfe effect [161]. From (4.7) it is obvious 
that any difference in the tensor sector of the theory as compared to GR will lead to a modification of 
the measured spectral shift. 



4.3.1 Frequency shift in GR 



In GR, the vector and scalar potentials are given by relations (4.3), (4.4) and (4.5). The tensor potential 
Hij is a solution to eq. (3.3) with m| = 0. Therefore, this potential is given by eq. (4.2) where the 
Green function is replaced by the Green function G+, since when m2 0. With the 

account of these four relations, the frequency shift measured by the observer in GR reads 



5u:^^ {t, x) 



A^x' 



G+ (t,x,t',x') 



bTkk (t',x')-<5Too(t',x') 



5Tij {t',x') 



dr, 



ph 



dVG+ {t,x,t',x') 



_d_ 
dx'' 



n 



+26%j{t',x')] . 



{5Tjk {t',x') + 6jk5Too {t',x')) 

(4.8) 



The instantaneous contributions to the potentials cancel each others in this relation. Indeed, this 
expression involves only localized sources, i.e., the components of the localized energy-momentum 
tensor ST^i,. The effects of these sources are retarded because of the retarded Green function G'^. As 
a consequence, in GR the frequency shift cancels outside of the light cone of the localized source of 
gravitational perturbations, and there is no instantaneous interaction. 



4.3.2 Frequency shift in massive gravity 

The situation is different in massive gravity models described by the action (3.16) for which m| / 0. 
Indeed, the vector and scalar potentials are still given by relations (4.3), (4.4) and (4.5). But the tensor 
potential is given by eq. (4.2) with m"^ ^ 0. Therefore, the spectral shift measured by an observer of 
velocity is expressed trough 

— = + / dTph—H,^, (4.9) 

where = Hij — -ffijloR- Hence, there is only one extra term in relation (4.9) as compared to GR. 
This term extends over the whole space and is responsible for an instantaneous shift of frequency. 
Therefore, the shift is observable outside of the light cone of the source. At this point, let us stress that 
this extra term looks like the extra electric field (4.1) obtained in the Lorentz- violating electrodynamics 
model discussed before. Indeed, these two extra terms involve the same function AG^ with sources 
which extend over the whole space 

{t, x) = -2 [ ^AG+ {t, X, t', x') TTi, {t', x') . (4.10) 
J Mpj 

To demonstrate that this term is responsible for an instantaneous shift of frequency, we will focus 
for simplicity reasons on the derivative of the spectral shift with respect to the affine parameters Tph 
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of the Ught geodesic 




5lo 5uj^^ 



(4.11) 



This quantity will now be determined for a concrete example of gravitational source. 
4.3.3 Example of instantaneous interaction 

As an explicit example of physical instantaneous interaction, suppose we have a source which appears 
at i = 0, described by the following energy-momentum tensor 



where Xg = (0, 0, d). The presence of the Heaviside function in these expressions guarantees that 
^T^j^u = before t = 0. Let the observer be located at = and the light's geodesic be along the x 
direction, with n* = —5\. 



Since the distance between the observer and the source is d, and because the speed of light is one in 
our units c = 1, in GR the observer will not measure any spectral shift before t = d. 

In the model (3.16), there is an instantaneous interaction and (4. 11) is not zero for t < d. In order 
to determine the frequency shift, the integral (4.10) has to be evaluated for the source (4.12). If the 
source is chosen such that = 82162 j, the only relevant component of the transverse and traceless 
tensor vTjj is given by^ 



Hence, using relations (4.9) and (4.10), the derivative of the spectral shift with respect to the affine 
parameters of the light geodesic is expressed through 



^For an arbitrary energy-momentum tensor ST/j,^, the transverse and traceless anisotropic stress tensor tt^ can be deter- 
mined using transverse and longitudinal operators. See appendix B.3. 




2fiijt^e (t) didj6^ {x - xs) 
Anijte (t) dj6'^ {x - x,) , 



(4.12) 
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Since AG^ involves the Bessel function J'l {m2x) divided by its argument and since m2X <ti 1, the 
following Taylor expansion can be used to find an approximation of this last expression 

J, {m,x) ^ 1 _ m|r^ 3^3 
m2X 2 16 V 2 ; 

With the account of all these relations, the measured variation of the frequency, caused by the source 
located at distance d of the observer, is simply given for t < d (outside of the light cone of the source) 
by 

= -«74i* + ^H), ^<t<d. 
647rMp[ 

This last expression clearly shows a measurable instantaneous interaction, proportional to the mass of 
the graviton to the fourth power. 



4.4 Origin of the instantaneous interaction 

The origin of this physical instantaneous interaction can be traced back to the three scalar fields 0* 
which couple to the metric through W^'-' in F (X, W^^) . For the sake of argument, let us redefine W^^ 
as 

W'^ = i^g^^'' - e ^ ) d^^'d,<l>>, (4.13) 

where e < 1 is a constant. The standard definition of W^^ correspond to e = 1. Only for this value 
of e are massive gravity models invariant under the symmetry (2.28) which guarantees the absence of 
instabilities. It follows from (4.13) that the three fields 0* live, and therefore propagate, in the effective 
metric 

JC 

Consequently, in the vacuum (2. 1 8) these three fields propagate at a velocity given by 

u! = k 



1-e ./l^e 

Hence, for e < these fields have a propagation velocity smaller than the speed of light while for 
< e < 1 they propagate faster than light. For the models of interest, e = 1 and these fields have an 
infinite velocity. 

The presence of fields with an infinite propagation velocity is not enough to account for an instan- 
taneous interaction. For instance, there is no modifications of the scalar and vector sectors of massive 
gravity models possessing the dilatation symmetry (3.15), with consequence that only the tensor part 
of the metric couples to these instantaneous fields. Hence, it is not surprising that one has to look 
for an observable involving tensor modes in order to see an instantaneous interaction. Moreover, if 
one gives up the requirement of the invariance under the dilatation symmetry, Newton's potential is 
modified because of a direct coupling between the scalar part of the metric and these instantaneous 
fields. Hence, an instantaneous interaction is also expected in the scalar sector for those models. 

Before going to the conclusions of this chapter, let us consider the physical instantaneous interac- 
tion from a "fluid mechanics"point of view for arbitrary function (X, W"^). The generic notations 
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for the energy-momentum tensor introduced in section 3.2 can also be used for the energy-momentum 
tensor of the Goldstone fields 0'^, 



t 



111/ 



+ 



I v,,v. 



li'Jv 



The energy, pressure and anisotropic stress densities for the Goldstone fluid are given by 
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2 ~ ^ dW^i ) 



(4.15) 



respectively. We know from OR that the tensor modes only couple to matter fields through the 
anisotropic stress tensor. Hence, these relations make clear that it is necessary to have a function 
T depending on W^^ in order for the scalar fields to give a mass to gravitons. 

The perturbations of the Goldstone fluid belong to two categories: isotropic sound waves which 
are perturbations of the energy and pressure densities and anisotropic sound waves originated in small 
modifications of vrj^i,. In the vacuum, the energy and pressure perturbations 6p^ = dt^^v^v^ and 
= ^t/iuj'^'^ are given by 



M 



pi 



M 



pi 



mi { <^ -Eq] + (m: 



3mq 



1 



After substituting in these relations the solutions (3.11), (3.12), (3.13) and (3.14) of the EoM, one 
finds that the energy and pressure density perturbations of the Goldstone fluid are 



0, 



3p — 2 



mo 



3/. "''df 



Therefore, a perturbation 5p of the matter energy density produces an isotropic sound wave in the 
Goldstone fluid which propagates at an infinite velocity 



oo. 



Hence, there is an instantaneous interaction in the scalar sector. Note that massive gravity models 
possessing the dilatation symmetry (3.15) are an exception. For these models, 6p(j) = with conse- 
quences that there is no isotropic sound wave at all. 

If the matter energy-momentum tensor has a non-trivial anisotropic stress (^vr^j^ = STx^j^f^p'^ — 
j^^ 5T\^j^^ /?,, it will generate an anisotropic sound wave in the Goldstone fluid (i^ttq^ = 0) 

Ml.ml 



djUi + diOj + Hij + - [Widj - Sijdl) E 



which for the class of models considered here becomes 



ml I dVG+ {t,x,t',x') vr,, {t',x') + '^-^mf-^^%J^6p. 



6p 



53 



There are two terms in this relation, both responsible for an instantaneous anisotropic sound wave. 
While the second term vanishes for massive gravity models possessing the dilatation symmetry (3. 15), 
the first term is the origin of the physical instantaneous interaction observed in the previous sections. 
Indeed, this anisotropic wave propagates at an infinite velocity because of the unlocalized transverse 
and traceless tensor vTjj which implies that Sirfj has a non-vanishing value in the entire space. 

4.5 Summary and prospects 

Let us conclude this chapter by summarizing the previous discussion. It has then been argued that 
physical instantaneous interactions are present in massive gravity models described by the function 
f (Z*-^). The origin of this interaction lies in an anisotropic sound wave which propagates at an in- 
finite velocity in the Goldstone fluid. The existence of such interactions is related to the spontaneous 
breaking of the Lorentz-symmetry induced by space-time dependent vacuum expectation values of 
the four Goldstone scalar fields (p'^. Indeed, this symmetry breaking allows for modifications of the 
dispersion relation of the tensor modes as compared to OR without affecting the dispersion relations 
of other sectors. As a consequence, the instantaneous contributions to the potentials do not can- 
cel in the graviton propagator, unlike in GR. This is analogous to the situation in Lorentz-violating 
electrodynamics, where the instantaneous contributions to the potentials do not cancel in the photon 
propagator. 

It has been demonstrated here that a gravitational source localized in space is responsible for an 
instantaneous frequency shift of light beams seen by an observer. It is worth noting that, at the lin- 
earized level, the amplitude of this instantaneous spectral shift is proportional to the graviton mass 
to the fourth power, and therefore is very small as compared to the usual retarded interaction be- 
tween gravitational sources. Regardless of the strength, there is no causal paradox associated with 
superluminal propagation in this model as the breaking of the Lorentz symmetry implies a preferred 
frame. 

The action (3.16) is a low-energy effective action. One should then expect corrections contain- 
ing higher-derivative terms to be present in massive gravity models. These corrections could not be 
responsible for modifications of the dispersion relations of the scalar and vector fields proportional 
to m^. Therefore, the instantaneous contributions coming from the mass of the graviton cannot be 
canceled by higher-derivative terms, and the inclusion of such terms will not affect the conclusion of 
this work. 

Finally, the presence of the physical instantaneous interaction is supposed to be responsible for 
the violation of the black hole "no-hair"theorem present in these models [98] (the physics of black 
holes will be discussed in chapter 5). Indeed, a crucial feature of black holes in GR is the absence of 
"hair", which simply means that a black hole is entirely characterized by its mass, charge and angular 
momentum. This "no-hair"theorem is a consequence of the existence of an horizon which divides the 
space-time in two regions: the outside and the inside of the black hole. The causal structure of the 
metric implies then that no signal with velocity smaller or equal to the velocity of light can escape from 
the inside of the black hole. It seems reasonable to expect that an instantaneous interaction such as the 
one discussed here could carry information through the black hole horizon while interactions which 
propagate at finite velocities v < c cannot. Hence, one does not expect the "no-hair"theorem to hold 
in Lorentz-violating massive gravity. It is worth noting that this instantaneous interaction should also 
allow to look behind the cosmological horizon, since the cosmological horizon is somehow similar to 
the black hole horizon. It would be interesting to study this last issue in more detail. 
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5 

"BiackjXofzs 



In GR, a crucial role is played by the spherically symmetric vacuum solution to the Einstein equations 
- the Schwarzschild solution. This role is twofold. First, this solution describes the metric outside 
of spherical non-rotating bodies and gives rise, in the weak field limit, to Newtonian gravity. It pro- 
vides therefore a useful approximation in many astrophysical situations. Second, the Schwarzschild 
solution describes the result of the gravitational collapse of massive stars, black holes [162, 163, 164]. 
Although the existence of black holes has not yet been directly confirmed, there exists indirect evi- 
dences that some of binary stellar systems contain black holes as one of the companion [165, 166], 
and that many galaxies, including the Milky Way, harbor super-massive black holes in their centers 
[167, 168]. It is conceivable that black holes will be directly observed in the near future, and that their 
properties, including the metric configuration near the horizon, will be quantitatively tested [ 1 69, 1 70] , 
thus providing a probe of GR in a fully non-linear regime. 

5.1 Hairy black holes 

Because of the horizon enclosing each black holes that acts as a "one-way membrane", the properties 
of black holes are extremely resistant to any modifications. Particles and radiation can enter the 
horizon from outside, but escaping to the exterior is impossible. Consequently, information cannot 
get away from a black hole and most of the physical properties of an object (or field) falling into 
a black hole become unobservable once the object (or field) crosses the horizon. For instance, the 
baryon number of a molecule falling into a black hole is not transferred to the black hole once the 
molecule enters the horizon: it simply becomes unobservable [171, 172]. 

The uniqueness theorems for the Reissner-Nordstrom [173] and Kerr [174, 175] black holes fol- 
low from this particularity of the horizon, as the statement by Wheeler that ''Black holes have no 
hair"{\16\. These theorems claim that any collapsing object will result in a stationary black hole 
completely characterized by its mass, electric charge and angular momentum, regardless of the details 
of the collapse or the properties of the collapsing object [177, 178, 179] It is natural that these 
quantities remain observable since they are all conserved quantities subject to a Gauss law, implying 
that they can be measured by a distant observer. Hence, if one defines hair as those free parameters of 
the black hole which are not subject to a Gauss law [181], one indeed concludes that black holes have 
no hair in GR. 

The Schwarzschild metric, together with properly arranged scalar fields, is a solution to Einstein 

'For a recent discussion of stationary black holes with global charges, see [1 80] and references therein. 
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equations in massive gravity as well [98]. The following solution for the scalar fields 



It 

= t + 2Jrr~s + rs\n 



x\ 



where is the Schwarzschild radius of the black hole, together with the Schwarzschild metric im- 
ply that the energy-momentum tensor for the Goldstone fields vanishes outside of the black hole-. 
Hence, this field configuration is a solution to the massive gravitational field equations. Still, the 
properties of black holes are expected to be different in massive gravity with spontaneous break- 
ing of Lorentz-invariance. For instance, rotating black holes carry a long-range tensor component 
(quadrupole moment) which is certainly modified when tensor modes acquire a mass [98]. 

More generally, black holes are expected to have hair [98] in massive gravity theories because 
of the physical instantaneous interaction present in these models. As already discussed, this instanta- 
neous interaction is supported by the three fields 0* which propagate in the effective metric (4. 14) with 
e = 1. Hence, these fields have an infinite propagation velocity in the vacuum. In order to understand 
the dynamics of these fields in the Schwarzschild space-time, let us write the Schwarzschild solution 
of massive gravity in a coordinate system regular at the horizon : 

ds^ = dr^ -di?2 _ r-2 (r, i?, r,) dO^ 

r (r, R, Ts) 



r (r, R, Ts) X (cos if sin 0, sin sin 9, cos 9) , 



where 



r (r, R, rs 



2/3 



This coordinate system, known as the Lemaitre reference frame, corresponds to the frame of free 
falling observers [164]. In this frame, the effective metric G^'^ in which the fields 0* propagate is 
simply given by 

dsl = (1 - e) dr2 - . . dR^ - r2 (r, R, r,) 

r (r, R, rs) 

Although we are interested in the case e = 1, it is instructive to keep e general for the sake of argument. 
The coefficient in front of dr^ can be absorbed by the following rescaling 

^ R ^ 



with consequence that the effective metric reduces to the Lemaitre one with a different Schwarzschild 
radius fs 

ds% = dr^ - _ dfi^ - (r, i?, f,) dO^ = (1 - e) r,. 

r{T,R,rs) 



"For this Schwarzschild solution of massive gravity, it is straightforward to show that X = 1 while Tr {W") — (— )"3 as 
in Minkowski space-time. Hence, the energy-momentum tensor for the Goldstone fields reduces to its value in the vacuum 
which is zero. 
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This effective metric has an horizon at r {t, R,rs) = r^. Consequently, black holes appear larger for 
sub-luminal fields (e < 0) and smaller for superluminal (0 < e < 1) ones. Although the previous 
coordinates redefinition becomes singular in the limit e ^ 1, one may wonder what is the value of 
rs in this limit which corresponds to the instantaneous fields of interest: the effective Schwarzschild 
radius fs goes to zero when e ^ 1. Hence, one may conclude that there is no horizon for these fields, 
and therefore that they can escape from black holes. 

As expected, this short discussion implies that instantaneous interaction can carry information 
outside of black holes, which should therefore have a large number of hair. This conclusion, based 
perhaps on naive arguments, is reinforced by direct calculations made in the context of GR with the 
model of Lorentz- violating electrodynamics introduced in section 4. 1 for which there is an instanta- 
neous electromagnetic interaction. Indeed, by considering neutral non-spherically symmetric black 
holes, the authors of Ref. [98] have demonstrated the existence of hair which can be interpreted as the 
electric dipole moment of the black holes. 

There are also thermodynamical arguments in favor of black hole hair in massive gravity mod- 
els with spontaneous breaking of Lorentz symmetry. The fact that the black hole horizon appears 
larger for sub-luminal particles and smaller for superluminal ones implies that the temperature of the 
Hawking radiation emitted by the black hole depends on the propagation velocity of the fields emitted 
[151]. Hence, the temperature of this radiation is not universal any longer, and one may think of 
processes such that the black hole mass and angular momentum remain constant while the entropy 
outside decreases [151, 182], invalidating the second law of thermodynamics. Still, on top of reduc- 
ing the entropy outside of the black hole, such processes also produce change inside the black hole. 
Consequently, the validity of the second law could be restore if these changes are observable from 
outside, i.e., if black holes have hair on top of their mass and angular momentum [98]. 

The possible existence of hairy black hole in massive gravity models suggests that there might 
exist spherically symmetric solutions other than the Schwarzschild one. This supposition has been 
confirmed by finding of exact spherically symmetric vacuum solutions in massive gravity theories 
[183] corresponding to modified black holes. This chapter is devoted to these solutions. 

5.2 Static spherically symmetric ansatz and equations 

The first step in the determination of exact static spherically symmetric solutions consist in finding 
a good ansatz to describe them. The ansatz for the metric is the same as in GR. Hence, the only 
subtlety consist in getting the correct configuration for the four scalar fields. With the knowledge of 
the vacuum solution (2.18), it is possible to convince oneself that any static spherically symmetric 
configuration in massive gravity models can be written in the following form, 

ds^ = a (r ) dt^ + 2(5 (r) dtdr - /3 (r ) dr^ - k (r) dQ^ , 
/ = t + h{r), 

r 

This field configuration is invariant under two residual coordinate transformations. The first one is an 
arbitrary change of the radial coordinate r r' = r' (r) which allows to set either k, = r"^ or 4> = r. 
The second one consist in redefining the time variable t — > t' = t + r (r). This last transformation 
allows one to cancel either 5 (r) or h (r). We choose the conditions k = r"^ and 5 = 0. Thus, we get 
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the following ansatz, 

/ = t + h{r), 

= ^{r)-. (5.1) 
r 

As compared to GR, this configuration contains two additional radial functions, h (r) and cf) (r). 

As has been pointed out in section 2.2.2, the rotational invariance of the vacuum (and likewise, of 
the ansatz (5.1)) requires that the function (X, W^^) depends on W^^ through three scalar combi- 
nations Wn = Tr {W^). Hence, the function JT depends on four scalar variables which are expressed 
in terms of the radial functions and their derivatives as follows. 



where the two functions /i and /2 are given by 

a/3X'' 



/l — ^ ov ^ h — -2 



In these expressions and in what follows, the prime denotes the derivative with respect to the radial 
coordinate r. 

With the ansatz (5.1), six of the ten Einstein equations are identically satisfied. Consequently, 
there is only four equations left to be solved for four unknown radial functions, 

CO _ J_^0 nr _ J_.r 



Qti _ 1 .e .r _ n 



pi 

The explicit expressions for the components of and are given in the appendix A.4. Consider 
first the equation = 0. Assuming h! ^ 0, this equation gives 

= XJ^x + fi (^1 - 2/1.F2 + 3/^.^3) , (5.2) 

where J^x = dJ^/dX as previously while = dJ^/dwi. Furthermore, the time and radial com- 
ponents of the energy-momentum tensor differ by the quantity proportional to eq. (5.2). Therefore, 
when this equation holds one has = t^. This implies that = or, equivalently, 

a (r) (j (r) = 1, 

in full analogy with the Schwarzschild solution in GR. Hence, the ten Einstein equations reduce to the 
following four equations, 

1 = a(3, (5.3) 

= - + "^-'^{J^-2XTx), (5.4) 
/ // 2 

= - + ^-'^{J^ + XJ^x-WiJ^i-2w2J'2-^W3J^3), (5.5) 

r 2 2 

= + /i (^1 - 2/1.F2 + 3/^.^3) , (5.6) 
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For a generic function T, this system of equations is well defined. Indeed, since the function h{r) 
enters eqs. (5.3 - 5.6) only through the variable X, one may consider X as an independent variable 
instead of h{r). Then the fourth equation allows to find in terms of X, while the first equation gives 
/9 in terms of a. The second equation then gives X in terms of a and the third equation allows to 
determine a as a function of r. 



5.3 Analytical example 

Finding analytical solutions of the non-linear system of equations like (5.3-5.6) is impossible for a 
generic function T . So, in order to get some insight into the behavior of the solutions, let us choose 
the function T in such a way that the resulting equations are solvable analytically. 
Consider the function T of the following form, 

^ = Co + + ci {w\ - 3wiW2 - 6wi + - 12) , (5.7) 

where co is an arbitrary dimensionless constant and ci = ±1 (the numerical value of ci can be 
absorbed into the constant A). The coefficients inside the parentheses are chosen in such a way that 
the vacuum (2.18) is a solution to the Einstein equations. Our example contains, therefore, a single 
continuous free parameter co. Two additional constraints should be imposed on cq. The first one 
comes from the requirement that the graviton is non-tachyonic, rn^ > 0. This translates into the 
inequality 

Co - 6ci > 0. (5.8) 

The second condition is necessary to ensure that scalar modes with pathological behavior do not 
reappear upon addition of higher-derivative terms (recall that the model (2.29) is understood as the 
low-energy effective theory, so such terms are generically present). This condition is given by the 
inequality (2.30) which in the case at hand implies 

CO > 0. (5.9) 

Before solving the Einstein equations for the particular model of massive gravity based on the 
function (5.7), we should stress that the choice of this functional form is by no means unique. This 
particular form has been chosen in order to simplify the solution of the field equations, as will become 
clear below. 

5.3.1 The static spherically symmetric solutions 

Let us start with eq. (5.6) which has motivated our choice of the particular functional form (5.7). This 
equation reads 



1 



-Co + (t)"^ ( 6ci ^ Co - 6c] 



Because of our choice of the function T, this equation contains the variable X as an overall factor 
only. Hence, it reduces to a closed differential equation for cj). The solution to this equation is 

= 6r, 
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where the constant b satisfies the equation 



= (62 _ 1) 66^ + Qb^ + 



Co 

Cl 



(5.10) 



We are interested in real positive values^^ of b. For cq/ci > there is only one such solution, 

6 = 1, 



while for co/ci < there exists another one 

& = ^ (-1 + - 2co/3ci j . 

Thus, at negative cq/ci we have two different branches of solutions. Then, the remaining two equa- 
tions (5.4) and (5.5) can be written as follows 








a a 
— + — 
r 1 

a" + X- 



1 



3Ac + m^co 1 



+ 



a 



X 



3Ac (2 + A) , 



(5.11) 
(5.12) 



where 



A 



-126' 



Co' 



A, 



2m^ci (b^ 



1) 



Recall that, according to our normalization, the constant ci only takes two values, ci = ±1. It is 
worth noting that Ac corresponds to the cosmological constant present in the model. 

Eq. (5 . 1 2) is a linear inhomogeneous equation for a. Its general solution can be found analytically. 
Making use of this solution and integrating the remaining equations one finally obtains 



adt 

t± 
bx\ 



dr 

a 



r^dO^ 
S A 



a 



+ 1 



1/2 



(5.13) 



where 



a (r) 



1 



rs_ 
r 



S A 2 



(5.14) 



Here and S are two integration constants: is the usual Schwarzschild radius while 5" is a scalar 
charge whose presence reflects the modification of the gravitational interaction as compared to GR. 
At S" = this solution reduces to the conventional Schwarzschild solution describing a black hole of 
mass M = Ts (2G)^^. The behavior of this class of solutions is determined by the two integration 
constants and S, and the value of the parameter A. Let us discussed their properties. 



The case 6 < can be reduced to 6 > by the inversion of coordinates. 
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Asymptotic behavior 



The behavior of these solutions at spatial infinity depends on the constants cq and ci = ±1. Models 
with ci = 1 require cq > 6 to ensure that the graviton is non-tachyonic. Then, the only solution of 
eq. (5.10) is 6 = 1. Hence, there is no cosmological constant Ac = and A < so that the metric is 
growing at spatial infinity as Sr^'^^. Such solutions do not describe asymptotically flat space-time and 
will therefore not be discussed further. 

Models characterized by ci = — 1 are more interesting, provided that cq > in order to satisfy 
eqs. (5.8) and (5.9). In that case A > and two branches of solutions exist. The first one posses a 
non-zero Ac, which can be positive or negative depending on the numerical value of cq. These solu- 
tions asymptote to the generalization of the vacuum solution in presence of a non-zero cosmological 
constant. The second branch of solutions are characterized by Ac = (and therefore by 6 = 1). They 
asymptote to Minkowski space-time. We will only consider solutions with Ac = in what follows. 

The black-hole mass 

The total mass of the static, asymptotically flat space-time under consideration is given by the Komar 
mass [9] which takes the following form 



• If < A < 1, the third term on the right-hand side of eq. (5.14) dominates at large distances 
with consequence that the Komar mass of these solutions grows with distances. 

• If A > 1, the standard Schwarzschild term dominates at infinity implying that the Komar mass 
tends to M = r s (2G)^^ when r ^ oo. The solutions with positive (negative) M have attrac- 
tive (repulsive) behavior at infinity. 

Behavior of the solutions 

At the origin r = both terms proportional to and S are singular, so the metric always possesses a 
singularity unless Vg = S = 0. This singularity may or may not be hidden by the horizon depending 
on the signs and values of and S. The solutions possessing the horizon are candidates for modified 
black holes. 

The horizon is always present if both r<j and S are positive. Such black holes have attractive 
gravitational potential at all distances, which is stronger than for a conventional black hole of mass 
M. The horizon size of the modified black hole is larger than = 2GM. Solution where both 
and S are negative imply to naked singularities and will therefore be disregarded. Let us study models 
with A > 1 and < A < 1 separately. 

• For models characterized by A > 1, the standard Schwarzschild term dominates at infinity. If 

> and < 0, the presence of the horizon depends on the relative values of S and r^. 
It exists for sufficiently small l^l. Defining, on dimensional grounds, the mass parameter s 
associated with S by the relation \S\ = s^'^, the existence of an horizon requires that 




A- 1 




(5.15) 
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Figure 5.1: The deviation of goo from one (proportional to the Newtonian potential goo — 1 = 2$) for 
tree different situations. For models characterized by A > 1, Figs, (a) and (b) correspond to > 
and S < with the numerical values satisfying (a) and not satisfying (b) eq. (5.15), while Fig. (c) 
represents solution with < and S* > 0. For models characterized by < A < 1, Figs, (a) and 
(b) correspond to < and S > with the numerical values satisfying (a) and not satisfying (b) 
eq. (5.16), while Fig. (c) represents solution with Vg > and S < 



The Newtonian potentials for solutions satisfying and not satisfying the condition (5.15) are 
shown in Fig. 5.1 - (a) and (b), respectively. When the horizon exists, the gravitational field 
is attractive all the way to the horizon. The attraction is weaker than in the case of the usual 
Schwarzschild black hole of mass r^, and the horizon size is smaller. The behavior of the 
gravitational force with distance mimics that of the smaller-mass black hole plus a continuous 
distribution of "dark matter", with the total mass enclosed within the radius r approaching M 
as r ^ oo. 

At < and S > 0, the modified black hole anti-gravitates at large distances and gravitates 
close to the horizon. The attraction changes to repulsion at 

1 

A5 



A-1 



The corresponding Newtonian potential is shown in Fig. 5.1 - (c). 

• For models characterized by < A < 1, the new contribution proportional to the scalar charge 
S dominates at infinity. If > and S < 0, the modified black hole anti-gravitates at large 
distances and gravitates close to the horizon. As before, the attraction changes to repulsion at 
r = and the corresponding Newtonian potential is shown in Fig. 5.1- (c). 

Finally, at < and 5 > 0, the presence of the horizon depends on the relative values of S 
and Ts- The existence of an horizon requires that 

A-1 



The Newtonian potentials for solutions satisfying and not satisfying the condition (5.16) are 
shown in Fig. 5.1 - (a) and (b), respectively. As previously, when the horizon exists, the gravi- 
tational field is attractive all the way to the horizon. 
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A remark is in order at this point. In the conventional GR the constant or, equivalently, the black 
hole mass M, is also a free parameter which can, in principle, be positive or negative. In GR, however, 
only positive values make sense for the following reasons. First, negative-mass Schwarzschild solu- 
tions possess naked singularity at the origin, which is physically unacceptable. Second, the conven- 
tional matter satisfies the null energy condition which ensures that any compact spherically- symmetric 
matter distribution has a positive mass [184]. None of these arguments go through in the case of mas- 
sive gravity. Fig. 5.1 -(c) gives an example of solution with repulsive behavior at large distances and 
without naked singularity: as for a conventional black hole, the singularity of this solution is hidden 
behind the horizon. The positivity of energy is also not expected in massive gravity. This is related 
to the fact that the background (2.18) breaks time translations, and only the combination of the time 
translations with the shifts of by a constant remains unbroken. In this respect the massive gravity 
model is exactly analogous to the Ghost Condensate model, where the negative-energy states have 
been constructed explicitely [143]. 

5.3.2 Correspondence with linear analysis 

The solutions found in the previous section have the asymptotic behavior different from that obtained 
in the linear perturbation theory discussed in chapter 3. In order to compare the exact solution with the 
results of the linear perturbation theory, let us discuss the exact solution in the gauge where h (r) = 0, 
6 {r) (cf. section 5.2). In this gauge the perturbation theory corresponds to assuming that the 
variations of all metric components are of the same order. In other words, they are formally assigned a 
small parameter e to the first power. The solutions described in chapter 3 satisfy the Einstein equations 
expanded to the linear order in e. 

The solution (5.13) is not of this type. Transforming it into the gauge /i (r) = one finds that 
gor = 6 (r) does not decay as fast as the perturbations of other components, for instance, as q — 1. In 
fact, in the equations expanded in powers of perturbations the terms of order 5'^ balance those linear in 
a — 1. In other words, in the formal expansion of the solution (5.13) in powers of the small parameter 
e the perturbation 6 should be assigned the order ^/e rather than e. Hence, the solution (5.13) is non- 
linear even at large distances from the center. A similar phenomenon has been observed in the context 
of bigravity models in Ref. [185]. 

Another difference between the solution (5.13) and the solution to the linearized equations is that 
the former is static, while in the latter only metric components are static (in the gauge goi = 0). 
Indeed, the scalar fields Eq have time dependence which may be viewed as an accretion of a fluid with 
zero energy-momentum tensor (see eq. (3.14)). 

5.3.3 Gravitational field of a star 

In GR, one may relate the mass of a star to an integration constant of the vacuum solution in the 
exterior space by matching the interior and exterior solutions at the star surface (see, e.g., Ref. [8]). In 
massive gravity, one may try to use the same approach to determine the scalar charge S of an ordinary 
star. The analytical solution in the interior region is required for the matching procedure. 

The star is described, to a good approximation, by a diagonal energy-momentum tensor = 
{p, —p, —p, —p), where p and p are the energy density and pressure inside the star, respectively. This 
energy-momentum tensor is assumed to be responsible for the external gravitational field described by 
eqs. (5.13). Since there is no direct coupling between the ordinary matter and the Goldstone fields, 7^^ 
must be conserved separately, V^TJ' = 0. For simplicity, we take the energy density to be constant 
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a.tr<R, where r = i? is the surface of the star, and zero outside. The pressure p cannot be chosen 
independently; it is determined by the conservation of T^. 

Because of the spherical symmetry, the ansatz (5.1) holds. The Einstein equations in the interior 
of the star are obtained from eqs. (5.3 - 5.5) by adding the contributions of the energy-momentum of 
the star, while eq. (5.6) remains unchanged. The resulting set of equations can be solved analytically. 
The solution reads 



a (r) 
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For simplicity, we have expanded the first equation in powers of p, while the other relations are exact. 
Since the geometry inside the star is regular, the integration constants si and S2 must be set to zero. 

The interior solution has to be matched with the solution (5. 13) at r = R. It is convenient to match 
the variable X which equals 1 in the interior region. In the gauge h{r) = this variable is nothing 
but the g^^ component of the metric. Hence, it must be continuous. Making use of eqs. (5.13) one can 
see that the continuity of X atr = R requires that S = 0. Therefore, the scalar charge of an ordinary 
star is zero. 

It remains an open question how objects (e.g., black holes) with S / can be created. The 
argument given above does not apply to time-dependent configurations, so it is possible that a non- 
zero scalar charge may be acquired during the gravitational collapse. 



5.4 J=-{Z'^) models 

As mentioned earlier, models characterized by the function JT of a single variable = X^W^^ 
are of a particular interest. We discuss in this section the exact static spherically symmetric solutions 
in these models. Our goal is to demonstrate that the solutions found earlier are not specific to the 
particular form of the function (5.7) and exist also in models obeying the dilatation symmetry (3.15). 

In section 5.3 the analytical solutions of eqs. (5.3 - 5.6) were obtained by choosing the function 
in such a way that the dependence on X factors out in eq. (5.6). Since now JT has only one argument, 
the derivatives of T with respect to X and W^^ are no more independent. For this reason we did not 
succeed in constructing non-trivial examples where the Einstein equations are solvable analytically. 
Hence, to demonstrate the existence of unusual solutions we have to use numerical methods. 

Consider the following function T, 

JT = Co I 21 + 2 + ^ I + 2ci ( 4 - 2^2 - 1 - 2-^ ^ I , (5.17) 

were Zn = Tr(Z") are tree independent scalars made out of Z*-', with z„ = (X'^)^ Wn- The co- 
efficients in front of individual terms have been adjusted so that the flat metric and the scalar fields 
given by eq. (2.18) solve the field equations. We are interested in solutions to the field equations that 
asymptote to this vacuum state. 
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Figure 5.2: The deviation of goo from one for three different cases: the usual Schwarzschild solution 
(solid line) and two solutions corresponding to different values of the parameters of the function (5.17) 
(long-dashed and short-dashed lines). The integration constants of these solutions have been chosen 
such that the external horizon is located at r = 1. 

In addition to the adjustments already made, the following inequality should be imposed on the 
coefficients cq and ci to ensure that the graviton is non-tachyonic, 



This guarantees that the square of the graviton mass is non-negative. Moreover, this inequality is 
sufficient for the absence of pathological scalar modes which may appear upon addition of higher- 
derivative terms. As in the previous example, the overall scale of the coefficients cq and ci can be 
absorbed in the parameter A, so without loss of generality we may set ci = ±1. 

For this class of models, the field eqs. (5.3 - 5.6) may be viewed as equations for a{r), (3{r), 
■d{r) = X"^ fi — 1 and g{r) = f2 — 1. Then eq. (5.3) gives /3 in terms of a, while eq. (5.6) enables 
to express g in terms of ??. The two remaining equations form a coupled set of non-linear equations 
for a and they have to be solved numerically. 

The numerical solutions are shown in Fig. 5.2 for different value of the parameters cq and ci. For 
all these graphs, we have assumed that the external horizon is located at r = 1 and that = 100 at the 
horizon. The large value of i? is chosen in order to make the difference between the modified solution 
and the Schwarzschild solution visible on the plot (large values of correspond to large scalar charge 
S of the previous section). 

The plots show the behavior qualitatively similar to that discussed in section 5.3. In particular: 

• if Co = ci = —1, the Newtonian potential 2<I> = a — 1 is attractive at short distances and 
becomes repulsive at larger distances; 

• if Co = 3 and ci = 1, the Newtonian potential is attractive outside of the horizon, and becomes 
repulsive close to the singularity. 

The deviations from the Schwarzschild metric are larger for larger values of the "scalar charge" (pa- 
rameterized by the value of i9 at the horizon). The Schwarzschild solution is recovered at i9 ^ 0. 



Co - 2ci > 0. 



(5.18) 
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5.5 Summary and prospects 



To summarize, there exist spherically symmetric vacuum solutions in massive gravity models which 
depend on two integration constants, the Schwarzschild radius and an extra parameter S which can 
be called the "scalar charge". At zero value of the scalar charge the standard Schwarzschild solution 
is recovered, while at non-zero S the metric is modified with respect to the Schwarzschild case. 

The solutions having non-zero scalar charge exhibit much reacher behavior than the Schwarzschild 
solution in GR. As can be seen from the explicit example of section 5.3, both the short and long 
distance behavior may be modified depending on the parameters of the model. 

Unlike in GR, the solutions may have a negative Komar mass. Such solutions have repulsive grav- 
itational interaction at large distances. At short distances the repulsion may change to attraction and 
give rise to the horizon, hiding the singularity at the origin. Such solutions represent anti-gravitating 
black holes. 

In the case of a positive Komar mass, the S'-dependent contributions may make the gravitational 
attraction weaker at short distances (cf. fig. 5.1 - (a)). In this case the gravitational force decays 
with distance slower than thus mimicking the presence of dark matter. Interestingly, solutions 
with the same value of M but different scalar charge S have different behavior, which corresponds 
to different amount of the apparent "dark matter". This is in contrast with other models possessing 
modifications of the gravitational potential [61, 32, 45], where the modification of the gravitational 
force is determined by the parameters of the model. 

It is currently an open question how objects with non-zero scalar charge may be created. As has 
been argued in Sect. 5.3.3, the absence of direct coupling between the Goldstone fields and ordinary 
matter results in zero scalar charges of static matter distributions. Thus, the gravitational field of 
ordinary stars is described by the S = solutions, i.e., by the standard Schwarzschild metric. This 
may be not the case for black holes, especially the super-massive black holes in the centers of galaxies, 
which may be of primordial origin [186, 1 87] . In any case, this question requires further investigation. 

Another open question is the stability of the modified black hole solutions. Several kinds of 
instabilities may be present. Among perturbations of the solutions there may exist unstable modes 
with the characteristic time scale of order of the horizon size; in this case the interpretation in terms 
of black holes is not possible. Second potential source of problems is generic presence of the higher- 
derivative terms not included in the action (2.29). One has to check that the 2-parameter family of 
modified black holes survives their inclusion. By analogy with the ghost condensate case, one may 
expect that these terms produce at least a slow Jeans-type instability [48], which is not, however, 
dangerous for the black hole interpretation (see discussion of section 2.2.6). Finally, the presence of 
negative mass solutions may lead to instabilities of the quantum-mechanical nature similar to those 
found in Ref. [146]. 

To conclude this list, let us mention also the solutions satisfying h = which were not considered 
here. In this case the ten Einstein equations reduce to three equations for a, j3 and which form a 
(generically) well-defined system. It remains to be seen whether this system has asymptotically flat 
solutions. In any case, the Schwarzschild solution does not belong to this class which is characterized 

by a(3 / 1. 
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6 

Structure formation 



Given that massive gravity theories with spontaneous breaking of Lorentz-invariance pass the most 
obvious constraints, one may wonder if they reproduce correctly more subtle parts of modem cos- 
mology, in particular the formation of structures. In the standard cosmological model based on 
GR, the formation of the different structures observed in our Universe (galaxies, clusters of galax- 
ies, super-clusters, voids, . . . ) is understood through the theory of cosmological perturbations (see 
[102, 152, 188, 189, 190, 191] and references therein) which describes the evolution of the primordial 
inhomogeneities of the cosmic fluids filling the early Universe. 

These primordial inhomogeneities grow with the expansion of the Universe because of the low- 
energy gravitational instability present in GR (the Jeans instability controlled by curvature). This in- 
stability implies that a small over-density will grow by attracting surrounding matter. The gravitational 
force exert by this over-density will then increase with consequence that it will attract more and more 
matter. In an expanding Universe, the growth of the over-density is partially counteracted by the ex- 
pansion which tends to dilute it. Still, with some assumptions about cold dark matter [192, 193, 194], 
the general picture which emerges from this seems to be in agreement with the most recent cosmo- 
logical observations [193, 195]. 

It is not obvious a priori that massive gravity theories could reproduce correctly the mechanism 
of structure formation since the background contains the condensates of four Goldstone fields whose 
perturbations mix with the matter density perturbations. Still, we will show in this chapter that cosmo- 
logical perturbation theory in massive gravity models have predictions similar to those of GR [ 1 60] . 

This chapter is organized as follows. First, we review the cosmological solutions in massive 
gravity with Lorentz-symmetry breaking. The aim of this discussion is to show that models invariant 
under the dilatation symmetry (3.15) are attractors of the cosmological evolution [125]. Then, we 
move to the original part and study perturbations of the cosmological solutions found in such models. 

6.1 Cosmological models 

GR has provided a complete framework for cosmologists to work with. Hence, it is of common 
belief that the laws of gravity rule the Universe's dynamic. Consequently, cosmological solutions 
have played a central role in any theory of gravity since the formulation of GR in the beginning of the 
twenty century. Massive gravity models with spontaneous breaking of the Lorentz symmetry are no 
exceptions. Indeed, cosmological solutions were the first exact solutions of the massive gravitational 
field equations to be studied (after the vacuum solution, of course). 

The Standard Model of cosmology, known as the ACDM model, is based on a few assumptions. 
One of them, the cosmological principle, states that the Universe is homogeneous and isotropic on 
cosmological distances. This assumption fixes the geometry of the Universe to be described by the 
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Friedmann-Lemaitre-Robertson-Walker (FLRW) metric 

I 

ds^ = ai-nf (dv^ - dxidx') with dxidx' = — .dr^ + r^dO^ (6.1) 

1 — Kr^ 

where k = —1,0,1 for an open, flat or closed Universe, respectively. Hence, the cosmological 
principle reduces the dynamics of the Universe to the dynamics of a single variable, the scale factor 
a [t]) . The cosmological principle also imposes constraints on the energy-momentum tensor of any 
cosmological fluid, since only perfect fluids are homogeneous and isotropic. Hence, the energy- 
momentum tensor (3.1) of any cosmological fluid reads 

TfiU = {P + V) VfMViy - Pg^,y, 

where p = T^j^^v^^v'^ and p = T^vj^'^ /'i are the energy and pressure densities measured by an ob- 
server of four- velocity which is comoving with the fluid, respectively. This four- velocity obeys the 
geodesic equation v'^VyV^ = with v^^v^ = 1, implying that the affine parameter of the observer can 
be chosen such that = (a, 0, 0, 0). 

Although the cosmological principle was first based on philosophical considerations, its is now 
supported by several high-precisions cosmological observations, such as Large Sky Survey [42, 196] 
or studies of the CMB anisotropics [43]. Consequently, theories of massive gravity have to contain this 
cosmological principle in their cosmological solutions. For this reason, let us assume the following 
spatially-flat homogeneous and isotropic ansatz [89] for the Goldstone scalar fields, 

/ = 0(r?), 0^ = r(r/)x^ (6.2) 

together with the spatially-flat FLRW metric (A; = in (6.1)). This ansatz cannot be directly gener- 
alized to the case of an open or closed Universe, and it is not known whether massive gravity models 
based on (2.29) admit solutions of these types. With this ansatz, the Einstein equations for the mas- 
sive gravitational field reduce to the usual Friedmann and Raychaudhuri equations of GR with extra 
contributions consisting in an energy and pressure densities for the Goldstone fields: 

2 2 

3W' = ^(/9 + /9</,), 2n' + = -^{p + p^), (6.3) 

where prime denotes the derivative with respect to the conformal time rj. Here Ti = a' /a is the 
conformal Hubble parameter so that H = H/a is the physical Hubble parameter, p and p are the 
total energy and pressure densities of usual matter fields (baryons, photons, . . . ) respectively, while 
/O0 and p^ are the energy and pressure densities for the Goldstone scalar fields (section 4.4). Note that 
the Goldstone anisotropic stress tensor irf,^, is identically zero for this ansatz in agreement with the 
homogeneity and isotropy of the background. 

For an arbitrary function of X and W^^ , the four Goldstone equations obtained by varying the 
action with respect to 0'* reduce to 

= X^/'^d^ [a^X^/^J^x) + 3a^-WJ^w (6.4) 

As before, dT/dW'^ = hjTw while W = -8i^W^^ !?> = a^V^ and X = 0^2^/2 j^^jg jj^^j 
the Goldstone equations are equivalent to the equations of energy-momentum conservation for the 
Goldstone fluid, V ^t^^ = 0. Consequently, those equations can be formulated in the same way as 
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the conservation equations of usual matter fields. For any cosmological fluid including the Goldstone 
fluid, the equations for energy and momentum conservation read 

p[ + 3W {pi + Pi) , where i = baryons, photons, (j)'^,. . ■ . 

Since the conservation of total energy-momentum is a consequence of the Einstein equations, if there 
is N non-interacting usual matter fields in the Universe, there is + 2 independent equations describ- 
ing the Universe's dynamic which may be chosen as 

m^ = :^{p + p^), = p'^ + 3n{p^ + p^), = p'^ + 3n{pn + Pn), (6.5) 

where the indices n = 1, . . . , N labels the usual matter fields and p = J2n Pn- 

Given N equations of state pn = WnPn, one can determined the energy densities of those fields 
as functions of the scale factor /?„ oc a^^^^^^"\ Similarly, the Goldstone equation gives p^ as a 
function of the scale factor. Finally, the Friedmann equation gives the scale factor as a function of 
the conformal time r/. Yet, there are two unknown functions in the Goldstone ansatz, (/){r]) and T(r/), 
with only one Goldstone equation to constrain them. Hence, there is a freedom in choosing one of 
those two functions. This freedom is a direct consequence of the fact that the action is supposed to be 
invariant under the symmetry (2.28). Indeed, the ansatz for is still homogeneous as a shift of the 
spatial coordinates by a constant can be compensate by a (j)^ dependent shift of the fields (p^. 

Note however that this ambiguity can be fixed by specifying the boundary conditions for the fields 
(j)^ at spatial infinity [125]. For instance, imagine that the space is a torus of size L. Then, (p^ would 
have to satisfy some periodicity condition, 

implying r = const. Other boundary conditions may lead to time-depend t^t]). Therefore, the 
ambiguity in choosing different functions r (77) is analogous to the ambiguity in choosing the vacuum 
in theories with flat directions. 



6.1.1 Cosmological attractors 

One choice which seems quite natural, as it looks like the vacuum solution, consist in requiring that 
only varies with time, so that r = 1. Then, in cosmological background the function depends 
only on X = a^'^cf)''^ and W = a~^, and the Goldstone equation (6.4) reduces to 

X^l-T, = (6.6) 

This equation implies that at late time either X or J^x goes to zero. Since the mass parameters mf are 
linear combination of the derivatives of J^, one may wonder whether they go to zero or remain finite 
in the limit a — > cxd. Indeed, mass parameters going to zero will imply that the cutoff scale decreases 
with the scale factor. Similarly, if X goes to infinity, the validity of the low energy effective field 
theory becomes questionable. 

There is a broad class of massive gravity theories for which the mass parameters are finite in the 
limit a — > cxD, so that the effective field theory description remains valid [125]. Let us assume that 
X asymptotes to some power of a at late time. This requirement is not that restrictive since eq. (6.6) 
implies that any algebraic function JT (X, a^^) satisfies it. Then, there is a constant 7 such that the 
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combination a '^X^ goes to a non-zero constant Zq in the limit of growing scale factor a —>■ oo. Let 
us replace X by a new variable Z = a^'^X"' . Then, (X, a~^) becomes a function of Z and a 



If one assumes further that the function JT (Z, a ^) is regular in the limit of growing scale factor, one 
finally find that 

P {Z, a"^) ^ J^o (Zo) when 

This discussion implies that a wide class of functions (X, W^^) depend only on the combination 
Z^^ = X^W^^ in the limit when the scale factor goes to infinity. Hence, there is another motivation 
(beside the correct Newtonian limit) for studying theories which depend on a single argument Z*^ and 
satisfy the dilation symmetry (3.15). 

Let us discuss the cosmological solution of massive gravity possessing this dilatation symmetry. 
For these models, the derivatives of with respect to X and are proportional to one another, 

= ^—^Z, J^W = 77}^Z, 

X W 

where dT/dZ^^ = 6ijTz and Z = —6ijZ'^^/3 = (\)''^'^ j a^^^'^ . The energy and pressure densities for 
the Goldstone fluid read 



with 



ZT7 

^1 + 2 (37 - 1) % . (6.7) 



The Goldstone equation (6.6) reduces to 



a3-i/7- 

For a given function T, this equation determines the dependence of the variable Z^^ on the scale 
factor. One has to specify JT in order to solve it. However, this relation gives some hints about 
the Goldstone dynamics without need of an exact solution. For example, models characterized by 
7 = 1/3 are particular. The Goldstone equation implies then that Z^^ is constant without constraints 
on J^z, while the equation of state (6.7) implies that behaves as a cosmological constant whose 
value is determined by the initial conditions. 

One can also determined some features of models characterized by 7 > 1/3 or 7 < without 
need of an exact solution. Indeed, for those value of 7 the Goldstone equation implies that 

^i-i/27_^^ ^ Q ^j^gj^ a ^ 00. 

Hence, either Z goes to zero so that the mass parameters decrease with the scale factor, or J^z goes 
to zero. The case of interest is when Z goes to a constant Zq in this limit such that J^z (-^o) = 0. 
In this case, the graviton mass remains finite at o — > cxd and the effective field theory description 
remains valid. Moreover, for models satisfying Tz (Zq) = 0, Z = Zq is a. natural solution to the 
Goldstone equation for which 0' oc a^^^/"'. For those models, the Goldstone fluid also corresponds 
to a cosmological constant whose value is fixed by the initial conditions. 

Finally, it is worth noting that for < 7 < 1 /3 and regular function J^, Z grows with the scale 
factor implying that the effective field theory description breaks at late time. 
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6.1.2 Other example 



Other massive gravity theories may imply different cosmological models. For the sake of illustration, 
one may consider solutions with r / 1 before going to the original part of this chapter. Let us have a 
look the following function of X and W"^^ 

T = {w'' - w'm'^) p {X) , 

where P {X) is a function of X only. For this particular model, the mass of the graviton is given 
by ml = 2m?W {1 + P {X). Let us choose the ansatz t = a such that W = 1. Then the 
homogeneous and isotropic Goldstone fluid is pressure-less, = 0, so that the Goldstone equation 
(6.4) reduces to 

P'4, = -3W/)<^ P4, = (6.8) 

where pdm is a constant corresponding to the initial Goldstone energy density. Hence, depending on 
the particular model considered, massive gravity theories could provide a cold dark matter candidate 
in the right hand side of the Friedmann equation 

M^j V J 

Moreover, if P {X) is simply given by 

P{X) = 1 + 

then eq. (6.8) implies that X = 1 — pdm/ (3a^) so that P ^ 2 in the limit a ^ 00. Hence, the graviton 
mass tends to a non-zero constant in the limit of growing scale factor. 



6.2 Cosmological perturbations 

To address the question of structure formation, we focus on massive gravity theories possessing the 
dilatation symmetry (3.15) and neglect possible deviations from the point Z = Zq by considering 
models for which Fz {Zq) = 0. For those models, Z^^ = —ZqS''^ is a solution of the Goldstone equa- 
tions so that the mass parameters mf are constants while = Z^'^''^ a}^^/^ . Then, corresponds to 
a cosmological constant and the EoM (6.5) for the background reduce to 

^2 

^'^^ = ^P + P^) ^ = p; + 3W(p„+p„). (6.9) 
pi 

We also assume that the total cosmological constant p\, which includes p^, is of the order of the 
present-day cosmological constant. Its contribution to the Friedmann equation at the epoch of struc- 
ture formation is therefore negligible. With all these considerations in mind, the solution to eqs. (6.9) 
can be parameterized as 

2/(l+3u;) 

-I , /9 = /5oa-3(i+-), (6.10) 



where po and 770 are constants related through 



- 2 _ 



pi 
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Transverse and longitudinal perturbations 

In order to describe perturbations of a flat FLRW space-time, one has to generalize the decomposition 
(2.3) into transverse and longitudinal fields [102], which then reads 

hoQ = 2o?ip, 
hoi = a2 {Si + diB) , 

hij = {2^Sij - 2didjE - d,Fj - djFi + Hij) , 

where, as in Minkowski space-time, the vector perturbations Si, Fi and are transverse, while the 
tensor perturbation Hij is transverse and traceless. Since FLRW background is time-dependent, the 
gauge-invariant scalar perturbations introduced in section 2.1.1 have to be generalized as follows 

^ = ^-a-^[a{E' + 3)]' , ^ = ^ + H{E' + B), 
= e -^'{B + E'), E = C-E. 

The two gauge-invariant vector fields are still given by the same combinations as in flat space-time, 

= Si + F/, cji = - Fi, 

and Hij is still gauge-invariant as expected from its tensorial structure. 

Perturbations of the energy-momentum tensor 

The notations introduced in section 3.2 in order to describe the energy-momentum perturbations above 
Minkowski space-time need also to be generalized to FLRW background. The main difference as 
compared to the flat space-time originates in the non-zero value of the energy and pressure densities 
of the background, so that the most general perturbations of the energy-momentum tensor of ordinary 
matter fields are parameterized in the following way, 

STij,y = {5p + 5p)Vfj,Vu - g^uSp + (p + p) {vu6v^ + v^6vu) - phfj_„ 

As previously, 5p and 5p are the perturbations of the matter and pressure densities measured by a 
comoving observer, 5q^ is the perturbation of the energy flux perpendicular to and is the per- 
turbation of the anisotropic stress tensor. As in Minkowski space-time, the fact that = (a, 0, 0, 0) 
implies that Sq^ = Sttq^ = 0. The observer's four-velocity is a time-like vector. Then, it satisfies 
Vfj_v^ = 1 which at the linearized level implies that 

8vq = aip. 

Therefore, the linearized energy-momentum tensor of any cosmological fluid is given by 

a^{5p + 2p^), (6.12) 
a{p + p) 6vi - a^p {Si + diB) + a6qi, 

[6ij6p - p {2i;6^j - 2didjE - d,Fj - d^Fi + Hi^)] + 6^^. 

Let us introduce the following parametrization 

a{p + p) 6vi + a6qi = (Ci + diC) , 

5-Kij = [{^didj - Sijdl) it + diiTj + djin + TTij] , 



(6.11) 



5%i 
5Ti 
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where the vector perturbations d and vTj are transverse while the tensor perturbation iTij is trans- 
verse and traceless. At this point, let us comment about the interpretation of Q + diC- In this 
three-dimensional vector field are hidden both the three-dimensional velocity perturbation 6vi and 
the perturbation of the energy flux 6qi. The equations for the gravitational field are not able to differ- 
entiate between these two three-dimensional vector fields, with consequence that only Q + diC will 
be determined by these equations. To differentiate between 6vi and 5qi, one must had an equation for 
6qi. In other words, one has to suppose some physics which will give a non-zero 5qi. The discussion 
of such physics is beyond this work. Therefore, we will assume 6qi = 0. 

In Minkowski space-time, the perturbations of the energy-momentum tensor are all gauge - in- 
variant perturbations. This is a direct consequence of the fact that this tensor is identically zero in 
the vacuum. Since the energy-momentum tensor of usual matter is different from zero in the FLRW 
background, some of the perturbations introduced above are no more gauge-invariant. Indeed, the 
perturbations vTjj, vTj and vr of the anisotropic stress tensor are still gauge invariant as well as the trans- 
verse vector The other perturbations, namely 6p, 6p and ^ are not gauge-invariant. Hence, one 
needs to introduce three gauge-invariant scalars 

6p=[6p-p' {E' + B)]/p, 6p = [6p - p' {E' + B)] /p, 5^ = ( - [p + p) {e' + B) . 

With all these notations, the equations for the conservation of energy and momentum read 

= {pdpY + 3n{p6p+pSp)-3{p + p)^' -df6(;, 
= S[ + m6<: - p6p - 2dl7r - {p + p) , 

= C + '^nCi-dlTTi. (6.13) 
6.2.1 Jeans instability 

Contrary to what have been done in chapter 3, the goal here is not to determine the gravitational 
interaction created by an arbitrary source, but rather to understand how the inhomogeneities described 
by ST^,^ evolve under their own gravitational interaction. 

As a warming up, let us go back to Minkowski space-time and restrict ourselves to Newtonian 
gravity. In the Euclidian space, the usual equations of hydrodynamics [197] for a fluid of density p 
and pressure p are the conservation and Euler equations, which read 

= dtp + di (pvi) , 

= dtVi + VjdjVi + -dip + di^, 

p 

where Vi is the velocity of the fluid, while $ is Newton's potential satisfying the usual Poisson equa- 
tion. A static and uniform fluid is described by the following solution to these two equations: 

p = constant, p = constant, Vi = 0, = 0. 

Any perturbation of the fluid is characterized by an energy 5p, a pressure 6p, a velocity 6vi and a 
potential ^ which are solutions of the following three equations 

= dt5p + pdiSvi, 

= dt5v^ + -di5p + dii^, 
P 

= df^-ATiG5p. 
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Combining these three equations gives a differential equation for 5p known as the sound equation in 
the presence of gravity 



Q={d^,-cldf-A^Gp)5p, cl 



dp/ s 



6p 



where Cg is the speed of sound. The solution of this wave equation is simply given by a plane wave 
whose time dependence is 



6p cx exp. 




1 \t 



where kj = ^ Ai^Gpj is the Jeans wave number. One concludes from this calculation that there 
are two opposite effects in competition here: while gravity tends to increase over-densities, random 
thermal motion described by bp acts as to dilute the over-dense regions. The amplitude of the Jeans 
wave number signals the transition between small perturbations k < kj for which gravity dominates 
and large perturbations k > kj for which gravity is negligible. In the last case, the perturbations are 
just sound waves. Indeed, for k < kj the time dependence of 6p is given by an exponential which 
characterizes the gravitational instability, while for k > kj its time dependence is just an oscillating 
function. In an expending Universe both contributions are diluted by the expansion so that matter 
over-densities grow as a power of the scale factor. 



6.2.2 Evolution of the inhomogeneities 

Theories of structure formation involve two distinct parts. The first one is a theory for the production 
of the primordial inhomogeneities while the second is a model to describe the evolution of these 
primordial inhomogeneities in the expanding Universe. Models of inflation (see [152, 189, 190, 191, 
198, 199] and references therein) are perhaps the most popular candidates as theory for the primordial 
inhomogeneities, while their evolution is understood in the Standard Model of Cosmology through 
GR [152, 189, 190, 191, 193, 195, 200]. 

The present discussion will only concern the evolution of the primordial inhomogeneities. For 
this purpose, one has to consider a system of fourteen independent equations made of ten linearized 
Einstein equations and four linearized Goldstone equations'. In massive gravity models, these equa- 
tions reduce to those of GR in the limit m ^ 0. As in the flat space-time, the equations for the tensor, 
vector and scalar perturbations decouple in the linear regime allowing for a separate study of these 
three sectors. The equations governing each of these sectors are derived in appendix A.3. 



Tensor perturbations 

Although tensor perturbations are not generated by density fluctuations, they are produced in many 
theories of the primordial inhomogeneities in addition to scalar perturbations. For example, gravita- 
tional waves are generated in most models of inflation [201, 202]. If they were produced in the early 
Universe, tensor modes would have left an imprint in the CMB polarization [203, 204, 205]-. Hence, 
it is important to understand how they evolve during the expansion of the Universe. 

'We could also discuss the perturbations by considering the ten linearized Einstein equations along with the four equa- 
tions (6.13) which describe the conservation of the energy-momentum tensor of usual matter Indeed, the Einstein equations 
imply the conservation of the total energy-momentum tensor with consequence that there is effectively fourteen independent 
equations. 

"For a discussion about the influence of the graviton mass on the CMB polarization, see [206, 207]. 
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Figure 6.1: Evolution of gravitational waves in a flat FLRW Universe dominated by usual matter with 
w = as a function of conformal time. On the left, the predictions of GR for three different modes 
labeled by their comoving momenta. On the right, the same modes in massive gravity theories. The 
graviton mass has a negligible influence on small scale while it strongly affects large scale modes. 



The tensor modes are given by an inhomogeneous equation which reads 

= H'lj + 2nH[^ - dfHij + a2 ^m^Fi,- + . (6.14) 

If the mass of the graviton is larger than the Hubble constant, m2 ^ TL/a, which we assume to be the 
case in what follows, eq. (6.14) describes massive gravitational waves. For simplicity, let us neglect 
the anisotropic stress vrjj = 0^^. Then this equation is identical to the equation for a free massive 
scalar field in a flat FLRW background. In the three-dimensional Fourier space, the previous equation 
reduces to 

= H'/j + 2nHlj + k'^Hij, p2»^2^ 

or 

= H^j + 2nH'ij + a^mlHij, < mj, 

in the relativistic and non-relativistic limits, respectively. Note that = a^'^k'^ is the physical mo- 
mentum. These equations state that the amplitude of the gravitational waves scales like Hij oc a^^ 
and Hij oc a~^/^ in the relativistic and non-relativistic limits, respectively. The exact solution of 
eq. (6.14) for a matter dominated Universe are plotted in Fig. 6.1 for both GR and massive gravity 
models. As expected, the graviton mass strongly affects the dynamics of large scale modes (recall 
that X = 27r/k) which are the latest to enter the cosmological horizon, while small scale modes which 
enter the horizon and decay first are practically not affected by the graviton mass. 

Before concluding this discussion about tensor modes, let us mention the possibility of non - 
relativistic gravitational waves to be a dark matter candidate [97], although this possibility seems 
strongly constrained by pulsar observations [209]. 

"^We then ignore, among other effects, the damping provided by cosmological neutrinos [208]. 
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Vector perturbations 

Primordial vector modes describe vortical fluid perturbations in the early Universe. For example, 
primordial magnetic fields may generate vector perturbations through their anisotropic stress [210, 
211] The three gauge-invariant equations describing vector perturbations are 



2 ^Ci 



0, ro' + 2nvji 



2 1 2 , 



0, m2<Ti = 0. 



The first of these equations allows to express Q in terms of wi, while the third equation gives cjj = 0. 
Therefore, the only non-trivial equation is the second one. It differs from the conventional one by 
the term proportional to the graviton mass m\ which cancels at cjj = 0. Thus, this equation is the 
conventional one 



0, 



(6.15) 



and there are no differences in the vector sector as compared to GR. If one neglects the anisotropic 
stress vTj, the previous equation describes a field with the amplitude decreasing as w-i oc a~^. 



Scalar perturbations 

Observations of the CMB anisotropy spectrum shows that primordial scalar perturbations were most 
certainly dominated by adiabatic perturbations [43, 213]. For such perturbations, 6p can be expressed 
in terms of 5p by means of the matter equation of state which reads 



6p = cl 6p 
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where Cs is the sound velocity (c^ = w for a perfect fluid). The scalar sector consists of 8 scalar 
perturbations: H, E^, dp, 6p, (5^ and vr. The behavior of these perturbations is governed by the 

previous equation of state along with four Einstein equations and two Goldstone equations 
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(6.18) 



(6.19) 
(6.20) 
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For a discussion about tiie influence of primordial vector modes on tlie CMB polarization, see [212]. 
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There is therefore 7 equations for 8 scalar perturbations. To close this system, we must either add an 
equation describing the physics behind tt, or we must set vr = 0. We choose the second option here. 
Since we are interested in massive gravity models possessing the dilatation symmetry (3.15), the 



relations mg 



877714 and 771-4 ~ 7 ("^"^a 
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between the mass parameters have been used in the 



previous system of equations. This system can be solved as follows. At 777^ 7^ 0, eq. (6.21) can be 
used to express <I> — r?' j<\)' in terms of the other fields so that eq. (6.20) becomes a closed equation 
for H 



= 5? 



The solution of this equation which does not grow at spatial infinity reads 



7772^ 



(6.22) 



where ^'o (x*) is the function of the space coordinates introduced in section 3.2. As in Minkowski 
space-time, this function arises as an integration constant which is related to in the presence of a mode 
with the dispersion relation uP' = 0. Then, eqs. (6.16), (6.17) and (6.21) can be used to express 5p, (5( 
and H° in terms of ^ and while eq. (6.18) reads 
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With the account of all these relations, the remaining equation (6.19) becomes a closed inhomoge- 
neous equation for 
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This relation is known as the Bardeen equation. Inverting the relation a 
this equations as 



a [t]) enables one to write 
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(6.24) 



Once the solution to this equation is found, the other variables are determined by eqs. (6.16), (6.17), 
(6.21) and (6.23). In particular, if ^' is a solution to eq. (6.24), the density contrast is given by 

t2 



2M 



pi 



- ^aT-L^da - m"^) ^ + a^/^-i 37^ 



,1 

7 
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The conventional cosmological perturbations are recovered by setting the mass parameters to zero, 
= 0. In this case, eq. (6.18) gives <I> — ^' = which implies ^'o = (cf. eqs. (6.22) and 
(6.23)). Then both eqs. (6.24) and (6.25) reduce to the standard equations describing cosmological 
perturbations in the Einstein theory. Note that the value of ^'o is determined essentially by the initial 
conditions. Setting ^'o = would eliminate the ^'o-dependent terms in eqs. (6.24) and (6.25) and 
bring these relations to the conventional form even in the case m\ ^ 0. 

In order to stress the differences between massive gravity models and GR, let us consider different 
situations. 
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Matter perturbations 

For matter perturbations, = and our previous discussion of the Jeans instability suggest that 
no random thermal motion will be able to dilute the over-dense regions since kj = oo. This naive 
expectation is confirmed by the exact solution to eq. (6.24). Indeed, the equation for ^' in a matter- 
dominated Universe {w = 0) reduces to 
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(6.26) 
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which differs from the standard one by the presence of an inhomogeneous term proportional to ^'q. 
The solution to this equation reads 

2 



2 + 37 

where Ci{x'^) are two integration constants. Substituting this solution into eq. (6.25) one finds the 
density contrast 
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The first two terms in this equation are precisely the ones which appear in the standard Einstein the- 
ory, the second term describing the linear growth of the perturbations, 6p oc a. The difference with 
the conventional case consists in the third term on the right hand side of eq. (6.27). The perturba- 
tions corresponding to this term grow proportionally to a^^^. For 7 > 1 or 7 < these "anoma- 
lous"perturbations grow slower than the standard ones. Another case of interest is given by 7 = 1. 
Indeed, in this case the ^'o-term in eq. (6.27) can be absorbed by a redefinition of C2(x*) so that the 
density contrast has the same value as in GR. 



Radiation perturbations 

The situation is similar for a relativistic fluid = 1/3, except for the presence of pressure to compen- 
sate the gravitational instability. This has for consequence that the equations for ^' is more difficult 
to solve analytically. Hence, we will only consider radiation perturbations in the radiation epoch 
w = 1/3. In that case, eq. (6.24) reduces to the following one. 
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In the three-dimensional momentum space, the solution to this equation is given by 
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where n = rjk/\/?, is proportional to the scale factor, Cj(A;*) are two integration constants and /(n, 7) 
is given by 
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The density contrast is obtained by substituting this solution into eq. (6.25) 
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For generic values of 7, this expression is cumbersome. Yet, one may suspect that for 7 < —1 and for 
< 7 the density contrast grows faster than in GR because of the third term in the right hand side of 
this expression. 

For simplicity let us concentrate on modes much smaller than the Hubble scale, fc^ » H^. For 
those modes, n S> 1 and the density contrast reduce to 
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As one may see from this expression, for —1 < 7 < the ^Q-dependent contribution to the density 
contrast decays with the scale factor so that only the standard contribution remains. Thus, in this range 
of 7 the perturbations behave just as predicted by GR in both matter and radiation-dominated epochs. 

The case 7 = 1 is also special here. Indeed, it is straightforward to show that /(n, 1) = ^'o 
plus a term which is absorbed by a redefinition of ci(A;*). Thus, the dependence on "i>o cancels out 
in the density contrast so that only the standard part remains. At other values of 7, the 4'o-dependent 
contributions to the density contrast grow in the radiation - dominated Universe. 



Dark energy perturbations 

Before going to the conclusions of this chapter, let us consider perturbations of a dark energy fluid 
characterized by = —1 in a dark energy-dominated Universe w = —1. The equation for \I' reduces 
to 
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In the three-dimensional momentum space, the solution to this equation is given by 

VI/ = ci {k') cosh (^^^ + C2 {k') sinh (^^^ + ^oa^/^-\ (6.29) 

where as before Ci{k^) are two integration constants. The density contrast is obtained by substituting 
this solution into eq. (6.25) 
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The density contrast as given by this relation should be the correct solution of our problem. However, 
things are more subtle here. For both matter and radiation perturbations, equation (6.17) gives 
once the potentials 4' and $ are determined. For dark energy perturbations, cx p {1 + w) vanishes 
identically so that equation (6.17) reduces to = —H^. Therefore, this is an extra equation which 
has to be satisfied by (6.29), and which implies that ci = C2 = ^'o = 0. Thus, the density contrast 
vanishes identically 5p = as in the conventional GR. It is worth noting that these results do also 
apply in an inflationary epoch. 

6.3 Summary and prospects 

Let us conclude this chapter by summarizing the previous discussion. After having discussed cosmo- 
logical models in massive gravity theories, we have studied the perturbations in one specific cosmo- 
logical model based on the sub-class of theories invariant under the dilatation symmetry (3.15). We 
have argued that cosmological perturbations contain two contributions, the "normal" and the "anoma- 
lous"one. The first, normal contribution has the behavior identical to that found in the conventional 
GR. It is therefore in agreement with observations to the same extent as the latter. In particular, the 
"normal"contribution can describe successfully at least the linear stage of the structure formation. 

The second, "anomalous"contribution of the cosmological perturbations is specific to the model 
of massive gravity with the dilatation symmetry (3.15). These perturbations originate from the con- 
densation of the four scalar fields present in the model and depend on the unknown function ^'o(a;*) 
already introduced in section 3.2. This function enters the solutions of the EoM as an integration 
constant. Hence, its value cannot be determined within the model (3.16) and has to be specified as an 
initial conditions. 

The behavior of the "anomalous"perturbations at different stages of the evolution of the Universe 
depends on the value of 7. At the matter-dominated stage the anomalous perturbations grow not 
faster than the standard ones for 7 < and 1 < 7. In the radiation-dominated epoch this occurs at 
— 1 < 7 < and 7 = 1. Thus, at — 1 < 7 < and 7 = 1 the normal perturbations dominate at both 
radiation and matter-dominated stages. During an inflationary stage or in a dark energy-dominated 
Universe, the "anomalous"perturbations cancels out in the density contrast. 

The appearance of this time-independent arbitrary function is not surprising since the same func- 
tion ^'o already enters the expression for the gravitational potential of an isolated massive body (3. 17). 
Its origin may be traced back to the existence of the scalar mode with the dispersion relation uj'^ = 0. 
This mode is not dynamical in the models considered. However, the action (3.16) is a low-energy 
effective action, so one should expect corrections containing higher-derivative terms to be present. 
In general, these corrections make ^'o a dynamical variable with the dispersion relation = ak^, 
where a is a small coefficient (see section 2.2.6). Therefore, ^'o becomes a slowly varying function of 
time. The slow evolution may drive "^q to a particular value at the inflationary epoch and thus prepare 
the initial conditions for the radiation-dominated stage. If this initial value of ^'o is small, then the 
growth of the "anomalous"part of perturbations may become irrelevant and corresponding values of 
7 phenomenologically acceptable. This question deserves furthers studies. 

Finally, it is worth noting that in the Standard Model of Cosmology, perturbations had to grown 
since recombination by a factor greater than the one predicted by GR alone. Indeed, GR fails at getting 
the baryons perturbations anywhere near the amplitude required to generate the observed structures 
of the Universe. For this reason, it is generally assumed that dark matter perturbations start growing 
before the emission of the CMB [193]. Things are different in massive gravity models since for 
7<— lorO<7<lorl<7, the density contrast grows faster than in GR. It would be interesting 
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to see if this could eventually eliminate the need for dark matter in the formation of structures. 
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7 

Condxisions & J^inaC comments 



In the first part of this thesis, we have reviewed the problems and issues of theories with massive 
gravitons, before introducing a particular class of theories for which the Lorentz symmetry is sponta- 
neously broken by the vacuum expectation value of four scalar fields. This discussion had two aims. 
The first one was to illustrate the difficulties which arise when trying to give a mass to gravitons. The 
second was to show how theories with Lorentz symmetry breaking go through all these difficulties 
and describe consistent effective field theories which reduce to GR in the limit of vanishing graviton 
mass. In this first part, we have also introduced a minimal class of theories possessing only two mas- 
sive propagating modes with helicities ±2, and have shown that for the sub-class of models invariant 
under the dilatation symmetry (3.15) there is no modification of the Newtonian potential. 

In the second part of this thesis, we have discussed the original contributions to the study of this 
minimal class of massive gravity theories. Let us summarize here the conclusions of these works. 

It has been argued in chapter 4 that physical instantaneous interactions are present in this minimal 
class of massive gravity theories because of the spontaneous breaking of the Lorentz-invariance. This 
symmetry breaking allows for independent modifications of the dispersion relations of the tensor, 
vector and scalar sectors so that the instantaneous contributions to the potentials do not cancel in 
the graviton propagator, unlike in GR. A concrete example of instantaneous interaction is given for 
theories invariant under the dilatation symmetry (3.15): we demonstrated by direct calculations that a 
gravitational source localized in space is responsible for an instantaneous frequency shift of a distant 
light beam. This interaction originates in an anisotropic sound wave which propagates at an infinite 
velocity in the Goldstone fluid. It is supposed to be responsible for the violation of the black hole 
"no-hair" theorem which motivates the work summarized in chapter 5. 

In chapter 5, the static vacuum spherically symmetric solutions in massive gravity have been 
obtained both analytically and numerically. The solutions depend on two integration constants, instead 
of one in GR: the Schwarzschild radius and an additional "scalar charge''^. At zero value of S 
and positive the standard Schwarzschild black hole solutions are recovered. Depending on the 
parameters of the model and the signs of and S, the solutions may or may not have horizon. Those 
with the horizon describe modified black holes provided they are stable against small perturbations. 
In the analytically solvable example, the modified black hole solutions may have both attractive and 
repulsive (anti-gravitating) behavior at large distances. At intermediate distances the gravitational 
potential of a modified black hole may mimic the presence of dark matter. Modified black hole 
solutions are also found numerically in more realistic massive gravity models which are attractors of 
the cosmological evolution. 
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Finally, in chapter 6 we have discussed more subde parts of modern cosmology, namely the mech- 
anism of structure formation. After having introduced cosmological models in the minimal class of 
massive gravity theories, we have studied the perturbations in one specific cosmological model based 
on the sub-class of theories invariant under the dilatation symmetry (3.15). The Friedman equation 
in these models acquires an unconventional term due to the Lorentz-breaking condensates which has 
the equation of state w = —1/ (87) with 7 being a free parameter taking values outside of the range 
[0, 1/3]. Apart from the standard contributions, the perturbations above the Friedmann background 
contain an extra piece which is proportional to an arbitrary function ^'(x*) of the space coordinates. 
This function appears as an integration constant and corresponds to a non-propagating scalar mode 
which may, however, become dynamical with the account of the higher-derivative corrections. For 
— 1 < 7 < and 7 = 1 the "anomalous"perturbations grow slower than the standard ones and thus 
the model is compatible with observations. Whether the model is experimentally acceptable at other 
values of 7 depends on the value of the function at the beginning of the radiation-dominated 

epoch. 

In this thesis, we have addressed only a few questions related to massive gravitons and possible 
modifications of GR at large scales. Many of the interesting questions still remain open. Among all 
possible further investigations, let us summerize some which may seem promising in light of what 
have been learned in this thesis. 

• It is still an open question how objects (stars, black holes, . . . ) with a non-zero scalar charge S 
may be created. Since a non-zero scalar charge will imply a modification of the Schwarzschild 
metric, this question surely deserves further investigations. 

• Another open question concern the stability of the spherically symmetric solutions. These so- 
lutions describe modified black hole provided they are stable against small perturbations. It is 
therefore important to study this issue in more detail. 

• The presence of a physical instantaneous interaction should allow to look behind the horizon of 
a black hole or behind the cosmological horizon. It would therefore be interesting to understand 
how information is carried through these horizons. 

• It would be interesting to compare the predictions of massive gravity theories concerning cos- 
mological perturbations directly to observations, in order to determine which values of 7 and 
^'o correspond to the best fit to observations. 

The work presented in this thesis correspond to a minimal class of models, which does not claim 
to be the only viable candidates for theories with massive gravitons. Even if some aspects of these 
models have been studied in this thesis and elsewhere, an intensive theoretical as well as experimental 
work has still to be done before one is able to conclude with certainty about their ability to describe 
correctly the gravitational interaction. Finally, we would like to stress that if GR correctly describes 
the gravitational interaction from cosmological scales up to the Planck scale, the study of massive 
gravity models enables one to comprehend which deviations from the standard laws of gravity are 
possible and how we may constrain them. 
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Part III 



T/te gravitationaijitidtqiiations 



The action for the massive gravitational field reads 

<S = <Sgr + with <5<^ = y" d^xv^A^J^ (X, F^ VF^^) , 
where X, and W^'-' are given by 



X ' 



and where are four scalar fields. The variation of this action with respect to the metric gives the 
Einstein equations 

which are ten equations determining the metric. 7^,^ is the energy-momentum tensor of standard 
matter fields while t,^^ is the energy momentum tensor of the scalar fields, obtained by varying <S^ 
with respect to the metric 



A' 



(A.1) 



The variation of the action with respect to the scalar fields give four Goldstone equations which read 



5, 



'/3 



-99 



V 



dX x^ " 



(A.2) 



It is worth noting that the Goldstone equations implies the conservation of energy-momentum ten- 
sor of the scalar fields. This is easily demonstrated after the introduction of more compact notations. 
Let x""" be 
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Then, the energy-momentum tensor of the scalar fields reads 




while the Goldstone equations are given by 



= 



5^ (y^gO^apg^'^dA^) 



With these new notations, it is straightforward to show that 




dx 



where V/j, (V^) is the covariant derivative with respect to x^^ (x^u). Consequently, the Goldstone equa- 
tions and the conservation equations for the scalar fields are identical provided that d^cj)" ^ 0. This 
discussion enable to identify the number of independent equations of the massive gravitational field. 
If there is no usual matter field in the right hand side of the Einstein equations, then the Goldstone 
equations are not independent equations since they are a consequence of the fact that the Einstein 
tensor is divergence-less. On the contrary, if one consider space-time with matter field in it, then 
the Einstein equations implies the conservation of the total energy-momentum tensor T^i, + t^i, with 
consequences that the Goldstone equations are independent equations. In these situations, one has to 
consider either the conservation equations for the matter fields or the Goldstone equations along with 
the Einstein equations. 

A.l Energy-momentum tensor 

The energy-momentum tensor associated to a matter field Lagrangian C is given by the variation of 
its Lagrangian with respect to the metric 



The tensor defined by this last relation is a symmetric tensor of rank two. It contains all the information 
concerning the energy and momentum of the matter fields described by this Lagrangian. 

Lets v^^ with V/^v'^ = 1 be the velocity of an observer comoving with the fluid described by the 
energy-momentum tensor. Then, the two following tensors are projection tensor along and perpendic- 
ular to 
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respectively. The energy-momentum tensor of the fluid can be decomposed as 
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where p = T^^v^v" and p = T^uj'^'^/'i are respectively the energy and pressure densities measured 
by the comoving observer, q^^ = —Tjxv'^j'^'^ is the energy flux perpendicular to v'' and ir'^'^ = 
'^X'^j^^P'^ + j^^'^XiP"' /"^ is the anisotropic pressure tensor 

With these notations, it is possible to determine the energy and pressure densities of the Goldstone 
scalar fields. Indeed, one has that 
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A.2 The mass parameters 

In backgrounds characterized by flat three-dimensional space gij oc 6ij, there are five mass parameters 
m? z = . . . 4, defined by the following relations 
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where W = —l/3SijW^^ and where the first and second non-zero derivatives of the function JT are 
denoted as follows 
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For function independent of as the ones considered in this thesis, the mass mf is proportional to a 
combination of the energy-momentum tensor of the scalar fields 

mj = -2 {p^ + p^) , 

with consequence that this mass vanishes in Minkowski space-time. For the class of models charac- 
terized by a function T = T (Z*-^) with Z*-^ = X'^W^^ , it is straightforward to show that the five 
mass parameters are 

-7[3(1 - 2-i)ZJ=z + 67^' (3J^zzi + 2J^zz2)], 
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where 



jr = J^Z^ij, 



For those particular models, one has 

mQ = 37 {m\ — m\/2) , m\ = 7 (3m| — m|) . 



A.3 Linearized gravity 

Since metric theories of gravity are generally non-linear, perturbation theory plays an important role 
in exploring such theories. Consider small fluctuations about a fixed background with metric 7^,^ 

With this parametrization, and with 7 = det (jnu) and h = h^^^^^^, the inverse of the metric and its 
determinant g = det {Qfiv) are given by 

- 1 1 

-h + - [ — 
2 4 V 2 

With these relations, the Ricci tensor is given up to quadratic order in the fluctuations by 

1 



'-9 



'-1 



1. 



where (V*) is the covariant derivative with respect to x'^ (x^). These relations allows to ex- 
press the Einstein-Hilbert action and the gravitational field equations up to the second order in the 
perturbations. 
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A.3.1 The gravity sector 

If the Ricci tensor is decomposed according to 

TZ^u + ^Tlnv + 5'^'Tl^lv + • ■ ■ , 

where STZ^y and are respectively the linearized and quadratic contributions to the Ricci tensor, 

then the Einstein-Hilbert action describing gravitational perturbations above a fixed space-time is 
given up to the quadratic order by the following relation 



5eh = -M^i / d^x V=7 



In this action, the zero order contribution fixes the background space-time. The first order contribu- 
tions are proportional to the EoM of the background and therefore vanish. The only relevant terms are 
the quadratic ones, which in Minkowski space-time reduce to 



EH 



_ 1 
Mink. 4 



There is two equivalent ways to get the linearized Einstein equations for the gravitational field. 
The first possibility is to linearize the action up to the second order, then varying this action with 
respect to the perturbations. This approach is used in the chapter 2 for Minkowski background. The 
second possibility is to linearize directly the Einstein equations. This second option is often easier 
since linear equations are simpler than quadratic actions. For this reason, this second procedure is 
used to derive the linearized EoM in FLRW background needed in chapter 6. 

The left-hand-side of Einstein's equations is given by the Einstein tensor 

In order to linearize this tensor, let us concentrate on the flat Friedmann-Lemaitre-Robertson-Walker 
solution described by the metric (6.1) with /c = 0. Then, with the notations introduced in section 6.2, 
the first order Einstein tensor reads 

SGio = d,[2{n^ + ^') + {2n' + 71^)3] + ]^d^^wi + s.,{2n' + re), 
5Gij = 5,j[2n{2i) + ^)' + 2iij" -dl{^ -^) + 2{ip + !f){n'^ + 2n')] 

+didj - $ - 2 {2n' + n^) E] - HHij + J {dlHij - h^j) + {2n' + n'^) Hij 



2 



where Tl is the conformal Hubble parameter and the prime denotes the derivative with respect to the 
conformal time t]. The Einstein tensor in Minkowski space-time is deduced from these relations by 
taking a = l,Tl = 0, iJj = and by replacing the prime by a dot which correspond to the derivative 
with respect to the time t. 
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A.3.2 The scalar fields sector 

As discussed in section 6.1, the homogeneous and isotropic ansatz for the Goldstone scalar fields 
reads 



(fp = (j) (rj) , 0* = X*. 



With this ansatz, it is possible to show that the Goldstone equation (6.4) can be put in the following 
form 



n = mmi 
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Hence, this equation relates the masses ttiq and m\. This relation will be used to express the linearized 
field equations in a gauge invariant formulation. But before writing the EoM, one need to determine 
the linearized energy-momentum tensor of the Goldstone fields. This tensor is given by 
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A.3.3 The Einstein equations 

The energy-momentum tensor of usual matter fields have been discussed in section 6.2 (see relations 
6.12). Therefore, with the relations introduced above one can write the Einstein equations 
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-didj - $ - 2 {2n' + e] - nH[^ + - {dlHij - h^j) + {2n' + h^) h, 
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Those equations correspond respectively to the 00, Oi and ij components of the Einstein equations. 
By making use of the Friedmann and Raychaudhuri equations (6.3) for the background along with the 
Goldstone equation in the form (A. 3), and with the gauge-invariant fields introduced in section 6.2, 
the linearized Einstein equations for the massive gravitational field consist of one tensor equation 



= Ff, - dfH^j + 2nH[, + I mlHi, + 



27r,; 



M 



pi 



two vector equations 



pi \ pi / 







w[ + 2Tlwi — I rn^cTi + 



2^ 



and four scalar equations 
= 25|^ - 6W + 



= 2di + n<^>) + a^di 



-0/ 



+ 



/ P<t> + P<l> 

V M^i 



+ ml {d]E + 3^) 



+ 



M 



pi 



Mi 



didj 





M 



pi. 



-2^'" - 2$ (W^ + 2n') + d] - - 2?^ (2^' + + 



2P(5p - dflT 



+a'^ml [ ^ 



^0/ 

A/ 



+ {ml - 3ml) ^■ 



91 



A.3.4 The Goldstone equations 

The variation of the action with respect to the scalar fields give four Goldstone equations (A.2). At 
the linearized level, those four equations give one vector equation 







dr, 



mf + 2 



P4,+P 



[ZD, 



and two scalar equations 

,4 







mo I ^ 



m, 



(3^ + df 



■=0 



di \ a ^do 



V 2 Kk 



+ d 



4 P^^/ 



+ (mi-mi)9|H 



A.4 Static spherically symmetric ansatz 

For the ansatz (5.1) discussed in section 5.2, the non-zero components of the Einstein tensor are given 
by the following relations 



1 



1 



r 

a + ra' 



1 

4r 



a' + ra" 



aj3 



+ 



2a + ra' 

a(3 



For functions JT which are invariant under rotations of the Goldstone fields 0* internal space, the 
derivatives of with respect to W^^ are given by 



where Ti = dT jdwi. Therefore, the components of the energy-momentum tensor (A.l) of the four 
Goldstone fields which are not identically zero are given by 



A^ 



2 a \^ dW'J X2 

I f dr y 1 dr 



A 



In 



2^ dw^j 



dr fV'V^ dr4'dr^ 

rx + 7^77777 I ..n + 



dW'i \ A2 
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Instantaneous interaction 



This appendix provided some details about the calculations presented in section 4 devoted to a physical 
instantaneous interactions in massive gravity models. 



B.l Transverse and longitudinal projections operators 

The key of the understanding of the instantaneous interaction present in Lorentz-violating electrody- 
namics and in Lorentz-violating massive gravity models lies in the concept of transverse and longitu- 
dinal sources. Defining the transverse projection operator by [156] 



5ij5^ (x — x') 



1 



1 



4tt dx'^dx'i \ \x — X 



and the longitudinal projections operator by 



Lij (x, X ) 



1 



1 52 



47r dx^dx'^ \ \x — x 



(B.l) 



(B.2) 



one can show that 



^ Tij (x,x') = -J^Tij (x,x') = ^t^jLki (x,x') = 0. 



d: 

The transverse and longitudinal components of any three-dimensional vector field (x) are given by 

(x) = j A^x'Tij (x, x') Vj (x) , (x) = j d^x'Lij (x, x') Vj (x) . 

One can easily verify the following identities 

V;'ix) + V/^ix)=Vix), -l-V^(x) = ^, (Vxy^)Jx) = 0. 



The longitudinal vector field can be expressed through 



Vl" = div{x) with V (x) 



1 

47r 



d'^x 



dx'^ V |x — x' 
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B.2 First order geodesies 

Let us first consider tlie geodesic motion of an observer of four- velocity f This observer follows a 
time-like geodesic x"^ (rob) = {t (Tob) , {'''oh)) with an affine parameter Tob- It obeys the following 
equation 

— + T„J'v^'v'' = with = — . 

dTob dTob 

The four velocity of the observer is such that v'^v^ = 1. A natural solution to the previous equation is 
Minkowski space-time is given by v'* = (1,0,0,0). Therefore, the geodesic followed by the observer 
reads {Tob) = {Tob,x''). 

Perturbations of the metric rj^i, + h^y imply perturbations of the geodesic {Tob) + ^x'^ {^ob) 
which in turn impUes perturbations of the four velocity v'^ {rob) + Sv"" {Tob). Since the observer is 
supposed to be massive, its perturbed time-Uke geodesic satisfies 

{v'' + Svn{Vt^. + Ku){v'' + Svn = l ^ <5^;° = -(^, 

implying that the zero component of the velocity perturbation is given by hoo/2. The linearized 
geodesic equation for the massive observer becomes 

ddv°' 

Under a infinitesimal gauge transformation x''^ = — of the coordinates, the geodesic 

perturbations transform as 

— > 6x'>' = Sxf" - y". 
Hence, the following quantities are gauge-invariant fields 

Sn^ = Sx'^ + B + E, Sn^ = Sx' - + d^E. 

It is then straightforward to show that the gauge-invariant four- velocity of the observer is fixed by the 
following equation 

with 

5v° = 5v^ + B + E, 5v' = Sv' -F' + d'E. 



Light follows null geodesies {Tph) = {t {Tph) , x* {Tph)) characterized by an affine parameter 
Tph and satisfying the following geodesic equation 

di/" dr" 

^Tph dTph 

with u^^Ufi = 0. A solution to this equation in flat space-time can always be formulated as = 
(l, ra*) with ra| = 1. Hence, the geodesic can be parameterized as x" {Tph.) = c" -I- Tphu" where c" is 
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a constant four vector. As previously, perturbations of the metric implies perturbations of the geodesic 
a^" (Tph) + ^x'^ {Tph) which in turn imply perturbations of the four-vector tangent to the geodesic of 
the photons ii" {Tph) + Su" {rph) ■ Therefore, one has 

{u^ + 5u^) {r}^„ + V) (^"^ + ^u") = ^ 2Su° + /iqo - ^n'Su' + 2n'hoi + n^n'hij = 0. 

This relation impUes that SvP is determined by (5u* and the metric perturbations. Since Su°' transforms 
as Sv"" under a four dimensional coordinates reparametrization and since 



dr, 



ph 



the gauge-invariant four- velocity of Ught reads 

5u^ = 5u^ + B + E + n'di [b + , 6u' = Su' - F' + d'E + n^dj {F' + d'E) . 
Hence, the geodesic of light is given to the first order by 



dr, 



ph 



n n- 



$ + * + 2n'5i$ + ( diWj + djWi - H, 



^3 



- \5^^ [Hij 



+ diWj — djWi 



{6rdj + Sfdi) ^ - -.5"^ {djHu + diHji - diHij) 



B.3 Frequency shift 

If an observer of four- velocity sees photons with a four-momentum = ujqu'^, he will conclude 
that they have a frequency fixed by the following relation 

Perturbations of the metric imply a Doppler frequency shift because of the previous relation. Indeed, 
the frequence of the Ught beam is given at the linearized level by 



00 



After substituting in this relation the solutions to the geodesic equations 



-2$ + / dTph 



n n-' 



(Hij 



diWj — djWi 



one finally found that the spectral shift is given by 



5uj 

CO 



dr, 



ph 



^' - n'9i$ + n'-n^ ( -Hij - diWj 



1 



+ / dTobu'di'^. 
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Substituting the solution of the linearized EoM of GR into this relation gives eq. (4.8) after some 
manipulations. 

In Lorentz-violating massive gravity models based on the action (3.1 6), there is an extra contribu- 
tion as compared to GR due to the presence of a mass term in eq. (3.3) for the gravitational waves. For 
simplicity, we will determined the derivative of the spectral shift with respect to the affine parameters 
Tph of the hght geodesic 

d fSuj 5uj^^\ 
dTph \uj w y ' 

and demonstrate that this term give rise to a physical instantaneous interaction. Since 

= ^ H^j with Hj^j = Hij — -ffijloR, 
one has to solve the two following equations 

(□+mi)ij,, = -^, □ij,,.|GR = -^. 



Sources of tensorial modes 

The transverse (B.l) and longitudinal (B.2) projections operators can be used to determined the trans- 
verse and traceless anisotropic stress tensor vTjj. Indeed, one can show that for an arbitrary energy- 
momentum tensor 6T^,^, the transverse and traceless anisotropic stress tensor is given by 



djdi 



dtdk Sij_ dkch 1 didk djdi 
dl +2 31^2 dl dl 



1 



djdj 



with for an arbitrary functions / 

1 



dl 



fix) 



1 

47r 



\x — X 



For the example of section 4.3.3 consisting in a source which appears at t = and described by 



5%o 

5%, 



4fiijtQ (t) dj6^ {x - Xs) , 
4fj.ije{t)5^ {x-xs), 



with Xs = (0, 0, d) and |Ujj = 62iS2j, the only component of the transverse and traceless anisotropic 
stress tensor Hij which enters the equation of the frequency shift through the gravitational waves Hij 
is given by 



7TU = 29 (t) 
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dfdj df 



(x - Xs) 



because the light's geodesic is supposed to be along the x direction with 



A- 
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Derivative of tiie spectral sliift 

The derivative of the spectral shift with respect to the affine parameters Tph of the light geodesic 
is given by the differences between the solutions of the tensor equations in GR and in the class of 
massive gravity models considered here. For the example of section 4.3.3 summarized above, on has 
that 

1 /d^d^ + 
1 fd'^d'^ a2\ ^ 



implying that 



Note that if x = t'^ — {x — Xg)'^, the difference AG+ between the Green functions of the d'Alembert 
and Klein-Gordon equations is given by 



AG+ {t, 0, X, X,) ^ G+ (t, 0, X, xs) - G+ (t, 0, x, x,) = -9 (x) ""'^^ ^"^^-^ 



where ^1 is a Bessel functions of the first kind. For m^x very small, the following Taylor expansion 
gives a good approximation of 



Ji {m2^/x) _ 1 



,2, 



m2^yx 2 16 



""2^ + O {ml) . (B.3) 



After some manipulations, the derivative of the spectral shift with respect to the affine parameters 
Tph of the light geodesic reads 

/ d'.'dV '"7''fa,- '''""'-f''fa;-'^'' AG-^ (*,0,."..,) . (B4) 

f iJ <-^' M-* 



6Mi 



The first term in the right-hand-side of this equation vanishes outside of the light cone of the source 
since it involves the Heaviside function (i2 — d"^) . Let us then concentrate on the second term. 
With the Taylor expansion (B.3), it is possible to show that the second contribution to a^j can be 
approximate at the origin x = by 



m2(i 
T28 



r+oc r+oo 

/ dr e (t^ -r^ -(f - 2rd) - dr & {t"^ - - cf + 2rd 
Jo Jo 



+ O (ml) . 



For < t < d which is outside of the light cone of the source, this term clearly becomes an instanta- 
neous contribution to the derivative of the spectral shift a^j 

^^^t + Oiml). 
64 ^ ^' 
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Since the third term in (B.4) is of order O (ml) , the conclusion to this calculation is that the derivative 
of the spectral shift with respect to the affine parameters r^^ of the light geodesic 



is instantaneously affected by the appearance of the source ST^j^y characterized by a non-trivial tenso- 
rial anisotropic stress. 
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